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SYNTHESIS  OF  OSCILLATING  ADAPTIVE  SYSTEMS 


Abstract 

Oscillating  adaptive  systems  are  non-linear  feedback  systems  in 
which  there  circulates  in  the  feedback  loop  (or  loops) , a signal  whose 
frequency  components  are  large  relative  to  the  bandwidths  of  the  control 
(command  and  disturbance)  signal  components.  Under  certain  conditions 
the  system  response  to  these  control  components  is  quasi-linear , 
permitting  linear  analysis  accurate  within  stringent  engineering 
tolerances.  The  system  response  to  these  control  inputs  can  also  then 
have  zero  sensitivity  to  the  (linear  time-invariant)  plant  high-frequency 
gain  factor.  If  a linear  time-invariant  feedback  system  is  designed  to 
cope  with  such  plants,  the  uncertainty  in  this  gain  factor,  if  very 
large,  can  lead  to  enormous,  possibly  impractical,  loop  bandwidths. 

Hence,  the  zero  sensitivity  property  of  the  oscillating  system  makes 
them  very  attractive  as  adaptive  systems. 

The  quasi-linearity  conditions  permit  the  development  of  quanti- 
tative frequency-response  synthesis  techniques,  whereby  one  may  proceed 
from  a statement  of  quantitative  performance  specifications  and  plant 
uncertainty  bounds,  step  by  step  towards  an  optimum  design.  One  signi- 
ficant result  is  that  the  gain  factor  uncertainty,  inherently  banished 
as  a direct  adaptive  problem  by  the  oscillating  condition,  reappears  in 
different  forms  as  a significant  design  constraint.  Six  oscillating 
adaptive  structures,  four  of  them  new,  are  presented  each  with  its  own 
quantitative  design  procedure.  Each  has  its  own  package  of  feedback 
benefits,  costs  and  trade-offs,  with  the  gain  factor  uncertainty 
appearing  in  different  forms.  Numerical  design  examples  are  presented 
for  each  structure,  including  verification  by  computer  simulation. 
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CHAPTER  ONE 
PRELIMINARY  BACKGROUND 

1.1  Introduction 

One  of  the  most  important  reasons  for  using  feedback  is  to  reduce 
system  sensitivity  to  uncertain  parameters  of  a constrained  part  (plant) . 

In  the  theoretical  development,  the  plant  is  assumed  linear  time-invariant, 
but  with  uncertain  parameters.  Rate  of  parameters  variation  and  ability 
of  the  system  to  adjust  to  fast  rates  of  variation,  are  discussed  in 
Chapters  4-7,  where  pertinent.  If  linear  time-invariant  feedback  com- 
pensation is  used,  giving  a totally  linear  time- invariant  feedback  system, 
it  has  been  shown  ' ' that  the  principal  "cost  of  feedback"  is  in 
the  bandwidth  of  the  loop  transmission,  which  may  make  the  system  very 
sensitive  to  sensor  noise  and  high  frequency  parasitics.  The  principal 
factor  causing  the  loop  transmission  bandwidth  to  be  large,  is  the  un- 
certainty in  the  high  frequency  gain  factor  K of  the  plant  transfer 

function  P(s)  = KP,  (s)  , where  lim  P,  = s P , p being  the  excess  of 

n 9-+"0  n 

plant  poles  over  zeros.  The  above  factors  have  been  discussed  in 
considerable  detail  in  the  literature  (4, 5, 6, 7]  and  are  therefore  not 
further  elaborated  here. 

One  of  the  important  motivations  for  non-linear  (so-called  adaptive) 

systems,  is  to  reduce  this  "cost  of  linear  feedback",  i.e.,  to  obtain 

the  same  benefits  with  effectively  less  loop  transmission  bandwidth. 

However,  with  the  exception  of  the  Non-linear  Oscillating  Adaptive 
18  9 101 

Systems  ' ' , there  has  hardly  appeared  any  quantitative  synthesis 

theory  for  these  adaptive  systems,  i.e.,  wherein  one  can  design  to 
achieve  assigned  quantitative  performance  bounds  over  a given  quanti- 


tative range  of  plant  parameter  values. 
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The  essential  feature  of  an  oscillating  adaptive  system  is  the 
introduction  of  a non-linear  element  and  a signal  whose  frequency  compo- 
nents are  sufficiently  large  relative  to  those  of  the  control  components. 
This  high  frequency  signal  may  be  due  to  self-oscillation  or  it  may  be 
externally  applied.  There  are  obvious  practical  reasons  for  it  to  be  a 
pure  sinusoidal  (so  assumed  in  this  work) , although  theoretically  even 
aperiodic  signals  may  be  used  . One  vital  element  in  the  system 
operation  is  its  quasi-linearity  i.e.,  under  certain  "quasi-linearity" 
conditions,  the  non-linear  system  behaves  effectively  like  a linear 
system  towards  the  "control"  signal  components  of  major  interest,  which 
are  those  due  to  the  useful  "command"  inputs  and  to  the  unwanted 
"disturbance"  inputs.  Quasi-linearity  permits  one  to  characterize  these 
low-frequency  components  as  distinct  and  separate  from  the  "high-frequency 
carrier"  signal  in  the  system  - just  as  if  linear  superposition  is 
applicable. 

The  errors  involved  in  assuming  such  linear  superposition  are  finite 

r 12 1 

but  can  be  made  as  small  as  desired  by  designing  the  system  that 

at  the  input  to  the  non-linearity,  the  low-frequency  control  components 
are  commesurately  sufficiently  "small  and  slow"  relative  to  the  high- 
frequency  signal.  This  quasi-linearity  phenomena  has  long  been 
known  t 13 ' 14 ' 15 ' and  has  been  studied  extensively,  with  some 

interesting  results  recently  obtained  with  functional  analysis 

. u • till 
techniques 

This  present  research  concentrates  on  practical  engineering 
structures  and  design  techniques  for  the  adaptive  problem,  so  that  in 
view  of  the  extensive  literature  available,  the  quasi-linearity  aspect 
of  the  problem  is  not  treated  in  detail  in  this  work,  but  the  results 


I 


are  summarized  as  needed.  From  the  point  of  •. iew  of  this  work  the 
important  point  is  that  as  previously  noted,  the  error  involved  in  using 
superposition  is  finite  but  can  be  made  as  small  as  desired.  It  will  be 
noted  in  subsequent  chapters  that  it  is  simply  a matter  of  making  certain 
design  parameters  sufficiently  large. 

The  second  vital  element  in  the  operation  of  the  oscillating 

adaptive  system  is  its  zero  sensitivity  to  the  plant  high-frequency  gain 

factor.  This  is  related  to  the  quasi-linear  phenomenon.  The  latter 

permits  the  characterization  of  the  non-linear  element  N by  two  (or 

if  necessary,  more)  "describing  functions",  one  N for  the  high- 

o 

frequency  component,  the  second  , for  the  control  components.  Des- 
cribing function  representation  is  very  accurate  for  , for  it  is 

easy  and  otherwise  also  desirable,  to  design  the  feedback  loop  to  be  very 

highly  "low-pass".  Also  the  non-linear  element  needs  not  be  particularly 

1 171 

rich  in  higher  harmonics,  and  it  has  been  shown  that  under  these 

conditions  the  Nq  characterization  can  be  extremely  accurate.  The  use 
of  is  justified  by  the  quasi-linear  conditions  previously  discussed. 

Now,  turns  out  to  be  inversely  proportional  to  the  amplitude  of  the 

high-frequency  component,  which  in  turn  is  shown  later  to  be  directly 
proportional  to  the  high-frequency  plant  gain  factor  K , so  that  their 
product  KNf  is  independent  of  K 

History 

The  quasi-linear  feature  of  the  system  under  consideration  has  been 
well  known  and  practically  exploited  for  some  time  t 13 ' 14 ' 15 ' in  order 
to  "linearize"  unavoidable  non-linear  elements  with  undesirable  features, 
sometimes  present  in  practical  control  systems.  However,  the  zero 
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sensitivity  property  was  recognized  only  comparatively  recently, 

r 18  19] 

apparently  first  by  Minneapolis-Honeywell  workers  1 ' . They  promoted 

this  approach  in  the  competition  among  various  aerospace  companies  in 

the  late  1950 's  and  early  1960 's,  to  apply  adaptive  methods  to  the  flight 

control  problem.  In  fact,  the  oscillating  adaptive  approach  was  probably 

the  most  successful  of  these,  and  was  applied  to  the  X-15,  X-20 
[20  21  22] 

vehicles  ' ' among  others.  However,  Minneapolis-Honeywell  workers , 

although  they  reported  frequently  in  a qualitative  manner  on  their  work, 
never  did  produce  any  quantitative  design  theory.  It  seems  quite  certain 
that  their  design  was  by  ad-hoc,  engineering  cut  and  try,  because  to 
date,  long  after  they  have  ceased  promoting  these  systems,  they  have  not 
produced  any  such  design  theory.  It  is  interesting  that  the  final 
Minneapolis-Honeywell  design  is  a two-loop  system  with  an  objective 
similar  to  that  sought  in  the  two-loop  system  of  Chapter  4.  It  is 
superior  in  that  sense  (as  will  be  seen)  to  the  single-loop  systems  of 
Chapters  2-3.  However,  it  will  be  shown  that  there  are  important  advantages 
in  the  two-loop  system  of  Chapter  4,  over  that  used  by  Minneapolis- 
Honeywell. 

[22  23] 

Gelb  and  Vander  Velde  ' were  apparently  the  first  to  do  some 

analytic  work  on  the  adaptive  properties  of  oscillating  adaptive  systems. 

However,  their  major  attention  was  focused  on  the  quasi-linear  properties, 

and  no  attempt  was  made  to  develop  a quantitative  design  procedure. 

Consequently,  one  vital  design  parameter  (later  noted)  was  completely 

ignored,  without  which  a quantitative  design  theory  is  impossible  to 

( 8 1 

establish.  Horowitz  was  the  first  to  develop  a quantitative  design 
procedure  for  one  (of  the  six)  oscillating  adaptive  structures,  but  only 
for  one  of  the  two  possible  design  categories  (discussed  later)  for 


it* 
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this  one  structure.  This  category  is  restricted  to  designs  whose  loop 
transmissions  are  smooth  in  a certain  sense,  a constraint  which  is 
dropped  in  the  present  work. 

It  will  be  seen  that  each  of  the  six  structures  has  its  own  package 
of  trade-offs  between  feedback  costs  and  benefits,  which  offers  important 
flexibility  to  the  designer,  for  he  can  choose  that  package  which  is 
best  suited  to  his  particular  problem  package.  Some  of  the  important 
factors  which  appear  in  these  trade-off  packages  are: 

Kmax 

1.  The  amount  of  high-frequency  gain  factor  uncertainty  — 

min 

2.  The  ratio  of  largest  (low-frequency)  control  signal  input  which  the 
system  is  to  process  quasi-linearly , to  the  smallest  signal  parameter, 
which  is  the  maximum  tolerable  amplitude  at  the  plant,  of  the 
circulating  high-frequency  component. 

3.  The  maximum  rate  of  variation  of  the  high-frequency  gain  factor  K . 

4.  The  presence  of  higher  order  modes,  such  as  elastic  bending  modes 
in  an  air-frame. 

The  work  is  divided  into  six  main  parts: 

a.  SOAS  - Single  loop  Self-oscillating  Adaptive  Systems.  Chapter  2. 

b.  EEAS  - Single  loop  Externally  Excited  Adaptive  Systems.  Chapter  3. 

c.  SOAL  - A two  loop  SOAS  with  adaptive  saturating  level.  Chapter  4. 

d.  SOANL  - SOAL  with  non-linear  compensation  in  the  secondary  loop. 

Chapter  5. 

«•  EEAL  - A two  loop  EEAS  with  adaptive  saturating  level.  Chapter  6. 
f.  EEANL  - EEAL  with  non-linear  compensation  in  the  secondary  loop. 
Chapter  7. 


1 


The  first  two  systems  (SOAS  and  EEAS)  involve  known  single-loop 


structures  , but  with  hitherto  unknown  or 


theory.  The  last  four  systems  (SOAL,  SOANL,  EEAL,  EEANL)  involve  new 


(Self  Oscillating  Adaptive  Systems) 

2.1  General 

An  important  motivation  for  using  nonlinear  compensation  (i.e.  an 
adaptive  system)  is  to  reduce  the  loop  transmission  bandwidth,  and  still 

[o] 

get  the  desired  benefits  of  feedback.  The  SOAS  is  one  of  the 
simplest  Quasilinear  Adaptive  Systems  and  its  structure  is  shown  in 
Figure  2.1-1.  By  quasilinear,  is  meant  here  a nonlinear  system  which 
acts  very  closely  like  a linear  system  towards  the  control  signals  of 
interest,  if  some  conditions  (quasilinear  conditions)  are  met.  So,  a 
vital  prerequisite  for  the  SOAS  and  for  all  the  other  structures 
presented  here  is  non-violation  of  quasilinear  conditions  as  expressed 


in  Section  2.2. 
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G^(s)  , G2  (s)  and  G^fs)  , are  the  transfer  functions  of  the 
linear,  time-invariant  compensations  that  have  to  be  found.  P(s)  is 
the  constrained  uncertain  plant  transfer  function.  N is  the  introduced 
non  linearity.  The  quasi-linearity  conditions  apply  to  its  input. 

r(t)  is  the  applied  command  input  to  the  system  with  Laplace 

transform  R(s)  ■ J£(r(t)]  . It  is  important  to  define  the  "extreme" 

expected  command  input  r (t) 

© 

d (t)  represents  the  disturbance  referred  to  the  plant  output.  In 

case  the  real  disturbance  is  applied  at  a different  point  in  the  plant, 

the  equivalent  disturbance  at  the  plant  output  should  be  calculated, 

with  .its  Laplace  transform  D(s)  = £[d(t)]  . As  with  the  command  input, 

it  is  important  to  define  its  "extreme"  value  d (t) 

e 

n (t)  is  the  effective  noise  due  to  sensor  or  other  elements.  It 
is  assumed  its  characteristics  are  stationary,  gaussian  with  known 
power  spectrum. 

x(t)  is  the  input  to  the  nonlinear  element  N . If  quasilinear 
conditions  are  satisfied,  x(t)  has  four  distinct  components. 

x(t)  = x (t)  + x (t)  + x (t)  + x (t)  (2.1-1) 

o r d ti 

The  nonlinear  element  N is  chosen  so  there  is  self-oscillation  with 
frequency  in  the  loop.  If  the  other  signals  satisfy  quasi  linear- 

ity conditions,  there  is  no  possibility  for  "quenching"  this  self- 
oscillation. The  sub  o refers  to  the  oscillatory  component  of  x(t) 
with  frequency  . The  sub  r refers  to  the  component  of  x(t)  due 

to  the  applied  command  input  r(t)  . The  sub  d refers  to  the 
component  of  x(t)  due  to  d(t)  and  the  sub  n to  that  due  to  the 
noise  n(t)  . In  this  work  the  objective  is  to  find  a design  procedure 
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for  satisfying  system  specifications  while  minimizing  the  effect  of  the 
noise  at  the  plant  input  as  defined  in  Appendix  I. 

y(t)  is  the  non-linearity  output  signal.  Under  quasi-linearity 
conditions  y(t)  is  also  composed  of  four  distinct  components. 

y(t)  = y (t)  + y (t)  + y (t)  + y (t)  (2.1-2) 

o r a rj 

where  the  subs  have  the  same  meaning  as  in  (2.1-1). 

w(t)  is  the  plant  input  signal.  Again  under  quasi-linearity 
conditions,  superposition  is  valid,  giving 

w(t)  = w (t)  + w (t)  + w . (t)  + w (t)  (2.1-3) 

o r d n 

The  synthesis  technique  will  minimize  w (t)  in  the  sense  described  in 

n 

Appendix  I . 

z(t)  is  the  plant  output.  If  2.1-3  holds,  then 

z(t)  = z (t)  + z (t)  + z (t)  + z (t)  (2.1-4) 

o r d n 

Similarly  the  system  output  c(t)  is  then 

c(t)  = c (t)  + c (t)  + c _ ( t)  + c (t)  (2.1-5) 

o r d n 

Note  that  z (t)  = c (t)  , z (t)  = c (t)  , z (t)  = c (t)  , while 

o o r r n n 

cd(t)  = d (t)  + zd(t) 

The  validity  of  the  approximations  used  for  the  non-linear  element 
N , depends  on  "how  sinusoidal"  is  XQ(t)  around  its  zero  value. 
XQ(t),  as  explained,  is  the  oscillating  component  at  the  non-linearity 
input.  It  will  be  seen  that  there  is  no  problem  in  having  iG^G^fjcj)! 
decrease  very  rapidly  vs.  o>  , for  u)  > . Hence  the  higher 

harmonics  of  y^  (the  oscillating  component  of  y ) are  greatly 
attenuated  and  x (t)  is  very  well  approximated  by  a sinusoid,  in  a 


,r_ 


to 


proper  design. 


x (t)  = A sin  (w  t)  (2.1-6) 

o o 

For  the  above  reason,  describing  function  representation  of  N is 
highly  valid. 

2-2  Model  of  the  Non-linear  Element  and  the  Resulting  Constraints. ^ 12 * 
The  results  obtained  in  Reference  [12]  will  be  used  and  are  repeated 

here. 


Figure  2.2-1  - General  time-invariant  approximator 
for  a non-linear  operator 


The  model  for  the  non-linearity  .is  obtained  by  approximating  the 
non-linearity  operation  for  each  specific  signal  by  a corresponding 
linear  operation.  The  input  to  the  non-linearity  is  composed  of  four 
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different  signals  as  per  2.1-1.  For  each  of  these  inputs,  and  in  place 
of  the  common  non-linearity  N , a linear  operator  is  placed  as  shown 
in  Figure  2.2-1.  The  criterion  used  to  define  each  of  the  linear 
operators,  is  to  minimize  the  mean-square  difference  between  the  real 
non-linearity  output  y(t)  and  its  approximation  y(t) 

" 2 

(y(t)  - y(t))  > minimum  (2.2-1) 

The  signals  considered  are  stationary  ( x^(t)  is  stationary)  and  the 
non-linearities  considered  are  of  the  time  invariant  type. 


I 


! 


I 


I 

I I 

I 


By  proper  handling  of  the  equations  due  to  2.2-1,  the  result  is  a 
set  of  equations  on  the  weighting  functions  h^  , representing  the 
linear  approximations  in  Figure  2.2-1,  which  are  also  time  invariant. 


/ h.(t)R  (t-r)dt  = R (t)  ; Ti 

J r i x . x . yx . 


* 0 


(2.2-2) 


l l 


R is  the  autocorrelation  function  for  each  input  x.  , and  R 

x . x . i yx . 

li  l 

is  the  cross  correlation  function  between  the  non-linearity  output  y(t) 
and  each  input  x^(t)  . The  left  side  of  Equation  2.2-2  is  the  cross 

correlation  between  the  input  x^(t)  and  its  equivalent  approximation 


y^t) 


Since  the  several  inputs  x^(t)  are  uncorrelated,  the  left 


side  of  2.2-2  represents  the  cross  correlation  between  the  input  x^(t) 
and  the  total  approximate  output  y(t)  . So,  Eouation  2.2-2  tells  us 
that  if  the  cross-correlation  function  between  each  input  and  the  non- 
linearity output  equals  the  cross-correlation  function  between  the  same 
input  and  the  approximate  output,  then  the  mean  square  error  between  the 
approximate  ( y ( t ) ) and  real  (y(t))  output  is  minimized. 
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The  linear  operators  shown  in  Figure  2.2-1  can  be  derived  for  three 
types  of  signals.  The  first  is  a sinusoidal  signal  and  suits  the 


definition  of  XQ ( t)  . The  second  type  is  of  a slow  signal  if  compared 


to  xQ(t)  , and  by  a proper  design  of  the  system,  includes  x^(t)  and 


x,  (t)  . The  third  is  a random  signal  and  suits  the  definition  of 

a 


x (t)  , by  a proper  system  design.  Equation  2.2-2  can  be  repeatedly 


U2] 


used  to  define  the  several  describing  functions  for  these  linear 
operators : 


Non-linearity  model  for  the  sinusoidal  component  (representing  x^ (t) ) 


In  the  case  the  non-linearity  is  static  and  single  valued,  the 


approximator  for  the  non-linearity  Nq  , becomes  a real  value 


N = — y(t)  sin 
o A 


(2.2-3) 


Where  XQ(t)  = A sin  ip  , and  y(t)  is  the  non-linearity  output  due  to 


x (t)  . 

o 


Non-linearity  model  for  the  "Slow  component"  (representing  x (t)  and 


*d<t>± 


In  the  case  x (t)  and  x. (t)  are  "slow"  compared  to  x (t)  , 

r d o 


the  approximator  for  the  slow  components  is: 


N = ^- 


r xr(t) 


N . = 


d Xd(t) 


N = N = N.  = 
f r 


d xf(t) 


(2.2-4) 


Where  N^  and  N^  represent  the  approximator  for  the  signals  due  to 


the  command  ( x (t)  ) and  disturbance  ( x,  (t)  ) inputs,  and  N.  is  the 
r at 


approximator  for  the  "forcing"  signals  due  to  either  the  command  or 


disturbance  signals  ( xf(t)  ) 


i 
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Equation  2.2-4  has  a very  direct  physical  interpretation.  y(t) 
is  the  d.c.  value  of  the  non-linearity  output.  xf(t)  is  the  value  of 
the  "slow"  input  to  the  non-linearity.  The  non-linearity  approximator 
for  this  signal  is  like  a d.c.  gain,  being  the  relation  between  the  d.c. 
value  of  the  output  and  the  value  of  the  input.  To  avoid  infinite 
values  for  Nf  , when  xf(t)  tends  to  zero,  only  non-linearities  that 
have  no  constant  outputs  when  the  input  is  zero,  are  considered.  When 
xf(t)  -M3  , then  y(o)-»-0 

Non-linearity  model  for  the  random  signal. 

The  random  signal  at  the  non-linearity  input  is  due  to 

noise  in  the  system  output  sensor  n(t)  . In  the  case  x (t)  is  a 

n 

random  noise  with  gaussian  distribution  and  zero  mean,  then  the 

approximator  N for  the  random  input  x (t)  is: 

n n 

Nn  = "T  y(t)xn(t)  (2.2-5) 

o 

Where  o is  the  standard  deviation  of  x (t)  . Thus,  for  the  random 

n 

input,  the  non-linearity  approximator  is  described  as  a pure  gain,  when 
the  conditions  stated  above  are  satisfied. 

Non-linearity  model  for  combined  inputs. 

If  the  input  to  the  non-linearity  is  composed  of  a sinusoidal  and 
a slow  signal,  x(t)  = XQ(t)  +x^(t)  , particular  approximators  are  used 

yielding  the  Dual-input  describing  function  (D.I.D.F.).  The  D.I.D.F. 
"quasi-linearization”  of  the  non-linear  element  will  allow  analysis  and 
synthesis  of  systems  where  limit  cycling  exists  and  the  command  input 
and  disturbance  responses  are  the  final  issue.  We  can  assume  the  "slow- 
signal"  xf(t)  * x^Ct)  +xr(t)  , so  x ( t)  = Xq ( t)  +xd(t)  +xr(t)  . By 


"slow  signal"  we  assume  x^ft)  and  x^(t)  vary  "little"  during  a cycle 

of  oscillation  of  xQ(t)  . The  non-linear  approximation  is  called 

"quasi-linearization"  and  not  "linearization"  for  two  reasons:  first, 

each  of  the  h^(t)  in  Figure  2.2-1  performs  as  "a  linear  system"  to  the 

particular  set  of  inputs  x_^(t)  for  which  it  is  derived,  and  some  of 

the  h^ft)  may  depend  on  some  parameters  of  its  own  of  other  inputs 

(like  the  magnitude  of  the  oscillation  x (t)  ).  Second,  as  a 

o 

consequence  of  the  approximator  dependence  on  these  parameters,  there 
will  be  some  constraints  on  the  parameter  uncertainty  region  which 
allow  the  same  linear  model. 

The  non-linearity  used  in  this  work  is  an  ideal  relay.  The  same 

results  can  be  generalized  for  any  saturating,  static,  single  valued, 

odd  non-linearity,  if  it  is  driven  "hard"  enough  into  saturation,  twice 

in  each  cycle  of  the  oscillating  signal  x (t)  =A  sin  (u  t)  . Figure 

o o 

2.2-2  represents  the  ideal  relay  representation,  where  x(t)  is  its 
input  and  y(t)  its  output. 

iy 

+B 


X 

y 

* 0 

+B 

< 0 

-B 

Figure  2.2-2  - Ideal  relay  representation 


Figure  2.2-3  illustrates  the  ideal  relay  operation  (input-output) 


when  the  input  x(t)  =x  (t)  + x (t) 


Figure  2.2-3  - Ideal  relay  operation 
Input  - x(t)  = sinus + slow  signal 


non-linearity 
/ output 


the  switching  point  is  not 


symmetrically  related  to  the  oscillation,  so  a bias  component  exists 


To  calculate  the  non-linearity  approximator  N for  the 


sinusoidal  component.  Equation  2.2-3  is  used,  and  the  result  is 


To  calculate  the  non-linearity  approximator  N for  the  slow  component 


Equation  2.2-4  is  used  and  the  result  is 


„ _ _2_B_  -1  /V 1 

f TTX  Sln  V A / 


(2.2-7) 


f 1 

If  — $—  , Equations  2.2-6  and  2.2-7  can  be  approximated  by 

ft  J 


» -4| 

o ttA 


( 2 . 2-8a) 


■«-s 


(2.2-8b) 


With  an  error  of  less  than  6%  and  2%  respectively. 

More  than  that,  considering  xf(t)  not  a constant  value,  but  a 

relatively  slow  sinusoidal  signal  compared  to  XQ(t)  , and  by  using  the 

two  sinusoidal-input  describing  function  for  x (t)  and  x^(t)  , the 

o r 

results  obtained  in  Equation  2.2-8  are  still  valid  with  an  extra  error 

of  less  than  5%,  if  the  frequencies  of  the  slow  signal  to  satisfy 

Xf 

the  relation  [12] 


0) 

*>  * -f 

Xf  3 


(2.2-9) 


Hence,  the  idea  of  quasi-linear  approximation  as  stated  above  is 
quite  clear.  The  quasi-linear  approximations  for  both  the  sinusoidal 
component  ( Nq  ) and  the  bias  component  ( ) , are  functions  of  para- 

meters defining  the  inputs  themselves  (Equations  2.2-6,  2.2-7).  If  some 
constraints  are  included  in  these  signals  (quasi-linear  conditions) , the 
quasi-linear  approximator  becomes  simpler,  but  is  valid  only  in  the 
region  defined  by  the  constraints.  Summarizing,  the  quasi-linear 
conditions  are  (for  reasonable  engineering  values  of  a , 6 ) : 


1 


I 


Xr(t)  1 

< — 

A a 


a = 3 , 


— 6 = 3 ' 

'i)  6 

i o 


Xd(t)  < 1 
A a 


%d  < 1 
o 


(2.2-10) 


If  the  conditions  are  satisfied,  the  non-linear  approximators  (for 
an  ideal  relay)  for  the  several  signals  are: 


. 


(for  x (t)  ) 
o 


2r  Mf 

N = N = N = — — = — (for  x (t)  and  x,  (t)  ) 
f d r ttA  A r d 


(2.2-11) 


An  important  property  to  point  out:  If  the  system  includes  a saturating 


non-linearity  (specifically  an  ideal  relay) , and  there  exists  an 
oscillation  at  its  input,  the  transmission  for  the  relatively  slow 


signals  due  to  the  command  or  disturbance  inputs  is  defined  by  the 

oscillation  amplitude.  This  property  is  going  to  be  used  and  is 

responsible  for  the  adaptive  properties  of  the  non-linear  adaptive 

systems.  In  Equation  2.2-10  u>br  and  represent  the  bandwidths 

of  the  signals  x^(t)  and  x^(t)  respectively.  The  physical  meaning 

of  2.2-10  is  that  x (t)  and  x,  (t)  must  be  "small"  and  "slow"  when 

r a 

compared  respectively  to  the  amplitude  A and  to  the  frequency  u>o  of 

the  oscillation  x (t) 
o 

Suppose  now  the  input  to  the  non-linearity  is  composed  of  a 

"fast"  sinusoidal  component  x (t)  =A  sin(w  t)  , of  a "slow"  sinusoidal 

o o 

component  = b sinfu^t)  and  gaussian  noise  x^(t)  , with  standard 

r 121 

deviation  a and  mean  zero.  It  is  shown  that  if 


mJ. , i 

w R 
o 


b 1 
— < — 


o 1 


a = B = 3 


(2.2-12) 


-18- 


Then,  the  approximators  for  the  non-linearity  regarding  the  random  noise 

x (t)  , the  forced  signal  x.(t)  and  the  oscillating  component  x (t) 

f|  r o 

are  (for  an  ideal  relay) 


2B 
IT  A 


4B 

TTA 


(2.2-13) 


While  precise  accuracies  for  combinations  of  these  inputs  are  not  avail- 
able, good  engineering  results  are  expected.  In  any  case  if  there  is 
such  a combination,  the  appropriate  values  of  a and  8 can  be 
determined. 

Two  points  should  be  mentioned.  First  - even  if  n (t)  has  a 


I 


normal  amplitude  distribution,  the  output  of  the  saturating  non- 
linearity does  not  have  the  same  properties,  and  as  per  Figure  2.1-1, 
the  saturating  non-linearity  output  is  fedback,  so  the  actual  input  to 


t 


il 

* '! 


the  non-linearity  x^  (t)  is  expected  to  be  non-gaussian . But  G^G.^Pfjw) 

as  mentioned  before,  has  low  pass  characteristics,  which  tends  to 
[241 

"normalize"  the  amplitude  distribution  of  the  noise  component  x^ (t) 
at  the  non-linearity  input.  Second,  Equation  2.2-12c  will  not  affect 
the  design,  because  the  approach  is  to  satisfy  the  r(t)  and  d(t) 
specifications,  with  effectively  minimum  loop  bandwidth. 

2 . 3 Mathematical  relations 

Some  relations  involving  the  functions  and  the  signals  in  the  loop 
are  developed  in  this  section.  In  Figure  2.1-1,  let: 


L ( ju>)  A G G2P(j<o)  (2.3-la) 

Lf  ( ju>)  A NfL(ju>)  (2.3-lb) 


Where  the  gain  K is  explicitly  displayed,  then 
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From  Equation  2.3-lb  the  forced  loop  transmission,  which  is  the  one 
that  defines  the  system  response  for  both  the  command  and  disturbance 
input,  is 

or 

Lf  f joj)  = NfL  ( jw ) = Nf.G1G2P(jw)  = — 


K 2B 

Lf(jw)  = A ' ~ • GlG2Ph(ja)) 


Substituting  2.3-7  in  2.3-8, 


T . _ wi 

Lf(3“  2lGiG2Ph(jWo) 


(2.3-8) 


(2.3-9) 


which  reveals  the  outstanding  adaptive  property  of  the  SOAS  - Inherent 

zero  sensitivity  to  the  plant  parameter  K - the  plant  gain  factor. 

If  there  is  no  uncertainty  in  P,  (and  obviously  negligible  uncertainty 

n 

in  the  inserted  "soft"  compensation  function  G^(joj)  ),  then  L^(jw)  , 

which  is  the  effective  loop  transmission  for  the  command  and  disturbance 

signals,  has  no  uncertainty.  The  reason  is  that  when  the  plant  gain 

factor  K changes  its  value,  the  limit  cycle  amplitude  A changes 

£ 

proportionally,  so  — is  constant  in  Equation  2.3-8,  i.e.  the  variation 
in  K is  cancelled  out  by  the  proportional  variation  in  A . Note  the 
strange  combination  of  the  two  loop  transmissions. 

4B 

L (jw)  = — G, G-KP,  ( jw)  , is  valid  only  for  the  signal  due  to  the 
o ttA  1 2 h 

self-oscillation,  and  determines  the  amplitude  A of  the  oscillating 

U) 

2B  o 

component  x (t)  . L-(ja>)  = — G, G„KP,  ( ju>)  , is  valid  only  for  go  £ -5- r 
o f ttA  1 2 h 8 

0) 

o 

and  its  gain  is  determined  by  A . In  the  frequency  regions  — £ a)  < o)q, 
and  0)  > cdo  , the  loop  is  not  well  defined,  so  some  "strange"  results 
can  be  expected  if  we  try  to  extrapolate  results  from  the  allowed  bounds. 


a 

■-a 
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For  example,  if  a command  or  disturbance  input  has  strong  components  in 

the  frequency  range  near  , it  will  affect  xq = A sin(wot)  , by 

changing  its  amplitude  (depending  on  the  phase  of  the  signals  relative 

to  that  of  xq  ) , so  the  real  gain  identification  and  compensation 

A 

accomplished  by  having  — constant  (2.3-7)  is  spoiled.  Similarly,  in 

K. 

Figure  2.1-1,  the  plant  output  Z ( j<^0)  is  determined  by 

i , 4M  I I 

|Z(ju)o)  I = — - |G2P(jwQ)  | , and  is  independent  of  both  the  command  and 

disturbance  input  components  at  u) 


2 . 4 SPAS  sensitivity  to  uncertainty  in  plant  dynamics  P^(joi) 

If  some  of  the  dynamics  characteristics  of  the  plant  (zeros  and 
poles)  are  unknown,  it  will  affect  the  loop  transmission  and  the 
oscillation  frequency,  since  the  amplitude  and  the  frequency  of  the 
oscillation  must  satisfy 


4B 


‘ -1 


(2.4-1) 


In  this  research  we  consider  only  uncertainty  in  the  plant  gain 
factor  K of  P(ju)  =KPh(jw)  . A generalization  of  the  methods 
developed  here  is  possible  by  using  the  technique  of  References  9,10. 


2 . 5 Synthesis  procedure  for  the  SPAS 

In  this  section  a synthesis  procedure  is  derived  for  the  SPAS  to 
satisfy  a certain  specifications  set.  The  synthesis  procedure  enables 
the  designer  to  proceed  systematically  step  by  step,  and  obtain  a design 
which  satisfies  the  specifications  and  closely  minimizes  the  effect 
of  sensor  noise  at  the  plant  input,  in  the  sense  defined  in  Appendix  I. 


1 
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The  non-linearity  used  is  an  ideal  relay,  but  the  results  are  valid  for 
any  saturating,  static,  single  valued,  odd  non-linearity,  if  it  is  driven 
"hard"  enough  into  saturation. 

2.5.1  Specifications  and  data 

It  is  assumed  that  the  set  of  specifications  and  data  listed  below 
are  known . 


Desired  Command  input  response 


The  desired  response  to  the  command  input  is  supposed  to  be  known, 
in  a frequency  response  form,  either  analytically  or  graphically.  Since 
only  the  plant  gain  factor  uncertainty  is  considered  here,  there  is  no 
need  for  providing  bounds  on  the  permissible  variation  in  the  system 
transfer  function.  The  transmission  for  the  command  input  is 
(Figure  2.1-1) 


C 

Vju,)  = T(ja,) 


Lf ( jw) 

G3  1+Lf  (jui) 


(2.5-1) 


This  expression  is  valid  with  good  accuracy  for 
is  valid  for  this  range,  where  u>Q  5 , and 

the  desired  transmission. 


u 

a)  $ — , because  L,  ( j to) 

o t 

is  the  bandwidth  of 


Disturbance  Attenuation 


The  maximum  permitted  effect  at  the  system  output  due  to  the 
disturbance  signal  is  assumed  given  as  a function  of  frequency,  analytic- 
ally or  graphically.  The  disturbance  transfer  function  is  (Figure  2.1-1) 
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Since  L.(ju)  is  independent  of  K in  P = KP  , so  is  T (joo) 


Equation  2.5-2  is  well  valid  for  u>  $ — , so  if  there  is  a special 

p 

interest  in  defining  the  disturbance  attenuation  up  to  a certain 

frequency  u , , then  the  oscillating  frequency  should  be  a>  i Bw, 

a o d 


The  plant  definition 


The  representation  of  the  plant  by  its  differential  equations  or 
Laplace  transform  is  assumed  known,  as  well  as  the  range  of  its  gain 
uncertainty. 


P(jw)  = KPh(jw) 


(2.5-3) 


K1  * K * K2 


Extreme  command  input 


(2.5-4) 


In  a typical  "purely  linear  system"  this  is  not  defined  unless  there 
are  some  expected  non-linear  problems  like  velocity  or  acceleration 
saturation,  etc.  For  the  SOAS,  and  for  all  the  other  oscillating  systems, 
where  quasi-linearity  is  a constraint  not  to  be  violated,  the  extreme 
signal  the  system  is  to  process  must  be  specified.  The  extreme  command 
input  is  denoted  r (t)  , with  Laplace  transform  R (s)  and  has  the 

S 6 

largest  magnitude  and/or  fastest  dynamics  input  to  be  applied  to  the 
system.  In  general,  a maximum  step  input  will  be  considered. 


Extreme  disturbance  input 


The  desired  disturbance  attenuation  has  already  been  defined.  For 
the  SOAS  and  for  the  other  adaptive  oscillating  systems  developed  in 
later  sections,  it  is  also  necessary  to  define  the  expected  maximum 
equivalent  disturbance  at  the  plant  output.  The  problem  due  to  the 
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V 

* 'l 

i" 


extreme  disturbance  is  more  complex  than  the  one  due  to  the  extreme 

command,  and  in  some  cases  there  is  no  possible  solution  for  certain 

extreme  disturbances  (Section  2.5-3).  The  extreme  disturbance  input  is 

denoted  d (t)  , with  its  Laplace  transform  D (s) 

e e 

2.5-2  The  Constraints 


The  constraints  are  the  quasi-linear  ones  (Section  2.2)  and  another 

rgi 

important  one  first  used  decisively  by  Horowitz  1 J , concerning  the 
maximum  oscillation  component  output,  tolerable  in  any  practical  system. 
The  conditions  developed  in  detail  in  Section  2.2,  Equations  2-10-2-12, 
are  listed  here  for  easier  reference. 


max 

t 


|x.(t)| 


i = r ,d 


co 


bi 


i = r,d 


a = 3 
B = 3 


(2 . 5-5a) 

(2 . 5-5b) 

(2.5-5c) 

(2.5-5c) 


If  these  conditions  are  satisfied,  the  representation  for  the 
non-linearity  is 

Mf  2B 

N = N,  = — (=  — for  ideal  relay)  (2.5-6a) 

r d A A 

M 2M 

N = — = (=  — for  ideal  relay)  (2.5-6b) 

o A A ttA 

B is  the  saturating  level  of  the  ideal  relay  (see  Figure  2.2-2). 


Maximum  plant  output  oscillation 

In  the  SOAS,  with  its  deliberate  oscillation  signal,  there  must 


obviously  be  a maximum  tolerable  oscillation  amplitude  level  at  the 


plant  output,  denoted  by  m 

Ic  (t) | S m 
1 o 1 max 


or  (Figure  2.2-1) 


2.5-3  Synthesis  Procedure  (Smooth  solution) 


The  basic  SOAS  structure  is  drawn  again  in  Figure  2.5-1,  with  N 


an  ideal  relay  with  saturating  level  B 


The  loop  transmissions  are 
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o 4B 
N = — = — 
o A A 


(2 . 5-7b) 


Some  relations  will  be  derived  based  on  the  desired  transmission 
Tr  ( j to)  for  the  command  input  and  its  extreme  value  Re(jw)  . From 
Equation  2.5-1  we  have: 


a Cr 

Tr«">  - g3  i^t(iM) 


Also: 


a C^(3w) 

X (ju)  £ — — — 

r v J NfG2KPh(  joj) 


(2.5-8) 


(2.5-9) 


Combining  Equations  2.5-6,  8,  9,  we  get: 


Xr  RTr(jto) 

T<!ia,)  = MfG2KPh(ju) 


(2.5-10) 


Equation  2.5-10  is  in  the  frequency  domain.  But,  the  quasi-linear 

i i A 

constraint  defined  in  Equation  2.5-5,  max  |x  (t) | 5—  , gives  the  time 

t r a 

domain  constraint: 


xr(t) 


,1=1 
a 3 


There  is  no  "direct"  correspondence  between  constraints  on 


(2.5-11) 

|xr(t) 


and  on 


Xr  ( jui) 


But  some  good  engineering  results  can  be  obtained. 


In  Appendix  II  some  approximate  properties  have  been  developed,  which 

*r 

give  conditions  for  — (jw)  to  satisfy  Equation  2.5-11. 

X 

We  assume  in  this  section  — (iw)  is  "fairly  smooth".  Definition: 

A 

A transfer  function  is  "fairly  smooth"  if  all  its  poles  and  zeros  are 
negative  real.  Actually,  we  would  like  to  consider  any  transfer  function 
F (s)  "fairly  smooth"  if  | F ( jto) | can  be  well-approximated  by  another 


If  (ju>)  | such  that  F (s)  has  only  negative  real  poles  and  zeros. 
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In  Appendix  II  it  is  shown  that  for  a "smooth"  — (ju>)  > if 

Equation  2.5-11  is  satisfied,  then 

lx 


a w 


V <D 


(2.5-13) 


It  is  most  difficult  to  satisfy  Equations  11  and  13,  when  the  oscillation 
amplitude  A has  its  minimum  value  A^  , and  the  extreme  command  input 
rg  is  applied.  From  Equation  2.3-7,  A * A^  , when  the  plant  has  its 
minimum  gain  . Substituting  these  extreme  values  in  Equations  5-10 

and  5-13 


X 

A1 

X 

A£(j“) 

A1 

£ - 

0 

R T (jo)) 
e r 


0) 

v 


0) 

0)  <-2. 
br  6 


V w 


(2.5-14) 


(2.5-15) 


so 


ReTr(ja>) 

^ | 

1 

aa) 

u> 


or 


K 

Mf  G2K2Ph  ( I 5 — • a • to  • | R^T^.  ( jto) 


e r 


V id  <? 


' 8 


(2.5-16) 


Inequality  (16)  assures  non-violation  of  quasi-linearity  constraints, 
due  to  command  input  signals  for  K { K2 

An  expression  similar  to  Equation  16  may  be  obtained  for  the 
disturbance  inputs.  In  Figure  2.1-1,  due  to  disturbance  inputs: 

L (jw)  w 

Vu,*- 


Zd(ju>) 


-D(ju)  • 


1+Lf ( jw) 


(2.5-17) 
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xd(t) 

„ 1 

x j 

A 

5 — , gives 

X(3a.) 

au 


V (o 


(2.5-20) 


Substituting  in  Equations  19  and  20  the  most  difficult  combination  of 


U) 

* V (2.5-21) 

aw  o 


parameters  to  satisty,  ( K = K. 

1 

L ' 

> 

II 

> 

!-■ 

|Xd 

Dg( jw) 

# 

Lf  (jw) 

MfG2KlPh(jW) 

1+Lf ( jw) 

or : 


1 1 K2 

De(jw)  Lf(jto) 

|MfG2K2Ph(jw)|  > jq--a  -«o  • 

1+Lf(jw) 

(2.5-22) 

p 


A difference  may  be  noticed  between  Equations  16  and  22.  In 
Equation  16  the  inequality  is  completely  defined  apart  from  "far-off" 
poles  and  zeros  that  may  perhaps  be  added  to  the  nominal  T^(jio)  . In 
Equation  22  Lf(jw)  is  unknown,  being  the  most  important  result  of  the 
system  synthesis.  In  the  meantime,  let  us  derive  additional  relations 
involving  the  disturbance  response.  The  specified  maximum  disturbance 
response  provides  bounds  on  L^(jw)  for  each  frequency.  Based  on 


Equation  2.5-2,  and  on  the  defined  specifications: 


C0 

|Td(jw)|  £ 

- 

frd(j“)|  * 

T . . ( joo  ) 1 

d min  | 

1+Lf  ( jo) ) 


0) 

V“‘f 


(2.5-23) 


A very  useful  technique  for  design  ^ is  to  use  the  Nichol’s  chart 
(Appendix  I ) and  to  plot  on  it  the  boundaries  for  a discrete  number  of 
frequencies  of  , such  that  Equation  23  is  satisfied.  In  the  most 

general  case  there  are  bounds  on  due  to  plant  parameters  uncertainty 

and  to  disturbance  attenuation,  but  here  the  limit  cycle  and  quasi- 
linearity wipe  out  the  parameter  uncertainty  so  the  bounds  are  due 
entirely  to  Equation  23.  It  has  been  shown  ^ that  the  optimum  loop 
transmission  in  the  sense  defined  in  [5] , is  the  one  that  "stays"  on  its 
respective  bounds,  for  each  frequency. 

For  plotting  the  bounds  on  L^(ju)  at  several  frequencies  in  the 

Nichol's  chart,  once  | T.  . ( jco ) 1 has  been  defined,  the  usual  technique 

| d min  J | 

is  to  make  a change  of  variables  Jl-(jco)  = — i — T » and  use  an 

f J Lf(:io) 

"inverted"  Nichol's  chart,  so  the  closed  loop  contours  with  constant 
gains  can  be  used  for  plotting  the  bounds  due  to  the  disturbance 
attenuation  (Appendix  I ),  i.e.  let 


1 

if (ju) 

|Td(ja,)|  " 1+Lf  (joo) 

1 + *f(jw) 
r 

l+£f ( jw) 

* |Td  min  ^ | 

U 

v»‘f 


(2.5-24) 


A typical  specification  for  Td  • in  graphical  form  is  shown 

in  Figure  2.5-3. 


Figure  2.5-3  - Typical  sketch  for  disturbance  attenuation  specifications 


In  Figure  2.5-3  u is  the  highest  frequency  for  which  it  is 

ri 

important  to  define  the  disturbance  attenuation.  As  per  Equation  2.4 
we  can  plot  the  bounds  for  £f(j(o)  for  several  frequencies  in 
0 < u>  s u)„  . For  example,  for  01  = 10,  , the  bounds  on  the  Nichol's  chart 
If ( jw) 

for  - — - — /■.  --  = -12  dBs  is  plotted  in  Figure  2.5-4  with  the  arrows 

l + <lf(]0)) 

indicating  the  allowed  regions  in  the  "inverted"  Nichol's  chart,  where 

the  minimum  disturbance  attenuation  is  satisfied.  The  same  technique 

is  used  for  o>  = <o  , <0  = o>,  , <0  £ u>  s uic  , <0,  S 10  s 01  , and  plotted  in 

Z 5 4 bo  H 

Figure  2.5-4.  Obviously,  Lf  is  easily  obtainable  from  and 


vice  versa. 
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By  plotting  the  bounds  of  Z^(ju)  symmetrically  to  the  origin  of 
the  Nychol's  chart  ( Log  |a|  =0  , /cp=-180°  , in  the  open  loop 
coordinates),  the  bounds  of  L^(ju)  are  obtained.  This  is  shown  in 
Figure  2.5-5,  where  the  bounds  on  the  several  frequencies  for  L^(ju) 
are  plotted.  If  Lf(ju)  satisfies  these  bounds,  then  the  disturbance 
attenuation  specified  in  Figure  2.5-3  is  achieved. 


The  second  constraint  to  be  satisfied  is  the  maximum  allowable 
plant  output  oscillation  amplitude,  |CQ(t)  | Sm  . From  Figure  2.5-1 

|Co(3V|  = Mo|G2KPh(J"o)|  (2-5_25) 

The  largest  amplitude  of  oscillation  occurs  at  K=K^  , since 
Mo  ' jG2^<0o^|  ' jPh^Wo^  | are  constants-  Hence 

|co(j“o’|™*  * hWh(’v|  ■ |2B£G2Vh«“o>|  * m (2-5-261 

We  can  now  obtain  a relation  expressing  a price  that  must  be  paid  by 
using  SOAS.  Recalling  Equations  2.5-14  and  2.5-21,  and  replacing  w = o>o 


X 

PS 

1 

xd 

r(i“o) 

PS 

1 

R T (jw  ) 
e r J o 


*VWh(>o> 


De(j"o’ 


Lf(j»0) 


1+Lf (j“o) 


(2.5-27) 


(2.5-28) 


Substituting  Equation  2.5-26  in  both  27  and  28,  and  since  L^(jo)Q)  = 


and 


Lf(ju,o) 


1+Lf (jwQ) 


=1  , we  get  respectively 


& 


Vf(j“o) 


MfG2K2Ph(j“o) 


5 2 


K2  Ve^o* 


m 


xh 

A^o1 


D (ju)  ) 
e o 


MfG2K2Ph(^o> 


K.  D ( ju)  ) 


(2.5-29) 


(2.5-30) 


■H 
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In  Equations  27  to  30,  instead  of  the  sign  = (equal) , there  appears  the 

sign  (similar) . The  reason  is  that  the  region  of  validity  for 

u 

Equations  2.5-14  and  2.5-21  is  — , and  we  are  extrapolating  these 

p 

equations  for  w = u)q  • This  extrapolation  is  valid  as  per  Appendix  in. 
Equations  29  and  30  provide  non-linear  functions,  which  define  the 

xr  xd 

minimum  value  of  the  oscillation  frequency  when  a " smooth"  — and  - — 

A1  A1 

are  being  considered. 

K 

It  is  very  interesting  to  see  how  the  factor  - — — , which  was  the 

min 

main  reason  for  using  an  adaptive  solution  (in  place  of  a linear  time 

invariant  design) , reappears  in  the  definition  of  wo  , which  defines 

the  bandwidth  of  the  loop  transmission,  since  | L^.  ( jwo)  j = . The 

reason  for  this  "reappearance"  is  due  to  the  fact  that  the  quasi-linear 

constraint  is  the  most  difficult  to  satisfy  when  the  plant  has  its 

minimum  gain  K = K . = K.  . On  the  other  hand,  the  maximum  oscillation 

mm  1 

level  at  the  plant  output  occurs  when  the  plant  has  its  maximum  gain, 

K = K = K,,  . These  two  constraints  are  such  that  the  SOAS  is  over- 
max 2 

designed  in  these  aspects  at  the  non-extreme  plant  gains  and  command 
and  disturbance  inputs. 

The  second  NB  factor  is  the  ratio  of  a large  signal 

( Z . = TR  or  D ) in  the  system  to  a small  signal  level  - m . The 
ei  r e e 

product  of  these  two  factors  gives,  so  to  speak,  the  order  of  magnitude 

of  the  SOAS  design  problem.  Also, 

Zo  (ju  ) = R (j<*>  ) .T  (ju  ) (2.5-31) 

er  o e o r o 

Where  Zgr  represents  the  extreme  plant  output  due  to  the  extreme 
command  input  rg(t) 


Nr-  ' 
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Equation  29  becomes 


2K2 

* K,  m 

1 

Similarly 

Z , ( jco  ) = D ( j<o  ) x 
ed  J o e J o 


Lf{jV 


= D (ju>  ) 
e o 


(2-5 


(2.5 


1+Lf (j<0Q) 

Where  Z , represents  the  extreme  plant  output  due  to  the  extreme 
ed 


disturbance  d (t) 
e 


Equation  30  becomes 


x j 

r‘(j“o) 

i 


> 2^2  JZed(jV 
' K1 


m 


(2.5 


Let  Z . represent  the  extreme  forced  plant  output  components  ( i 
ei 

Equations  31  and  33  may  be  combined  into 

lx 


7^(jw  ) 
o 


K Z (ju>  ) 

>2  - -■-= — 

K m 


(2.5 


Where  X is  the  extreme  forced  signal  at  the  non-linearity  input, 
e 

X 

G 

Still  under  the  assumption  of  a smooth  — (jw)  , we  can  write 


So 


X 

r(j“} 

Ai 


au>  3u> 


k2  Z (j«o) 


2 rr- 


m 


X 

i 


3w 


(2.5 


(2.5 


A typical  sketch  for  u)Q  determination  is  plotted  in  Figure  2 


-32) 

-33) 

-34) 

r,d  ). 

-35) 

-36) 

,-37) 

i.5-6. 


o' min 


Figure  2.5-6  - Typical  sketch  of  inequation  2.5-37 
for  a smooth  solution 


In  some  problems  there  is  no  possible  solution,  for  some  disturbance 


signals,  without  violating  quasi-linear  constraints.  This  can  be  seen 


From  Equation  2.5-21  at  w = u 


XH 

i 

i 

£ 

aa) 

= 

3u) 

i 

o 

o 

-38- 


K_  , 

_ , . . i m 1 1 

D (ju  ) £ — • — * — 

e o 6 K_  w 
1 2 o 


(2.5-40) 


Equation  40  gives  the  condition  Dg(juo)  has  to  satisfy,  to  enable  a 
solution  for  the  SOAS,  without  violating  quasi-linear  constraints. 

K2 

A typical  situation  is  shown  in  Figure  7,  ( —=100  , m=0.6  ) . 

K1 

Because  u>o  is  unknown,  this  equation  is  a function  of  w , i.e.  as  if 
u>  is  being  tested  if  it  can  be  used  for  u 


Suppose  dgl  (t)  = (1  - e fc).v(t)  so  Del<s)  = 


Del(jui)  =• 


. This  is  plotted  in  Figure  7,  giving  a minimum 


u)/l+u) 


oscillation  frequency  value  * u)q1 5 1000  rd/sec.  If  the  extreme 

disturbance  is  a step  function  d^  (t)  = . l.y(t)  , so  |De2^u^|  = » 

mKj^ 

D _(joj)  is  parallel  to  - — — , and  there  is  no  u which  satisfies 

e2  J c 6o)  K_  o 

o 2 

Equation  40.  So,  if  there  is  a possibility  of  having  a step  disturbance 
at  the  plant  output  (this  is  usually  not  the  case) , the  SOAS  cannot  be 
a recommended  system,  since  the  maximum  allowed  step  would  be 

- — - 1 
estep  ~ 6K^  1000  ' 


in  our  case,  a ridiculously  small  value,  compared  to  the  oscillation 
component  m , at  the  plant  output.  This  problem  does  not  exist  for 
command  inputs,  since  by  rewriting  Equations  29  and  21  at  u)  = o>o  , 
we  get: 


T R ( jo>  ) 
r e J o 


(2.5-41) 


t:  * 


Even  if  rg(t)  is  a very  large  step  function,  Tr(jw)  will  allow 
a solution  for  Equation  41.  An  example  shown  in  Figure  7,  with 


r (t)  = 100. u(t) 
e 


Re(jlD)  = 


Tr  <s>  = 7 =*  7 

(1+s)  l + 3s  + 3s  +s 


As  0)  -*■  00 


Tr  ( ju>)  - — J 
0) 


So  R T ( jco) 
e r 


, which  results  in  w > u -*33  rd/sec 

O o2 


Supposing  there  is  a solution  for  the  oscillation  frequency  (the 

limitation  exists  if  the  disturbance  input  is  a step),  then  Lf(jw)  is 

the  one  which  lies  on  the  bounds  due  to  the  desired  disturbance 

attenuation,  and  has  L.(jo)  ) = - . w is  the  minimum  frequency 

r o 2 o 

allowing  the  above  conditions,  and  satisfying  u>  ^ Bu>  , u>  * Bu>u 

o r On 

u>_  is  the  transmission  bandwidth  and  w is  the  highest  frequency  of 
r H 

importance  in  the  defined  (Figure  3) . 


2 . 6 Numerical  Example  - Smooth  Solution 


In  this  example  the  procedure  defined  in  2.5  is  used. 


2.6-1  Specifications  and  Constraints 


Desired  transmission  for  the  command  input: 


Tr(s)  = — — r 

1 + — — 'S  + -2-T 
-1  . 1^ 


* 4>(s) 


(2.6-1) 


Where  <5>(s)  includes  far-off  poles  and  zeros  that  can  be  added  to  the 


"nominal"  desired  transmission  without  significantly  affecting  the  time 


response. 
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1 


1 


I 

i 


Note  there  are  no  specifications  for  oo  > w =1.2  rd/sec. 

H 

- The  Plant  - 

P(s)  =|  1 = Kj  J K f Kj  = 100 

- Extreme  Command  Input  - 

r (t)  = 75  u (t)  — R (s)  = ^ 

6 6 O 

- Extreme  Disturbance  Input  - 

The  extreme  disturbance  is  a step  input.  In  this  problem  the 
maximum  possible  amplitude  of  this  step  should  be  found  (without 
violating  non-linearity  constraints) . Also  - the  disturbance 
response  is  supposed  to  be  non-oscillatory . 


- Quasi-linear  Constraints  - 

The  used  parameters  are  u=B  = 3 . So: 


I 

! 


4 


Xr(t)  1 

xd(t)  1 

0) 

0) 

A * 3 ' 

A ' 3 ' 

Wbr  * 3 ' 

“bd  * 3 

- Plant  Output  Oscillation  - 

|c  (t)  I $ m = 1 

I o | 

- Excess  of  Poles  over  Zeros  - 

L^(joi)  is  to  have  an  excess  of  5 poles  over  zeros. 

2.6-2  The  solution 

In  the  Nichol's  chart  of  Figure  2.6-2,  the  effective  bounds  on 
( jui ) are  plotted.  They  are  due  to: 

1.  Disturbance  attenuation  specifications  as  per  Figure  2.6-1,  and 

valid  for  io£u)  =1.2  rd/sec  . The  technique  used  is  the  one 
H 

explained  in  Section  2.5-3,  Figures  2.5-3,  2.5-4,  2.5-5. 


I 


3- 


r 


The  bounds  on  L^(ju)  are  plotted  for  a certain  number  of 


frequencies. 


to  (rd/sec) 


T.  (jco)  . (dB) 
d J min 


2.  L^(juo)  is  the  point  represented  by  A in  the  Nichol's  chart,  with 


amplitude  -6  dB  and  phase  -180°. 


3.  A non-oscillatory  step  disturbance  response  can  be  handled  by  a 
very  good  and  practical  engineering  solution  by  not  allowing 


Lf  ( jco) 


to  be  larger  than  1.3,  or 


Lf  (jco) 


1+Lf ( jco) 

which  is  also  plotted  in  the  same  Nichol's  chart. 


1+Lf  ( jco) 


$2.3  dB 


Next,  using  Equation  2.5-16 

K„ 


M G K P . (jco)  > — • a • co  • R T (jco) 
f22h  er 


Substituting  the  functions  as  defined  in  this  numerical  example 

x 4>  (jco) 


|MfG2K2Ph(ja,)|  5 

22.500 

1 + 10  s + 100  s 

i s=:<o 


(Note  22,500  *■  87  dB) 


For  a first  reference  assume  $(jco)  =1  , and  this  limit  is  plotted 
in  Figure  2.6-3.  Using  Equations  2.5-22  and  2.5-26  at  co  = co^ 

K_ 


— • a co  D (jco  ) 
K o e J o 


Lf(j<0o) 


1+Lf  (jco  ) 


m 


S «fG2K2Ph(j»0)  * 2 


.08 

-32 

. 

.15 

-26 

.29 

-20 

.57 

-14 

1.15 

- 6 

|H 
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Substituting  the  functions  as  defined  in  this  problem, 


300  u D (ju>  ) S — , since 

o e o 2 


or  u .D  ( j u>  ) $ 

o e J o 600 


l+Lf(]'JjQ) 


Assuming  the  disturbance  a step  function 


D (joo)  = -7— 

e J jcu 


De(j“o)  = ~ 

o 


So  the  maximum  step  disturbance  allowed  is  dg(t)  = 60C)  W ^ or 

D (s)  = —7; — , which  is  extremely  small, 

e 600  s 


Using  Equation  2.5-22  again,  at  w 


, k2  , , V^o5 

MfG2K2Dh(jV  ^ — -a  o»oDe(j%)  1+L~("j"a)' ') 
1 i 1 1 r o 


5 2 


(-6  dB) 


Since  a)  is  still  unknown,  this  limit  is  also  plotted  in  Figure  2. 
o 

as  a function  of  w , but  it  is  valid  only  for  w = 


Using  Equation  2.5-26  at  u) 


VWh(i"o>  S 2 


(-6  dB) 


The  intersection  of  the  limits  due  to  Equations  2.5-16,  2.5-22,  2.5-i 
when  <Mju)  = 1 , defines  u>ol»>20  rd/sec  . If  has  extra 

allowed  poles,  the  oscillation  frequency  can  be  decreased. 


Suppose  j <t>  ( jcu)  J = 


(1+s)  | s= jo)  f 


Lf  (ju>) 
1+L  ( jui) 


Then  in  Equation  2.6-16,  at  u> 


i . .i  22 . 500  „ 1 

M G I P i]U  ) > -T  x j 

I x * * n 1 Jl+j.0  s + 100  s Js=jio jU+s)  |_s=ju) 


I, in 
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-r* 


v 

* '! 

I 


The  intersection  of  all  the  limits  defines  to  =4.6  rd/sec 

o 


If  extra  poles  can  be  allowed  in  the  transmission,  u could 

o 


eventually  be  decreased.  It  is  assumed  that  no  more  delays  can  be  added 


to  the  response  time.  Note  that  w >3  to  =3.6  rd/sec  , where  w is 

O H H 


represented  in  Figure  2.6-1.  This  new  constraint  w =4.6  rd/sec  at 

o 


the  point  A (-6  dB,  -180°)  is  plotted  back  in  Figure  2.6-2. 


L^(jto)  has  to  satisfy  all  the  constraints  in  Figure  2.6-2.  It  has 
[5] 


been  shown  that  the  optimum  loop  transmission  in  a certain  sense, 
is  the  one  that  "stays"  on  its  respective  bounds  for  each  frequency. 
In  this  case 


Lf(s)  = 


2-7  (^7 75  * A) 


2 \2 


Which  "stays"  on  the  bounds  for  the  specific  frequencies  (including 
to  =4.6  rd/sec  at  -6  dB,  180“ ) and  it  has  an  excess  of  5 poles  over 


reros.  Lf(s)  being  defined,  Equation  2.5-16  is  used  again 


MfG2K2Ph{jto) 


R .G 


Lf (jto) 


e 3 1+Lf(jto) 


But 


T_  (S) 

g3(s>  ' TOiT 


(1  + 10  s + 100  S2)  (1  + s)2 


1+Lf (S) 


So 


MfG2K2Ph(jw)  5 


22.500 


(1  + 10  s + 100  s2)  (1  + s)2 


Lf (jto) 


1+Lf (jto) 


This  limit  is  plotted  in  Figure  2.6-4. 
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Next  Equation  2.5-22  is  used 


K2  Lf(jw) 

MfG2K2Ph(ju))  * ~ * “>l 


e l+Lf(j«o) 


, Lf(jw) 

WWW  * T 


Which  is  also  plotted  in  Figure  2.6-4. 


As  mentioned  before  (due  to  Equation  2.5-26) , at  , 


Mf°2K2Ph<i“o>  * 2 


We  choose 


3.200  (1 + 10  s) 


M G K P (s)  = — : 

s (1 + 10  s + 100  s ) (1  + s)' 


Since  Lf(S)  = ~ K2ph<s> 


G L (s) 

— (s)  = 

A2  G2MfK2Ph(S) 


-i(s)  = -i- 
A ' 1200 


(1  + 10  s) 


(1  + 10  s + 100  s2)  (1  + s)2  (l  +j|  +^r-j 


By  defining  B = 16  tt  , A2  - 50  (so  A^^.5  ) 


,,  2 B , 100 

M,  = — — = 32  , and  K_P,  (S)  = 

f tt  2ns 


(l  + 10s)  1+^.+^ 

G (s)  - v - ■ ■ . . . . / 

(1  + 10  s + 100  S2)  (1  + s)2 


(1  + 10  s + 100  S2)  (1  + S)2 


G1(S)  = 24 


(1  + 2s  _sl\ 

V1  15  225; 


2 \2 


G3 (s)  has  already  been  defined: 

G (s)  = 2 2 

(1  + 10  s + 100  s ) (1  + s) 

The  complete  system  is  now  fully  defined. 


2.6-3  Simulations  and  Results 


Computer  simulations  were  performed  on  the  CDC  6400  system  at  the 

University  of  Colorado,  using  "Mimic"  program.  Command  and  disturbance 

step  responses  are  calculated  for  three  different  plant  gains  K = = 1 , 

K = K = 10  , K = K = 100  . The  command  input  has  always  its  extreme 

mea  2. 

value  re(t)  = 75  y (t  - 10)  , being  applied  at  time  t = 10  seconds  , 

after  the  limit  cycle  has  achieved  its  steady  state.  The  experiments 
last  50  seconds  each. 


In  Figures  2.6-5  ( K = 100  ) , 2.6-6  ( K=10  ) and  2.6-7  ( K=1  ) , 
RN  is  the  applied  command  input,  X18  represents  the  input  to  the  non- 
linearity and  XC  the  system  output.  For  a better  presentation  of  the 
results,  there  are  only  two  variables  in  each  figure,  corresponding  to 
the  same  experiment  (subtitled  by  & and  b.  ) . 


The  applied  step  disturbance  at  the  plant  output  is  .012,  which  is 


larger  than  the  expected  (conservative)  calculated  dg(t)  , considering 


is  plotted  in  Figure  2.6-8  and 


xd 

jr(j“) 

has  a smooth  shape. 

Xd 

r(ju.) 

Ai 

1 

*d  % 

is  larger  than 

1 

3u> 

1 

shows  how 


in  a certain  region  of 


i 


is  smaller 


l 


u 

f M 

K At 


1 


« 


2 . 7 Synthesis  Procedure  (Non-smooth  solution) 


limit  was  established 


Suppose  now  a non-smooth 


In  this  case  (Appendix  II) 


max 


max 


which  is  much  larger  than 


, which  is  the  limit  for  the  smooth 


; 


Figure  2.7-1  - Comparison  of  two  signals 
with  the  same  "peak"  time  response 


1 
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It  seems  at  first  that  there  is  no  limit  on  how  small  it  is 

possible  to  define  the  oscillating  frequency.  But,  there  are  two  points 

[9] 

to  be  considered  here.  First  the  complexity  of  the  band  pass  filter 
increases  with  the  decrease  in  the  oscillating  frequency,  since  the 
narrowness  of  the  band  pass  filter  augments.  This  is  represented  in 


Figure  2.7-3,  where  X represents  the  filter  complexity. 


Figure  2.7-3  - Filter  complexity  vs.  oscillating  frequency 

T 


V 

% 

* '! 

I 


Second,  the  narrowness  of  the  band  pass  filter  defines  the  values 
of  the  damping  factors  used  in  the  filter,  and  the  tendency  here  is  to 
decrease  these  values.  If  there  are  filters  (poles  mainly)  with  low 
damping  factors,  and  the  command  or  disturbance  signals  have 
"spectrum"  components  at  wq  , the  "real  information"  existing  in 


: -63- 

xQ(t)  is  masked.  This  masking  is  due  to  the  long  decay  time  and  signal 
amplification  (at  ) of  band  pass  filters  with  low  damping  factors. 

The  techniques  used  in  Chapter  2.5  can  be  applied  to  the  "non- 
X 

smooth"  - — (jco)  , with  minor  modifications. 
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CHAPTER  THI 


THE  EXTERNALLY  EXCITED  OSCILLATING  ADAPTIVE  SYSTEM 


(E.E.A.S.) 


3.1  General 


The  limitations  of  the  SOAS  have  been  presented  in  Chapter  2.  To 
avoid  some  of  these  limitations  it  is  possible  to  "inject"  an  external 
periodic  signal  at  some  points  of  the  loop,  and  the  system  will 
still  have  its  zero-sensitivity  properties  to  the  plant  gain  factor. 
This  system  is  called  - The  Externally  Excited  Oscillating  Adaptive 
System  ( EEAS ) - and  its  structure  is  shown  in  Figure  3.1-1. 

"3  (t) 


Figure  3.1-1  - EEAS  Structure 


In  Figure  3.1,  G^s)  , G2(s)  and  G3(s)  are  the  transfer 

functions  of  the  linear,  time  invariant  compensations  available  to  the 

designer  and  as  usual  P(s)  is  the  plant  transfer  function.  N is  the 

introduced  non-linearity.  The  signals  x(t)  = x (t)  + x (t)  + x--(t)  , 

o r d 

where  xQ(t)  = A sin(ioot  + <p)  , y(t)  = yQ(t)  + yr(t)  + yd(t)  , where 


S&KV' 


y (t)  = M sin(o)  t + cp)  . w(t)  , z(t)  , c(t)  , d(t)  and  n(t)  have 
o o o 

the  same  meaning  as  in  Chapter  2 for  the  SOAS.  If (t)  is  the  injected 

periodic  signal  with  period  Tq  , which  has  also  to  be  determined, 

where  <jj  = —■  is  its  first  harmonic  equivalent  frequency,  and  Aq  i 
o 

first  harmonic  component  amplitude. 

The  validity  of  the  describing  function  approximation  depends  on 
"how  sinusoidal"  the  signal  is  at  the  input  to  the  non-linearity. 
G^G-jP (s)  has  tc  have  low-pass  characteristics,  so  the  high  harmonics 
components  of  y(t)  (due  to  non-linear  characteristics  of  N ),  and 
of  »/(t)  (in  case  lf(t)  is  not  a sinusoid) , are  filtered  out. 


3.2  Model  of  the  non-linear  element  and  the  constraints  involved 


As  in  the  SOAS  and  with  the  notation  of  Chapter  2,  the  non-linear 
element  can  be  represented  by  its  Dual  Input  Describing  Function 
(D.I.D.F.) . 


N = — for  the  oscillating  signal 
o A 


N = — for  the  signal  due  to  the  command  signal 

IT 


(3.2-1) 


(3.2-2) 


N , = — for  the  signals  due  to  the  disturbances 
d A 


(3.2-3) 


’ o A 

M = M = = M 

d r 2 f 


N _ £ N , = N 
f “ d r 


M = — for  an  ideal  relay  with  saturation  level = B 


(3.2-4) 


(3.2-5) 
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This  representation  is  valid  if  two  conditions  are  satisfied: 

- Quasi-linearity  constraints 

- Non-existence  of  a self-oscillation 

Quasi-linearity  Constraints 

These  constraints  are  identical  to  those  of  the  SOAS  (Section  2.2). 
Using  the  same  notation  we  recall  that: 


lx  (t)  I 
1 r 

£ 

A 

a 

(3.2-6) 

lxd(t) 

* 

A 

a 

(3.2-7) 

0) 

“b  * 1 

J 

(3.2-8) 

Non-existence  of  Self-oscillation 

Suppose  there  is  no  injected  signal  in  the  structure  represented 

in  Figure  3.1,  i.e.,  jf(t)  =0  . The  system  will  oscillate  at  the 

frequency  w , 3 |l  ( jco  ) =-180°  . When  we  apply  the  injected 
ir  | o n 

signal  at  frequency  uq  , the  self-oscillation  with  frequency  is 

not  desirable,  because  in  general  the  Triple  Input  Describing  Function 
(T.I.D.F.)  and  the  system  properties  are  more  complicated  if  both 
periodic  signals  are  present. 

The  injected  signal  must  "quench"  the  self-oscillation.  It  has 
[12  ] 

been  shown  that  for  saturating  non-linearities  the  quenching 
condition  is 

M 

-Lo(jlV  = |Lf(j“ii,|  < p 


(3.2-9) 


Which  actually  means  a minimum  gain  margin  condition  for  the  loop 
transmission.  For  the  ideal  relay  p=1.49  . Another  condition  for 

self-oscillation  quenching  is  presented  at  the  end  of  Section  3.4-2. 


3.3  Mathematical  Relations 


Some  relations  will  be  developed  in  this  section.  These  relations 
will  be  helpful  in  proofs,  synthesis  procedure  and  results  inter- 
pretations. Suppose  in  Figure  3.1-1,  d(t)  =n(t)  =r(t)  =0  , and 

L(jaj)  =G^G2P(jw)  defined  and  known.  We  assume  the  system  has  achieved 
its  steady  state  response  due  only  to  the  injected  signal  ^f(t)  , with 

first  harmonic  representation  CT(t)  =A  sin(o>  t)  . The  oscillating 

o o 


input  to  the  non-linearity  is  denoted  by  X(ja)Q)  , then: 

-I(jw  ) L( jio  ) 

„ , . . o o 

X(3“o)  1 + N L(ju)  ) 

o o 


(3.3-1) 


Where  the  phasor  I(jw  ) can  represent  A sin  w t , and  at  jw 


X (jw  ) = A e- 
o o 


(3.3-2) 


Then : 


-A  L( jw  ) 
jip  _ o J o 

J ~ 1 + N L(ju)  ) 
o o 


(3.3-3) 


A + A N L(jd)  ) = -A  e L(ju  ) 
0^0  o o 


(3.3-4) 


Replacing  Nq  by  its  value  in  Equation  3.2-1,  we  get: 


A + M h ( ju  ) = “A  e L(j» J 


~ + T^-  L(ju  ) = -L(jw  ) e ^ 
A A O O 

o o 


(3.3-5) 


If  L ( jto  ) =G,G„P(jw  ) , A and  M are  known,  A and  tp  can 

O 12  0 o o 

be  calculated,  which  define  the  oscillating  input  to  the  non-linearity. 
Equation  3.3-5  is  solved  graphically  in  Figure  3.3-1. 


OC  represents  the  complex  number  L(jw  ) , the  linear  part  of  the 

i 

loop  transmission  ■•'hich  is  assumed  known.  In  Equation  3.3-5  a real 

M __ 

positive  number  has  to  be  added  to  L(jw0)  • — ( = OA  in  Figure  3.3.1) 

ife 

■L..  ■ ||H  Ml  — . - _ _ , - - 
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to  get  a complex  number  with  amplitude  L(jo)Q)  j . Graphically  this 
is  done  by  tracing  from  A a parallel  line  to  the  real  axis,  until  it 
intersects  the  circle  with  radius  |L(jWQ)  at  B . OB  is  the 

— -i(0 

complex  number  - L(jco  )e  and  AB  is  the  real  number  — . In 

o 

this  way  tp  and  — are  graphically  calculated.  The  analytical 

o 

solution  is  also  developed.  Let: 


L(  ja)Q)  = L(jwo)  e' 


(3.3-6) 


So,  based  on  Equation  3.3-5 


K M 

r*r 


L(jw  ) cos  0 
o 


)}2  * fc 


L ( jto^)  sin  9j  = L(ju>o) 


(3.3-7) 


(f)2  * 


) cos  0 + 


2 2 

L ( ju)  ) cos  0 + 
o 


/M  \2 

| L(jw0)  1 2sin2  0 = |L(ja)o)|2 


(3.3-8) 


2 ^ 

(r)  + 2(^)*|L(ju)0)|cos  0'(r_) + (”^2  ■ 1)|L(juo)|  =0  (3-3_9) 

N N r\r  * 1 Q A 


From  which: 


r = |(_  fcos  e ± Z1  - {f)  3in/ 

o'  1 o f o' 


(3.3-10) 


— must  be  a real  positive  number, 
o 

We  will  discuss  the  several  possibilities  for  a real  positive 
A 

solution  for  in  Equation  3.3-10. 

o 


Case  1 


cos  050 


-tt/2  ^ 0 5 tt/2 


(3.3-11) 


Subcase  la 


M 

ir<  1 

o 


(3.3-12) 


In  this  subcase  we  can  write 


f\  - (-*)'  sin2  0 > /l  - sin2  6 


cos  9 > - — cos  0 
A 

o 


(3.3-13) 


3o  in  Equation  3.3-10 


/M  \ 

-fe) 


M \2  9 M 

sin  0 > — cos  0 
A / A 

o'  o 


and  there  is  always  a unique  real  positive  solution  for  — 


r-  L(juo) 

o 


( / ' © sin2  9 ' t cos  9) 


(3.3-14) 


which  is  the  solution  of  Equation  3.3-10  for  Subcase  la. 


Subcase  lb 


(3.3-15) 


In  this  subcase  we  can  write: 


/ /M  \2 j f 2 ' Mo 

fl  - (— ) sin  9 = /l  - sin  0 = cos  0 = — cos  9 
\A„/  o 


(3.3-16) 


-=  L(ja)o) 
o 


— ' COS  0 ± T~  cos  0 
\ A , 
o o 


(3.3-17) 


which  provide 


(3.3-18) 


■>7j> 


1 
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— = L(jo)Q)  (-2  cos  0) 
o 


(3.3-19) 


Since  cos  9 > 0 , the  solution  given  in  Equation  3.3-19  is  negative 


and  there  is  no  real  positive  solution  for 


in  Subcase  lb. 


Subcase  lc 


M 

it** 

o 


(3.3-20) 


Mo  /M  \ 2 

In  this  subcase,  if  — sin  9 > 1 , then  1 - sin  < 0 

o ' o ' 

M 

and  there  is  no  real  solution  for  Equation  3.3-10.  If  sin  0 £ 1 then 

o 


M \ 2 


2 / ~ — Mo 

sin  0 < /l  - sin^  0 = cos  0 < — cos  0 


In  Equation  3.3-10 


M 

- — cos  0 ± 
A 

o 


M \ 2 


sin  9 < 0 


(3.3-21) 


and  there  is  no  real  positive  solution  for  — in  Subcase  lc. 

o 

To  summarize  Case  1,  where  - /2  £ 0 s tt/2  , the  only  possible  real 

ft 

positive  solution  for  — in  Equation  3.3-10  exists  when  the  injected 

o 

signal  Aq  has  a higher  amplitude  than  the  non-linearity  output 

M 

oscillation  signal  M , i.e.  A > M or  — < 1 

o o o A 

o 

We  can  see  in  Figure  3.3-2  the  same  results  of  Case  1,  based  on 


the  graphical  solution.  In  this  figure  OB  represents  a typical  case 

i‘*‘y 


when  . Subcases  la,  lb  and  lc  are  represented  respectively 


by  OA  , OB  and  OC  . One  can  see  that  only  for  Subcase  la  there  is 

a possible  positive  solution,  represented  by  AD  , which  is  — 

o 


solution  for 


Imaginary  Axis 
S 


,1 


-73- 


So  Equation  3.3-10 

|m_| 


M 

/ 

in 

2 

(-cos  0)  = 

| O 

A 

o 

| COS  0 | £ 

\ 

/- 

> 

0 0 

2 

sin  0 

(3.3-24) 


and  there  will  always  be  only  one  real  positive  solution  for  — 

A 

o 


Subcase  2b 


M 

o 


— | sin  0 | < 1 
o 


(3.3-25) 


i 

In  this  subcase 

f\ 

M 

o 

2 • 2 
sin  8 

1 

o 

/ 

m 

o 

. 2 

1 - 

A 

°l 

sm  0 

and  there  will  be  two  positive  real  solutions  for  — in  Equation  10. 

o 

They  are  (since  -cos  0 = I cos  0|  ) : 


= L(j 


/ M 

’“°’i  (jr 


1 / M 2 „ > 

- I COS  e|  ± A - Sin2  0 

o / | o * 


(3.2-26) 


Subcase  2c 


M 

s2 


(3.2-27) 

~ |sin  0|  =1 
o 

/(  \m\2  “ 

1 sin  0=0  , and  Equation  10  provides 

o 

two  identical  positive  solutions  for  ~ . m this  subcase: 

o 


llT  = - L(j«oo)  | cos  0 

o o 


(3.2-28) 
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Subcase  2d 


M 

o 

IA 

I o 

> 1 

M ! 
° 1 
A 

|sin  9|  > 

° 

t 

M 

o 

21 

2 

. 2 „ 
sin  0 

o 

1 

(3.3-29) 


there  is  no  positive  real  solution  for 


in  Equation  10. 


In  Figure  3.3-2  the  results  of  Subcases  2a,  2b,  2c,  2d  are 

represented  by  OF  , OG  , OH  and  OL  respectively.  In  this  figure 
M 

the  limits  of  — L(jwo)  , which  provide  a real  positive  solution  for 
o 

ft 

— in  Equation  10,  are  also  plotted  and  defined  by  the  segments 
o 



QI  , IK  , KR  . 

The  phase  angle  "ip"  in  Equation  3.3-5  still  has  to  be  calculated. 
For  reference.  Equation  3.3-5  is  rewritten  here: 

r + Jr  L(j<V  = -L(j(V  ejip 

o o 


A 

M 

o 

A 

o 

+ 

A 

o 

-Mju> 


) e - jr  + /l( jo)  ) e^*ir  + 9-  tp 


Combining  Equations  3.3-30,  3.3-31,  3.3-32,  we  can  write 


tg(e-tp) 


Due  to  Equation  3.3-30,  we  can  also  write 


tg(9-<p) 


sin(9-4>) 


.it::on  3.3-32  gives  the  solution  for  (p  in  Equation  3.3-32.  As 


i,  due  to  the  saturating  non-linearity,  4)  is  not  a function  of 


Jfml” 

Li 

A 

A 

_£J 

+ 

M 

o 

A 

o 

L(jtOQ) 

Re 

A 

M 

o 

L(  juO 

A 

o 

+ 

A 

o 

/ 

A 

M 

o 

\ 

L 

1 A 

o 

+ 



A 

o 

L(jo)  ) 

° > 

r a 

M 

o 

L(jUQ) 

I A + 

A 

L o 

o 

M 

A 

o 

I L(  jco  ) 1 sin  0 
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S3. 4 Mathematical  Relations  Interpretation 


The  relations  obtained  in  Section  3.3  will  be  interpreted  in  this 

i 

section.  Basic  properties  will  be  pointed  out  as  a result  of  these 

delations.  When  necessary,  physical  interpretation  of  the  results  wil_, 

i 

ke  given  to  explain  some  "strange  results"  obtained. 


13.4-1  Zero  Sensitivity  Property  of  the  E.E.A.S.  to  Pure  Gain  Changes 
of  the  Plant  P 


We  recall  Equation  3.3-5: 


M 

A 

A 

+ L(jwo) 

o 

A 

o 

o 

-L(jo>o)  e 


-j<P 


put  L(jo)o)  = , P(»  A KPh(jw)  V o> 

ISO, 

W * |”o|[GlS"hll“o1]  - -[GXG21CPh(J"o>]e'W|' 


|Mo|] 


constant 


(3 .4-1)1 


Equation  3.4-1  shows  that  the  factor  — is  independent  of  K 


itself,  because:  Aq  , the  injected  signal  is  constant,  tp  is  define:] 


jby  ^L(  ja)Q)  , Mq  and  Aq  (Equation  3.3-32).  Mq  is  defined 


completely  by  the  non-linearity  saturating  level  used.  This  means  that 


•when  K = K = K_ 
max  2 


A will  have  its  maximum  value  A = A = A„ 

max  2 


[When  K = K . = K. 

mm  1 


A will  have  its  minimum  value  A = A . = A, 

mm  1 


The  loop  transmission  for  the  forced  and  disturbance  signals  is 
(Figure  3.1  and  Equations  3.2-1  - 3.2-5): 


Vi“>  -T-TWh'W-J 


(3.4-2 


I 

t 


r -77- 1 

L J 


i . e . , the  EEAS  properly 


and  is  consequently  independent  of  K , 
designed  also  has  zero  sensitivity  properties  for  pure  gain  changes  of 
the  plant.  This  property  exists  only  if  the  injected  signal  is  applied 
at  any  point  between  the  non-linearity  output  and  the  plant  input.  To 
prove  so,  we  will  make  another  possible  block  diagram  for  the  EEAS 
structure  in  Figure  3.4-1,  where 

G21(S)  X G22(s)  = G2(S)  (3.4-3) 

The  injected  signal  is  applied  at  the  input  to  G22(s)  ' The 

structure  represented  in  Figure  3.1-1  is  a special  case  of  the  structure 
represented  in  Figure  3.4-1. 

! 


,■ **W  *•* 


Figure  3.4-2  - Equivalent  System  for  the  Generalized  EE AS 


sensitivity  property  to  pure  gain  changes  is  kept 


Let  us  assume  now  that  the  injected  signal  is  applied  at  any  point 


in  the  loop  between  the  system  output  and  the  input  to  the  non-linearity 


The  block-diagram  of  this  configuration  is  represented  in  Figure  3.4-3 


are  linear  time  invariant  compensations,  such  that 


Figure  3.4-3  - Injected  signal  applied  between 
the  system  output  and  the  non-linearity 


oscillation  due  to  the  injected  signal 
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or 


or 


X(j“0)  = I(3Vg12(3V  - y(jU)o)*G2KPhGllG12(3U,o) 


j (w^t  + <P)  j (“0t)  j (w  t+<P) 

|A|e  = I Ao|e  *G12(j“0)  - I M0| e 


(3.4-9) 


I A|eDCP  = \ •G12(j03o)  - |M 


or 


A 

K 

1o|e3'P(Gl02Kph(j"o)) 
G12<J“o>e'j,,>  - |“o|(GlG2KPh(i“o>) 


(3.4-10) 

(3.4-11) 

(3.4-12) 


Equation  3.4-12  shows  that 


is  not  independent  of  K itself.  So 


the  loop  transmissions  for  the  forced  and  disturbances  signals 

Lf(jw)  - 2A  WV**0  = (l)fGlG2Ph(jW) 


(3.4-13) 


is  a function  of  the  gain  of  the  plant,  and  the  zero  sensitivity 
property  is  not  kept  when  the  injected  signal  is  applied  as  in 
Figure  3.4-3. 

A last  remark  - if  the  injected  signal  would  have  its  amplitude 

|a_  ’ 


variable  and  proportional  to  the  gain  of  the  plant  K , then 
in  Equation  3.4-12  would  be  a constant  value  and  zero  sensitivity 
properties  would  be  kept.  For  more  details  see  Section  6. 


o 

K 


3.4-2  Physical  Interpretation  of  the  Mathematical  Relations  Developed 

In  Section  3.3  some  relations  were  developed,  which  result  mainly 
from  the  solution  of  Equation  3.3-10.  The  solution  of  this  equation, 
depending  on  the  several  parameters,  sometimes  gives  "strange  results". 
If  the  signal  - = Aq  sin(u)ot)  is  injected  at  a certain  point  in  the 
loop,  we  "expect"  to  have  a well  defined  periodic  signal  at  any  point 


V 


-81 


'i 


The  discussion  of  Equation  3.3-10,  which  provides  the  solution  for 
the  amplitude  A of  the  oscillating  signal  at  the  non-linearity  input. 


(x  = A sin(u  t+tp))  shows  that  in  some  cases  we  have  one  defined  value 
o o 


for  A (as  "expected") . In  other  cases  there  are  two  possible  values 
for  A , which  mean  "two  oscillations"  at  the  same  frequency.  And 
finally,  in  some  cases  there  is  no  possible  value  for  A , which  means 
there  is  "no  oscillation"  component  at  the  non-linearity  input  even 
when  the  injected  signal  is  applied. 


The  two  last  possibilities  seem  rather  "strange",  and  we  will  give 
a physical  interpretation  for  these  results.  We  first  repeat 
Equation  3.3-10 


A 

A 

= |L(j«o)|  (- 

M 

o 

A 

cos  9 ± A - 

/ 

M 

O 

A 

o 

o 

o 

sin2  9 


(3.4-14) 


and  represent  in  Figure  3.4-4  the  results  of  the  several  possible  cases 
discussed  in  Section  3.3. 


The  independent  variable  in  Figure  3.4-4  is  the  angle  9_  as 
expressed  in  Equation  3.3-6  9 = t lG2P(jt°o)  ■L  (jWQ)  . In 

Figure  3.4-5,  the  several  possibilities  are  presented  as  a table,  with 
all  cases  involved. 


We  plot  the  same  results  , having 

Ini 


as  a dependent 


, in  Figure  3.4-6. 


variable,  instead  of 


~1® 
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Figure  3.4-6  - Graphical  Representation 
Discussion  of  Results  - Equation  3.3-10 

In  Figure  3.4-6  the  encircled  numbers  in  each  "region"  represent 
the  existing  number  of  solutions  in  the  several  cases,  resumed  in  the 
table  of  Figure  3.4-5. 

The  "strange"  results  can  now  be  physically  explained.  First 

A 

notice  in  Figure  3.4-6  that  for  0 = it  and  if  — =0  , there  is  one 

o 

A 1 

solution  for  Equation  3.3-10.  — =0  means  there  is  no  injected 

o 

signal,  and  the  system  self-oscillates  at  the  frequency  such  that 

0 = /l  (jw  ) = it  . In  this  case  we  are  back  to  the  SOAS.  For  any 
/ o ir 

|a  I 

0 * it  , when  I--;  =0  , there  is  no  solution  for  Equation  3.3-10,  which 

I M 

I °l 

defines  the  amplitude  and  phase  of  the  oscillating  component  due  to 
injected  signal.  If  there  is  no  injected  signal  ( Aq = 0 ) , the  system 
self-oscillates  at  w 

IT 

Suppose  now  we  inject  a very  small  amplitude  for  | Aq j . For 
any  0 * n , there  is  still  "no  solution"  for  Equation  3.3-7,  even  when 


(■ 

t ! 

i 


the  injected  signal  is  applied.  The  interpretation  here  is  that  the 
injected  signal  did  not  "quench"  the  self-oscillation,  and  the  D.I.D.F. 
for  the  non-linearity  is  not  valid,  because  there  are  two  signals  with 
different  frequencies  "circulating"  in  the  loop,  one  due  to  the  injected 
signal,  and  the  second  due  to  the  self-oscillating  frequency  w . If 

IT 

the  injected  signal  amplitude  Aq  is  augmented,  we  will  eventually 
"quench"  the  self-oscillation,  and  in  Figure  3.4-6  a region  with  one 
unique  solution  is  achieved,  as  expected.  When  ^-<0<4^-  ' if  the 
injected  signal  amplitude  is  augmented,  there  is  a region  in 
Figure  3.4-6  with  two  different  solutions  for  the  amplitude  signal  with 
the  same  frequency.  The  physical  interpretation  is  that  one  of  them  is 
unstable  and  only  the  signal  with  the  largest  amplitude  exists. 


The  same  results  could  be  obtained  (for  ip  ) by  using  the  Nichol's 
chart  as  in  Figure  3.4-7.  Based  on  the  EEAS  structure  shown  in 
Figure  3-1  we  can  write: 


Y(j“0)  = 

But 

And 

|I(j(0o)  I 


I ( ju)  ) L (j  0)  ) 
J o o J o 

1 + L ( jw  ) 
o o 


(3.4-15) 

(3.4-16) 

(3.4-17) 


Substituting  in  Equation  3.4-15,  we  get: 


M 

o 

Lo(j“o) 

A 

O 

1 + L (ju>  ) 
o o 

(3.4-18) 


In  the  Nichol's  chart,  shown  in  Figure  3.4-7,  the  straight  lines 
represent  L , the  open  loop  characteristics  (horizontal  lines  have 
constant  amplitude,  scaled  in  dB  ; vertical  lines  have  constant  phase 

*3 
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_ J 


angles,  scaled  in  degrees) , and  the  curved  lines  represent  constant 

amplitude  (scaled  in  dB)  and  constant  phase  (scaled  in  degrees)  of  the 

closed  loop  ■ ^ , 

(1  + L) 

| M_ 

We  notice  that  if 


Lo(i“o) 


1 + L ( ju  ) 
o J o 


< 1 , or  based  on  Equation  3.4-18, 


< 0 dBs  , which  are  represented  by  curved  contours,  there 


is  only  one  value  of  L (jw  ) for  each  value  of  0 
1 o J o 

TT  7T 

For  - — $ Q there  is  no  possibility  of  having 


L ( j(i)  ) 
o o 


1 + L ( ju)  ) 
o J o 


> 0 dBs  , which  means  that  for 


M 

o 

A 

o 

L (jw  ) 
o o 


1 + L ( jco  ) 
O J p 


> 1 


there  is  no  solution  for  Equation  3.4-18.  As  an  example,  in  Figure  3.4-7 
the  vertical  line  representing  0 = -60°  is  drawn  more  heavily.  This 


line  intercepts  the  closed  loop 

|m 


(rfc) 


contours  of 


with  values 


< 0 dBs , which  means 


For  j < 0 < -y-  , if 


< 1 


M 


>1  (>  0 dBs)  , for  each  0 there  are 


two  possible  solutions  for  Equation  3.4-18.  As  an  example,  in 

Figure  3.4-7  the  line  representing  0 = -160°  is  also  redrawn  and  if  we 
M 

consider  — = 2 (46  dB) , two  possible  solutions  are  represented  by 

o 


A and  B in  the  same  figure.  When 
only  one  solution,  represented  by  C 


M 


= -4  dB  (<  0 dB)  there  is 


Figure  3.4-8  represents  the  limit  cycle  build-up  of  a typical  EEAS . 
The  example  used  here  is  that  one  in  Section  3.5,  in  which  the 


synthesized  open  loop  transmission  hds  been  decreased  by  6 dB  . 

M 

In  this  case  0 5 -350°  and  ^<=*0.9  and  there  is  a unique  solution. 

o 

Figure  3.4-9  represents  the  limit  cycle  build-up  when  A^  is  half  of 
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Figure  3.4-8  - Limit  Cycle  Build-up  for  a Typical  EEAS 

XG^  is  the  input  to  the  non-linearity 
T is  time  in  seconds 


Figure  3.4-9  - EEAS  - Two  Oscillations  in  the  Loop  ( uq  and  ) 

XG1  is  the  input  to  the  non-linearity 
T is  time  in  seconds 
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o 

the  previous  value.  In  this  case  — «1.8  (> 0 dBs) , and  there  is  no 

o 

solution,  since  9sa-350°  . The  limit  cycle  "build-up"  is  composed  of 

two  different  frequencies,  one  due  to  the  injected  signal  and  the  second 
due  to  the  self-oscillation.  The  DIDF  model  for  the  non-linearity  is 
not  valid.  In  both  figures  the  non-linearity  input  is  represented. 

3.4-3  EEAS  Sensitivity  to  Dynamic  Changes  of  the  Plant 

It  has  been  shown  that  both  the  EEAS  and  the  SOAS  have  the  zero 
sensitivity  property  to  pure  gain  changes  of  the  plant.  It  was  also 
shown  (Section  2.4)  that  for  the  SOAS,  if  the  plant  P(jw)  = KP^(jio) 
changes  its  dynamics  characteristics  ( P^(jw)  changes),  then  the 
oscillating  frequency  changes  in  such  a manner  that 

/GlG2Ph(j“o)  = _180°  (3.4-19) 

is  satisfied. 


* 

I , 

i 


The  situation  is  different  for  the  EEAS,  since  the  oscillating 
frequency  wo  is  that  of  the  injected  signal,  if  "quenching"  is 
achieved.  What  happens  when  the  plant  phase  J P^ ( j(*)Q)  changes? 

The  non-linearity  used  here  is  such  that  its  representation  is  a pure 
gain  with  zero  phase. 


Repeating  Equation  3.4-18 


M 

o 

Lo(j“o) 

A 

= 

1 + L ( 

joj  ) 

O 

o 

o 

(3.4-20) 


and  A^  are  fixed.  Thus  as  the  plant  changes  its  phase  angle 


M 

o 

Lo(^o> 

A 

o 

- |l  + Lo(ja.o) 

= constant 


As  an  example, 


I 


contour  defined  by 


r -90-  r 


if  Lq ( jo)Q)  is  represented  in  Figure  3.4-7  at  D , 

/ Lq ( jtoQ)  \ 

( — — — = -4  dB  , /l  ( jto  ) = -60°  ) and  there  is  a change  in  the 

\ 1 + Lo(]Uq)  / oVJ  o'  ) 

phase  of  the  plant  at  ti>o  of  -100°,  the  new  L^jto^)  is  represented 


at  C 


L (ju)  \ 

- ■ ' - = -4  dB  , /L  (jto  ) = -160°  . 

/ ° °-  / 


Suppose  now  a situation  where  — > 1 , as  represented  by  B 


( j °0  ) \ 

— ; — — = +6  dB  , /l  ( j tu  ) =-160°  ) , and  there  is  a change  in  the 

+ Lo(]Wo)  / o o / 


phase  of  the  plant  at  of  +100°.  At  / Lo(ju)o)  =-160°  +100°  =-60° 


M 

o 

A 

o 

Lo(3aio) 

1 + Lo(^0) 


= +6  dB  , so  the 


system  would  eventually  begin  to  self-oscillate  and  the  DIDF  represent- 
ation for  the  non-linearity  is  not  valid  anymore.  This  is  true  for  a 
static  saturating  non-linearity,  where  the  change  in  phase  of  the  plant 
at  0)^  is  reflected  directly  as  a change  in  0 , the  phase  of  the  loop 

transmission  at  the  frequency  to 


This  is  not  necessarily  true  if  a dynamic  saturating  non-linearity 
would  be  used,  represented  by  its  amplitude  and  phase.  This  could  open 
new  possibilities  for  design  and  "insensitivity"  to  P(jt0Q)  changes 
(and  perhaps  also  to  P(jto)  changes).  No  work  has  been  done  on  this 
subject,  but  a possible  non-linearity  is  a relay  with  hysteresis,  as  in 

*5  , _ in  rtiv.  nrnn  ■ j.  [ 12  ] 


Figure  3.4-10.  The  DIDF  is 

N b r -i  /6  + xf 
Nf  = ^7  [sin  \-t~ 


(3.4-21) 


Probably  the  constraints  would  also  have  to  be  changed,  so  quasi- 
linearity conditions  are  kept.  From  now  on  we  will  consider  only 
changes  in  the  overall  gain  of  the  plant. 


Figure  3.4-10  - Relay  with  hysteresis  representation 
3. 5 Synthesis  Procedure  for  the  EEAS 

In  this  section  a synthesis  procedure  for  the  EEAS  will  be  derived. 

It  is  assumed  uncertaintly  only  in  plant  gain  K in  P(s)  =KP,  (s)  , 

h 

K?K2  . The  non-linearity  used  is  an  ideal  relay.  This  synthesis 
procedure  will  allow  the  designer  to  satisfy  the  specifications  and 
constraints,  and  optimize  the  noise  effect  as  defined  for  the  SOAS  and 
in  Appendix  I . 

3.5-1  Specifications  and  Data 

It  is  assumed  that  the  set  of  specifications  and  data  listed  below 
are  known  or  can  be  estimated. 

- Desired  Command  Input  Response  - Extreme  Command  Input 

- Disturbance  Attenuation  - Extreme  Disturbance  Input 


The  Plant 


F; 


i 


\ 

* 

* '! 

[( 


h- 
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result  is  the  condition: 


Lo(j%)  < p 


Lf(*V  ‘ 2 


(3 . 5-3a) 
(3.5-3b) 


p = 1.49  for  an  ideal  relay  (3.5-3c) 
This  constraint  is  given  clearly  in  terms  of  the  loop  transmission  gain 
margin,  since  is  such  that  ^Lq( jui^)  =^Lf(jaiTr)  = -180°  . The 


"quenching"  condition  automatically  takes  care  that 


Lo(3“o) 


1 + L ( jw  ) 
o J o 


has 


a "possible"  solution  as  mentioned  at  the  end  of  Section  3.4. 


Maximum  Plant  Output  Oscillation 


Finally 


[c  (t)  ] $ m 

o max 


( 3 . 5-4a) 


or  (Figure  3.5-1) 


A G,KP  (jw  ) 
o I 2 h o 


1 + L ( jto  ) 
o o 


S m 


( 3 . 5-4b) 


3.5-3  Synthesis  Procedure  (Smooth  Solution) 

The  basic  EEAS  structure  is  drawn  again  in  Figure  3.5-1,  with  non- 
linear element  and  ideal  relay  whose  saturation  level  is  B . As  in 
the  SOAS,  there  are  two  loop  transmissions.  Lq  , , if 


, x^  are  relatively  slow  aperiodic  signals  due  to  the  command  and 

r oj  u) 

o o 

disturbance  inputs.  Lf  is  valid  for  frequencies  less  than  — » — • 


Then 


Lf(jou)  = NfG1G2KPh(jw) 


2B 


N,  = — (for  the  ideal  relay) 
f itA 


(3.5-5a) 
( 3 . 5-5b) 


Figure  3.5-1  - Basic  EEAS  Structure 
Non-linearity  - Ideal  Relay 


First  some  relations  will  be  derived  based  on  the  desired  trans- 
mission for  the  command  input,  and  on  the  extreme  command  input.  All 
these  relations  are  derived  under  the  assumption  that  a "smooth" 
solution  is  used  as  discussed  in  Section  2.5-3  and  Appendix  II. 


Cr  Lf(ju>) 

Tr(j“)  = T(ja))  = G3(ja,)-1  + Lf(j^ 

Cr (ju) 

Xr(jU>  = NfG2KPh(ju)) 

Mf 

Substituting  C = RTf  , Nf = — in  Equation  3.5-7,  gives 


(3.5-6) 

(3.5-7) 


Xr  RTr(ju>) 

T(jw)  = MfG2KPh(joo) 


(3.5-8) 


► 

Equation  3.5-8  provides  a relation  defined  in  the  — ^'uency  domain. 
The  quasi-linearity  constraint  defined  in  Equation  3.5-la  gives  a 


relation  in  the  time  domain 


(3.5-9) 
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xr(t) 


,1=1 
a 3 


which  has  the  same  variables  in  the  left  side  as  in  Equation  3.5-8. 

Condition  3.5-9  is  the  most  difficult  to  satisfy  when: 

a)  A has  its  minimum  value.  A has  its  minimum  value  when  the  plant 

has  its  minimum  gain  K = K . = Kn  . This  has  already  been  shown  in 

mm  i 

Equation  3.4-1.  In  this  case  A = A . = A, 

mm  1 

b) _  x^(t)  has  its  maximum  value.  This  will  happen  when  the  command 

input  has  its  extreme  value  r (t)  . This  is  the  main  reason  the 

e 

extreme  command  input  has  to  be  defined,  as  mentioned  in  Section  3.5-1. 

Again  as  in  the  SOAS,  we  translate  the  time-domain  condition 
expressed  in  Equation  3.5-9  as  a frequency  domain  condition.  Assuming 
x^(t)  has  a "smooth"  shape,  and  using  the  results  of  Appendix  II,  if 


then 


xr  (t) 

< 1 

A 

a 

Xr  ( ju>)  1 

1 

A 

aco 

(3.5-10a) 


(3.5-10b) 


Returning  to  Equations  3.5-8  and  10b  and  replacing  the  parameters 


giving  the  most  difficult  conditions  ( A = A^  , K = , R(jw)  =Re(jw)  ) 


or 


Xr ( jw) 

ReTr(ju) 

0) 

a) 

o 

A1 

MfG2KlPh(ju) 

< — 

6 

Xr(jw) 

_1_ 

V a) 

A1 

£ 

aw 

xr(jw) 

*eVju)) 

I aa) 

! Ai 

MfG2KlPh(ju) 

(3.5-11) 


(3.5-12) 
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or 


K2  “q 

|MfG2K2Ph(ju)  | 5 — * a • “|ReTr(j“)  | V u<~ 


(3.5-13) 


A similar  expression  to  3.5-13  is  derived  for  the  extreme 
disturbance  at  the  plant  output.  In  Figure  3.5-1,  due  to  the 
disturbance  D , 


Lf (ju) 

Zd(ju»  = -D(ju)  • 1 + Lf  ('ju)" 


0)  £ 


6 


X (J.)  . V*!-  ■ 

dVJ  NfG2KPh(jai)  NfG2KPh(3w)  1 + Lf(ju) 


(3.5-14) 


ui-j-  (3.5-15) 


From  here  on  until  Equation  3.5-19  we  follow  the  same  reasoning 
used  in  the  development  of  Equation  3.5-13  for  the  command  input. 

Mf 

Substituting  N^  = — in  Equation  3.5-15 


-D(ioi) 


Lf (ju) 


d 

A (jw)  " MfG2KPh(jw)  1 + Lf(ju) 


(3.5-16) 


The  constraint  max|x,(t) I S-  (Equation  3.5-la)  has  to  be 
t I d | a 

satisfied,  and  the  most  difficult  conditions  to  satisfy  it  are  at 


A = A . = A,  , and  consequently  K=K  . = K,  , and  when  D(ju)  =D  (ju) 

min  1 min  1 e J 

Using  again  the  results  of  Appendix  II  for  smooth  x.  (t)  , if 


xd(t) 


i.i 


(3.5-17a) 


then 


Xd(juj) 


otu 


( 3 . 5-17b) 


Substituting  in  Equation  3.5-16 


X. 

— (j0)) 
A1 

Dg  ( j(*>) 

Lf  ( ju) 

MfG2K1Ph(jw) 

1 + Lf  (ju) 

au 


V 


(3.5-18) 


or 


I 
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& 


«fG2K2Ph(j») 


Lf ( jw) 


V3"'  ' l + L.Cjl») 


V w $ (3.5-19) 

p 


There  is  a difference  in  the  explicitness  of  Equation  3.5-19  if  compared 
with  3.5-13.  In  3.5-13  the  right  side  of  the  inequality  is  completely 
defined  apart  the  "far-off"  poles  and  zeros  that  can  be  added  to  the 
nominal  T^(jw)  . In  Equation  3.5-19  L^(jw)  is  unknown,  being  the 
most  important  result  of  the  system  synthesis.  Let  us  leave  this 
subject  for  the  time  being  and  let  us  work  out  some  useful  relations 
related  to  the  specified  disturbance  attenuation.  The  specified  minimum 
disturbance  attenuation  provides  bounds  on  L^(jw)  for  each  frequency. 
Based  on  Figure  3.5-1, 


Cd(jW) 

|Td(jw) | = 

D ( jw) 

1 + Lf ( jw) 


£ Tdmin  ( jw ) 


(3.5-20) 


[1] 


The  same  technique  used  for  the  SOAS  (Section  2.5-3)  will  be  used. 
Figures  2. 5- 3, 4, 5 are  valid  also  for  the  EEAS . 

The  second  constraint  that  has  to  be  satisfied  is  the  self- 
oscillation quenching  (Section  3.5-2,  Equation  3.5-3). 

(3.5-21) 


0 1.49 

L.(j w ) < £ = (-2.56  dB) 

f 7T  2 2 


This  is  equivalent  to  a minimum  gain  margin  for  L^(jw)  of  2.56  dB. 

In  general  the  disturbance  attenuation  requirement  is  more  "severe"  than 
the  "quenching"  demand,  since  the  disturbance  attenuation  due  to  the 
"quenching"  demand  is 

I T , ( jw  ) I S 10  dB  (3.5-22) 

which  is  not  very  common,  since  it  means  that  the  effect  of  disturbances 
at  w is  allowed  to  be  about  three  times  the  disturbance  itself  (If 

7T 

io(jw  ) I is  small  enough,  it  might  be  allowed). 

7T  I 


f -QR-  f 


The  third  constraint  to  be  satisfied  is  the  maximum  allowable  plant 
output  oscillation  amplitude  (Figure  3.5-1  and  Equation  3.5-4). 


G2KPh(jaio) 


Co(j“o)  Ao  ‘ 1 + L (ju>  ) 


(3.5-23) 


The  maximum  amplitude  of  oscillation  ocurrs  when  K=K  = K.  , giving 

max  / 


Co  ^o  max  1 + L (jo)  ) 

1 o J o 


(3.5-24) 


m 1 + L ( jio  ) 
o J o 


But  (Equation  3.4-20) 


i 

M 

o 

Lo(j“o) 

! 

A 

o 

1 + Lo(j“o) 

(3.5-25) 


1+Lo(3Wo>  Lo(3“o)  Lo(j°)o) 


•X  M x — 
f M. 


1 + Lo(j“o)  Lf(j“o) 


(3.5-26) 


Substituting  in  Equation  3.5-25 


MfG2K2Ph(j“o)  * m Lf(j“0) 


(3.5-27) 


Now  we  will  extrapolate  Equations  3.5-11  and  18  for  u)  = u>o  . This 
extrapolation  is  valid  as  per  Appendix  1 1 1, and  was  also  used  in  the  SOAS 
synthesis  (Section  2.5,  Equations  2.5-27,30).  The  same  is  valid  for 
Equations  3.5-32,35. 

X R T (jo)  ) 

— Mu  ) m 6 r_  - — (3.5-28) 


X 

ReTr(ju,o) 

VA'h'i-o1 

xd 

D (ju  ) 

A[,j“o) 

C VJ 

MfG2KlV3“o> 

+ Lf(jcoo) 


(3.5-29) 


Substituting  Equation  3.5-27  in  28  and  29 


X K„  T R 

r , . , 2 re,. 

1 i 1 f 


(3.5-30) 


X . K9  d0  <3“  ) 

d . 2_  e o 

A -,a)o)  " K.,  in  1 + L,  ( jcu  ) 

1 1 f J o 


(3.5-31) 


Td(^“o)  ” 1 + L (jw  ) 

f J o 


(3.5-32) 


Xd  K2  | 

^(3“o)  " |DeTd(^o) 


(3.5-33) 


Again  as  in  the  SOAS  , — , the  main  reason  for  using  an  adaptive 

K1 

system  instead  of  a linear  time  invariant  system,  reappears  as  a factor 
defining  the  oscillating  frequency.  The  reason  for  this  "reappearance" 
is  the  same  as  for  the  SOAS. 

Notice  that  Equations  3.5-30  and  33  are  similar  to  their  equivalent 
equations  obtained  for  the  SOAS  (Equations  2.5-29,30)  when  in  this  case 


Lf(j“o)  S0AS  = 2 and  W SOAS  = 1 + L,  ( jw  ) 

f J o 


= 2 . Here  lies  the 


main  advantage  of  the  EEAS  when  compared  to  the  SOAS,  since  the 
connection  between  the  oscillating  frequency  and  the  system  bandwidth 
has  been  broken. 

Equations  3.5-30  and  33  could  be  combined  in  one  equation,  which 
can  provide  a good  graphical  sketch  for  a better  understanding  of  the 
problem.  In  Equation  3.5-30 


Z (jw  ) = T R (jw  ) 
er  J o r e o 


(3.5-34) 


where  Z (jw  ) is  the  extreme  plant  output  due  to  the  extreme  command 
er  J o 

input  R^(jw)  at  the  frequency  w^ 


[-3 
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In  Equation  3.5-33,  and  Figure  3.5-1, 


e d J o 


D 

e 

LffjWo) 

Zed(ja)o) 

1 + Lf  (jwQ) 

Lf(j“0) 

Lf ( j^o’ 

(3.5-35) 


where  zetj(3t,l)0)  i-s  tlle  extreme  plant  output  due  to  the  extreme  disturb- 


ance signal  ( jco)  at  the  frequency  (jjq  . Substituting 


Equations  3.5-34  and  35  in  30  and  33 


Xf 

K2 

Z 

A“(ju,o) 

/.  j 

5 K.m  | 

1 

1 

f 

(3.5-36) 


where  represents  the  forced  signal  at  the  non-linearity  input  and 


Zg  is  the  extreme  plant  output,  due  to  either  the  extreme  command  or 


extreme  disturbance  signals. 


For  a smooth  shape  of  x(t)  , as  per  Equations  3.5-12  and  17, 
X, 


/««■»> 

A1 


aw 


(3.5-37) 


Combining  with  Equation  3.5-36 


i 

> h. 

aw 

5 K.m 

CJ 

1 

o 

Z 

“<jw,  0) 

Lf  o 


or 


2aKn 


Lo(jwo)  =2Lf(jwo)  >~ 


w Z (jw  ) 

O 0 J O 


(3.5-38) 


A typical  sketch  representing  Equation  3.5-38  is  shown  in 


Figure  3.5-2,  where  wq  , the  oscillation  frequency,  is  graphically 


determined,  as  the  intersection  of  curves  Lo(jw)  and 


2a  K„ 


— - — • w Z (jw  ) 
mK,  o e o 


wo  is  the  smallest  frequency  satisfying 


Equation  3.5-38. 


J 


Figure  3.5-2  - Typical  sketch  for  oscillation  frequency  determination 

Xf 

EEAS  - Smooth  case  of  — (ju>) 

A1 

This  graphical  solution  only  gives  an  insight  of  the  problem, 

since  Z (ju>  ) and  L (jw  ) are  still  unknown.  Some  more  relations 
e o o o 

a derived.  Comparing  Equation  3.5-19  extrapolated  for  w = uio  and 
Equation  3.5-27, 


!De(j“0,Lf(j“0)| 


— a ID  ' I'.  , r - : TT- ^ MG  K P (jii)  ) $ m L,  ( ju)  ) 

K1  ° 1 + Lf<3W0)  | f 2 2 h o | | f o 

or  comparing  the  two  extreme  terms  of  inequality  39 


(3.5-39) 


T,(ju>  ) = 

d o 


1 

mK1 

1 + L,  ( joa  ) 

1 r ° 

J. 

aK.o) 

2 o 

De(j“o) 

(3.5-40) 


Equation  3.5-40  gives  a constraint,  valid  only  for  u>  = uq  , and  it  is 


equivalent  to  the  one  given  by  the  specified  disturbance  attenuation. 


which  is  valid  for  frequencies  w < w as  per  Figures  2. 5-3, 4,5.  As  an 

H 

-t  Ki  i 

example  suppose  de(t)  =e  , m=l  , — = — , ot  = 3 

...  1 1 
e(]U)  S ( 1 + S)  jaid  + jaj) 

O — juj 


De(ju)  = 


0)V  1 + 0) 


and  using  Equation  3.5-40 


V^o*  * 


This  expression,  valid  only  for  ^ = u0  i is  plotted  in  Figure  3.5-3. 


lVj“> 


u (rd/sec) 


Figure  3.5-3  - Typical  "Disturbance  Attenuation"  at 


The  bounds  of  Lf(ja)Q)  satisfying  the  "disturbance  attenuation 


requirements  at  are  plotted  in  Figure  3.5-4,  using  the  same 


technique  as  explained  in  Section  2.5-3,  Figures  2. 5-3, 4, 5. 
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Notice  in  Figures  3.5-3  and  4,  bounds  have  been  plotted  for  several 

frequencies.  Only  one  of  these  bounds  has  to  be  satisfied  - the  one 

related  to  u . Now  L.(ju)  can  be  defined.  In  the  same  Nichol's 
o f 

chart  the  bounds  on  L^(ju)  , due  to  the  specified  disturbance  attenuation 
(equivalent  to  Figure  2.3-6),  is  superimposed  to  the  bounds  on 
due  to  Equation  3.5-40  (equivalent  to  Figure  3.5-4).  Lf(ju>)  is  the 
loop  transmission  which  is  on  the  boundary  at  each  w , based  on  the 
disturbance  attenuation  specifications,  and  at  u)o  , stays  on  the 
boundary  defined  by  Equation  3.5-40;  , the  injected  oscillating 

frequency,  is  the  minimum  value  which  allows  L^(jio)  to  satisfy  the 


above  conditions  and  to  satisfy  minimum  gain  margin  requirement  due  to 
self-oscillation  quenching  (Equation  3.2-8) . should  also  satisfy 

the  relations: 

u>0  5 8oit  (3.5-40a) 

r 

id  >.  Bwu  (3.5-40b) 

O H 

Where  ujt  is  the  desired  command  input  response  bandwidth  and  is 

r 

related  to  the  disturbance  attenuation  response  as  per  Figure  2.5-3. 


In  this  way  the  noise  transmission  T^(ju>)  is  minimized  in  the  sense 
defined  in  Appendix  I.  The  desired  excess  of  poles  over  zeros  in 
should  be  known. 

Lf (ju) 

Once  L,  ( ju>)  is  defined,  — — — — is  also  defined  and  the 

f 1 + Lf(3w) 

prefilter  G3(jw)  can  be  found  from 


G3( jw)  = T^_ ( jw)  x 


(1  + Lf (jw)) 
Lf  (ju>) 


(3.5-41) 


I 


Next  MfG2K2Ph ^ jU>  to  '3e  defined  _ satisfies  inequalities 

3.5-13  and  19  repeated  here: 


!MfG2K2Ph(-’u)  I * £-*  <*  ’ “ • | ReTr  ( j^)  I 

K2  Lf(ju) 

iMfG2K2Ph(jw)|  De(ju)  * 1 + LfTj^)' 


(3.5-42) 


(3.5-43) 


And  at  u)q  , based  on  Equation  3.5-27 
|MfG2K2Ph(jWo)|  * mlLf(jV| 


(3.5-44) 


Once  MfG2K2Ph<jw)  and  Lf(jto)  are  known,  it  is  possible  to  determine 


G1(ju) 


Lf  (jw) 

WA1*"1 


(3.5-451 


The  ideal  relay  saturation  level  B = — - — and  A 2 are  chosen 
available  voltages,  and  in  practice  should  be  chosen  as  large  as 

possible.  The  injected  signal  amplitude  tf  =Aosin(wot)  is  given  by 

1 + L ( jw  ) 
o J o 

A = M — - — — — as  per  Equation  3.4-18. 

o o L (jut  ) ^ H 

o J o 

Before  we  proceed  on  the  numerical  example,  Equation  3.5-38  is 
repeated : 


L (jw  ) * — — W Z (ju)  ) 

o J o | mK^  o e o 


(3.5-46) 


As  per  Figure  3.5-2,  since  w > w , we  assume  L (jw  ) < L (jw  ) 

o TT  o o o r 


Suppose  now  Lf(juo)  (and  also  LQ(jwo)  = 2Lf(jwo)  ) have  very 
large  magnitude.  In  this  case  Equation  3.5-46  is  much  easier  to  satisfy, 
reducing  the  possible  oscillation  frequency  wq  . This  is  shown  in 
Figure  3.5-5,  and  also  in  the  numerical  example,  Section  3.6.  In 
Figure  3.5-5  three  possible  solutions  for  Lo(jwQ)  are  plotted  with 
•heir  corresponding  oscillation  frequencies.  L . is  the  case  where 


•5 


I 


f ; 

k 


Where  0(s)  is  the  function  including  all  the  far-off  poles  and  zeros 
that  can  be  added  to  the  "nominal"  desired  transmission  without 
significantly  affecting  the  time  response. 

- Disturbance  Attenuation  - 

The  minimum  disturbance  attenuation  required  is  plotted  in 


- Extreme  Command  Input  - 


1 

i 


I 

l 

i 

I 


: 


r-io9- t 


r (t)  = 75  .u (t) 
e 


R (s)  - ^ 


- Extreme  Disturbance  Input  - 

d (t)  * (1 -e“t/3)y(t)  -*  D (s)  = 

e e s (1  + 3s) 

Find  also  the  maximum  allowable  step  disturbance  (which  is  supposed  not 

to  be  oscillatory) . 

- Quasi-linear  Constraints  - 

The  used  parameters  are  a = 8 = 3 , so 


xr(t) 

< I 

x (t) 
d 

< 1 

0)  to 

„ o . o 

A 

* 3 ' 

A 

* 3 ' 

10  « > <0.  , S ~r~ 

br  3 bd  3 

- Plant  Output  Oscillation 


C (t)  £ m = 1 

o 


- Limit  Cycle  Quenching  - 

In  Equation  3.2-8  we  take  p = l , so 

L (jto  ) I < 1 and  I Lc  ( jw  ) I < q- 

O J IT  f IT  2 


- Excess  of  Poles  over  Zeros 


Lf(jw)  is  to  have  an  excess  of  9 poles  over  zeros. 


The  Solution 

In  the  Nichol's  chart  (Figure  3.6-2)  all  the  existing  bounds  of 
Lf(joj)  are  plotted.  They  are  due  to: 

1.  Disturbance  attenuation  specifications  as  per  Figure  3.6-1,  and 
valid  for  u ,<  .3  rd/sec  . 


wilSSS^ 

!**«  wBSSMns 


m^Smgss^^Mgg, 


MfilESU 


t -in-  r 

i 

2.  The  condition  due  to  Equation  3.5-40  can  be  viewed  as  an  "equivalent 

disturbance  attenuation",  and  is  valid  only  for  the  oscillating 
frequency  . Since  is  unknown,  it  is  necessary  to  plot 

this  "equivalent  disturbance"  bound  for  several  frequencies. 

3.  Lf(jio^)  < j (to  avoid  self-oscillation).  It  is  a minimum  gain 
margin  demand. 


4.  The  step  response,  as  requested,  should  not  be  oscillatory.  A good 

Lf (jw) 

engineering  solution  is  not  to  allow  r — — — — > 2.3  dB 

1 + Lf(;ju>) 

All  these  bounds  are  plotted  in  Figure  3.6-2.  First,  the 
constraints  on  L^(ju)  due  to  the  disturbance  attenuation  specifications 


are  plotted  based  on  the  technique  explained  in  Section  2.5-3, 


i Figures  2. 5-3, 4, 5.  The  bounds  on  L^(jw)  are  plotted  for  a certain 

i 

! number  of  frequencies. 

i 


These  bounds  are  plotted,  and  for  each  bound  its  associated 


frequency  and  desired  disturbance  attenuation  are  pointed  out.  Next 


Equation  3.5-40  provides  a bound  on 


j w (rd/sec) 

| 

Td(jw)  min  (dB) 

.01 

-30 

.02 

-23 

.05 

-15 

.10 

- 8 

.20 

- 1 

.30 

+ 2 
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Since  u)Q  is  unknown,  we  will  have  bounds  for  several  frequencies 


but  need  to  be  satisfied  only  at  u 


Substituting  the  values  by  the 


given  data  ( m = 1 , a = 3 , ^=5  , V ju)  = jw  (1  ^ 3jIo)  ) ' 


1 + 9u 


T,  (ja>  ) . 

d o min 


which  provides  the  following  table: 


o)  (rd/sec) 
o 

Td(jWo)  (dB) 

0 

-23.5 

.37 

-20.0 

.72 

-16.0 

1.2 

-12.0 

1.7 

- 9.0 

2.25 

- 6.0 

3.52 

- 3.0 

These  bounds  are  plotted  with  their  associated  frequencies  and  desired 
disturbance  attenuation  marked.  Note  that  if  the  extreme  disturbance 


was  a unit  step,  d^(t)  =u(t)  , the  limit 


/mK  \ 

T^jw  ) . = -23.5  dB  — -=fH  i 

I d J o |min  \aK2  15/ 


independently  of  the  frequency 


(which  corresponds  in  the  present  case  for  u = 0 ).  Next  j ( j ) | < 

is  also  marked  as  the  point  A , with  6 dB  gain  margin  at  the 

Lf ( jw) 

anqle  -180  . Finally  r—— — ■■■  <2.3  dB  is  also  plotted  in  the 

^ J 1 + Lf(3w) 

Nichols  chart. 

Once  all  the  bounds  have  been  plotted  in  the  Nichols  chart, 

I 

L j.  ( )(>))  is  derived,  based  on  the  technique  developed  by  Horowitz  ^ , 
and  summarized  in  Appendix  I. 
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There  are  still  two  conditions  to  be  satisfied: 

w 5 Bid,  , = 3 x .1  = .3  rd/sec 
o or 

u i Bid,.  = 3 x .3  = .9  rd/sec 

O Du 

The  derived  L^(s)  in  our  case  is: 


which  is  also  plotted  in  Figure  3.6-2,  where  the  oscillation  frequency 
u>o=1.25  rd/sec.  Notice  that  one  of  the  poles  has  already  a very  small 
damping  factor  ( £=  .025  ).  In  case  the  extreme  disturbance  was  a unit 
step  input,  and  the  same  oscillation  frequency  is  kept,  a solution  for 
L^(jw)  is  "available"  by  keeping  the  same  Lf(s)  ; the  only  change  is 
in  substituting  the  already  low  damping  factor  for  the  complex  pole 
( <d  = 1.25  , £=  .025  ) by  an  even  lower  damping  factor  ( w = 1.25  , 

C=  .006  ) to  satisfy  the  demmands  of  Lf(j«o)  as  per  Equation  3.5-40. 
Notice  also  that  if  a smooth  loop  transmission  (relate  to  Figure  3.5-5) 
is  considered,  then  tne  oscillation  frequency  could  not  be  less  than 


MfG2K2Ph(ju>)  5 15 


1 1 + 9w 


-115-  I 
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Lf (jw) 

2 l + Lf(jco) 


Lf(jw)  15  Lf(j(D) 

We  first  draw  - — rr— r in  Figure  3.6-3.  Next,  x ; — r 

1 + L (310)  ’ / 1 + L (300) 

✓l  + 9a 1 

is  plotted.  Then,  the  "nominal"  desired  transmission  function 


T (s)  = — = , is  plotted.  Far-off  poles  and  zeros  can 

,!■.!  s'* 

1+— r-s  + ^2 

be  added  to  this  nominal  .transmission. 

Lf  (s) 

All  the  poles  and  zeros  of  — — are  considered  as  far-off 

1 + Lf (s) 

poles  and  zeros  of  T^Cs)  . This  is  valid  since  for  all 

Lf  (s) 

(1)  < 3(o  =3.75  rd/sec  , ■ . ~ «1  . This  makes: 

o 1 + Lf (s) 

1 Lf(s) 

r ' ' m n=  - 1 +Lf(s) 


/U)  = . 1\ 


G (s)  = s-r- 

3 /a>  = . 1\ 


/( O = . -L\ 

{%,) 


Next  the  real  transmission  T^ts)  is  plotted,  and  so  the  limit  for 
MfG2K2Ph^  aS  per  E<luat;‘-on  3 .5-42. 


|MfG2K2Ph^3w)  | i — • a • (o  • |ReTr(ju» 
Substituting  the  values  of  this  example 


MfG2K2Ph(jw)  5 1.125  Tr(jto) 


which  is  also  plotted  in  Figure  3.6-3. 


j MfG2K2Ph  ( j'**)  j has  to  have  a magnitude  (at  all  frequencies)  which 
satisfies  both  conditions,  due  to  the  extreme  disturbance  and  to 


r. 
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the  extreme  input.  | MfG2K2Ph ^as  to  sat^s^y  also  the  condition 
given  by  Equation  3.5-44,  which  guarantees  a maximum  oscillation  level 
at  the  plant  output  ( s m ) . 


MfG2K2Ph( j“o>  * m Lf(:’wo)  = 1 * 2.8  = 2-8 


This  constraint  is  automatically  satisfied  if 
15 


MfG,K  P ( ju)  ) = 

f 2 2 h o 


K_ 


/l  + 9u 


Lf(jo)o) 


1 + Lf (j“0) 


, since  based  on  Equation  3.5-39 


— • a • a)  • D (ju  ) 
K,  o 


o' | 1 + Lf ( juQ) 


= MfG2K2Ph(jWo)  = m Lf(j«o) 


3nd  1 + Lf  < jw0) 

On  the  other  hand,  the  limit  of  |m,G_K_P,  (jw  ) I due  to  the 

I i 2.  c.  fl  O | 

command  input,  is  lower  than  the  limit  of  |m,G„K„P,  (ju  ) I due  to  the 

| t 2 2 h O I 

disturbance  input,  so  all  the  conditions  are  satisfied.  In  the  opposite 
case  there  are  two  possibilities;  either  add  some  more  far-off  poles  to 
G3  (this  means  some  extra  far-off  poles  to  Tr  ) » until  | MfG2K2Ph ^uo^  | 
due  to  the  command  signal  equals  or  is  lower  than  the  limit  due  to  the 
disturbance,  or  in  case  no  more  extra  far-off  poles  are  allowed,  the 
oscillation  frequency  should  be  increased. 


aK_u  D (ju  ) 
2 o e o 


has  been  satisfied. 


To  avoid  the  complexity  of  the  functions,  we  try  to  use  as  many  as 
possible  zeros  and  poles  of  L,  . We  choose 


MfG2K2P(s)  = 


/u  = 1.5\ 

112. u = -5  ; 

■(-ns)  ill'? V*) 


so 


MfG2K2p(ju  ) = 2.6 


Now  we  can  define 


MfG2K2P(3) 


560 


Ao  = 1.25\  /to  = 1.5\; 
Vc  = -025 ; U = -3  ) 


2B 


If  we  take  B = 20  it  , then  M^  = — =40  , and  considering  the  extreme 


plant  K2P(s)  =- 


„ „ „ 200 

MfK2P  = “ ' We  get 


G2(s)  = 


/to  = 1.5' 
.56  U = -5  , 

) 

i (“:;•?( 

Choosing  A2 = 50  (so  A = 10  ) 


G,  = 


.091 


ki£L 


1 /to  = 1.25\  /to  = 1.5V 
U = -025 ) U=  -3  J 


(S)  has  already  been  determined 
1 


G3(s) 


, ,2*.5  s* 


The  injected  signal  3=  A sin  (to  t)  has  a frequency  to  =1.25rd/sec 


o o 

and  its  amplitude  is  available  from 


M 

o 

I 

o 

Lo(jWo) 


1 + L ( jto  ) 
o J o 


= (-1.4  dB)  = (.85) 


Observation  - 


Lo(j“o) 


1 + L (jto  ) 
o J o 


is  measured  in  the  Nichols  chart  when 


L ( jto  ) 
. , . . o J o 

Lf(jt0o)  = ^ 


SO 


4B  4 x 20  7T 


To  .85  ” .85 


80 


.85 


.85 


= 94.116 


Hence  the  injected  signal  is 


(t)  = 94  sin(l. 25  t) 
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Simulation  Results 

The  numerical  example  has  been  simulated  on  a CDC  computer. 

Figures  3.6-4  - 3.6-6  represent  the  signals  for  several  gains  of 
the  plant  ( K=l,  10,  100  respectively)  when  an  extreme  command  input 
is  applied.  Subscript  a represents  the  plant  output,  and  subscript 
b represents  the  non-linearity  input. 

Figures  3.6-7  - 3.6-9  represent  signals  for  several  gains  of  the 
plant  ( K = 1 , 10,  100  respectively)  when  the  extreme  disturbance  input 
is  applied.  Subscripts  a and  b represent  plant  output  and  non- 
linearity input  respectively. 

Finally  simulations  were  also  performed  when  a step  disturbance  is 
applied  (conditions  are  more  severe  than  prescribed  by  specifications) . 

Figures  3.6-10-3.6-12  represent  signals  for  several  gains  of  the 
plant  ( K = 1 , 10,  100  respectively)  when  this  step  disturbance  is 
applied.  Subscripts  a and  b represent  plant  output  and  non-linear- 


ity input  respectively. 


Extreme  Command  Input  Extreme  Command  Input 

K * 100  K = 100 

FX  - Non-linearity  Input  — X FZ  - Plant  Output  Z Z 

Y50  - Command  Input  — fg  Y50  - Command  Input  Z.  I 


CHAPTER  FOUR 


S • 0. A. L. 


[25] 


(SOAS  with  secondary  loop 

controlling  the  adaptive  non-linear  saturating  level) 

4.1  General 

Chapters  2 and  3 have  presented  the  SOAS  and  EEAS,  which  are  self- 
adaptive  systems,  but  have  constraints  which  tend  to  increase  the 
oscillating  frequency  and/or  the  open  loop  bandwidth.  One  constraint 
defined  by  the  maximum  oscillation  amplitude  at  the  plant  output 


MfG2KPh<ju,o)  * m Lf(j“o) 


(4.1-1) 


is  the  most  difficult  to  satisfy  when  K = K =K  for  both  the  EEAS 

max  & 

and  SOAS  ( L ( jw  ) A -i.  for  the  SOAS). 

t o z 

The  other  constraint,  due  to  the  quasi-linear  constraints 


f * k 


(4 . l-2a) 


d 1 

— (t  S - 
A a 


(4 . l-2b) 


will  be  the  most  difficult  to  satisfy,  when  K=K  . = K.  and  the 

mm  1 

oscillating  component  xq = A sin(wot)  at  x will  have  its  minimum 


amplitude  A = A , = A. 

r mm  1 


This  causes  the  oscillation  frequency  w 


and  consequently  the  open  loop  bandwidth  to  be  a non-linear  function  of 
K2 

— , by  means  of  relations  2.5-37  and  3.5-37,  repeated  here. 

K1 


Kjm 


(4.1-3) 


A_  I . 


valid  for  both  the  EEAS  and  SOAS  (for  the  SOAS  Lf(jwQ)  = 2)-  Thus 
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the  main  factor 


we  tried  to  eliminate,  reappears  in  the  adaptive 


systems  in  the  expression  defining  the  loop  transmission  bandwidth. 

This  reappearance  is  inherent  in  the  SOAS  and  EEAS,  because  the  adaptive 
property  exists  due  to  K being  tracked  by  the  amplitude  A , which 
in  turn  causes  the  non-linear  element  to  vary  its  gain  Nf  inversely 
with  A 


Note  that  constraints  4. 1-1, 2 are  such  that  the  system  (EEAS  or 
SOAS)  is  overdesigned  at  any  non-extreme  value  of  the  plant  gain. 

K2 

The  idea  of  the  SOAS  is  to  eliminate  the  factor  — in 

K1 

Equation  4.1-3.  This  factor  appeared  because  A in  xq  = A sin(o)ot) 

2b 

was  used  to  make  N£  = — inversely  proportional  to  K In  the  SOAL, 

A is  kept  (or  at  least  tried  to  be  kept)  constant,  and  instead,  the 

g 

saturating  level  B in  Nf = — tracks  K , i.e.,  B is  made  inversely 
proportional  to  the  plant  gain  K 

Since  the  value  of  K itself  is  unknown,  K cannot  be  directly 
tracked.  The  technique  used  in  the  SOAL  is  to  adjust  B so  as  to  force 
A to  be  constant,  as  shown  in  Figu.-e  4.1-1.  A is  measured  and  compared 
to  a reference  A^  . The  difference  (error)  passes  through  a com- 
pensation network  tjj(s)  , whose  output  controls  B 

The  point  chosen  for  measuring  the  oscillation  level,  which  in- 
directly gives  the  plant  gain,  is  x , the  input  to  N for  two  main 

reasons.  First,  the  system  is  designed  to  meet  quasi-linearity 

v. 

constraints  ut  x , so 

|xf (t) I < j and  (4.1-4) 

maxjx(t)|  = A + lx  (t)l  = (4.1-5) 

t I f I -s 


Figure  4.1-1  - SOAL  Structure'  - Non-linearity  - Ideal  Relay 


This  means  that  if  A is  measured  by  the  maximum  (peak)  value  of  x(t)  , 

the  maximum  possible  relative  error  is  j . In  reality  the  error  in 

measuring  A is  much  smaller  in  a peak  to  peak  detector  measurement, 

because  x^  is  slow  relative  to  x , so  the  x,  component  practically 
f or 

disappears,  except  for  the  change  in  x^  in  half  an  oscillation  period. 

A priori.  Equation  2.5-5  requires  this  change  to  be  small.  Also,  x(t) 
crosses  the  zero  value  twice  per  cycle,  which  is  useful  for  measurement 
purposes . 

Next,  due  to  the  low  pass  band  characteristics  of  Gj^KP  (joi)  , 
xQ(t)  is  supposed  to  have  a "very"  sinusoidal  shape,  which  also  helps 
in  the  amplitude  A measurement  and  consequently  in  the  plant  gain 
identification  accuracy. 
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i 


In  Figure  4.1-1  there  are  "two  loops".  The  main  loop  is  defined  by 
Lf(jw)  = G^G^N^KP^ ( jw)  , which  is  the  same  as  for  the  SOAS  and  EEAS, 
controlling  the  forced  signals.  The  secondary  loop,  which  adjusts  the 
relay  saturation  level,  is  not  "an  instantaneous"  loop,  and  its  dynamics 
and  stability  problems  must  be  analyzed.  This  is  done  in  Section  4.2. 


The  Stability  and  Dynamics  of  the  Secondary  loop  - SOAL 


The  dynamics  and  stability  of  the  secondary  loop  can  be  obtained  by 
[4] 

using  the  technique  of  the  "Fundamental  Feedback  Equation",  which 

incidentally  demonstrates  the  power  of  this  concept. 


i 


§<S>  ^ t (s)  + t (s)  t (s)  * - - 

J Ol  Cl  os  1 - ipt  (s) 

cs 

"Q  is  the  controlling  variable 

{jf  is  the  controlled  source  (controlled  variable) 


(4.2-1) 


In  Figure  4.1-1  consider  the  independent  input  as  a change  in  A , 
and  the  output  as  change  in  B . The  controlling  variable  ^ and  the 
controlled  source  tf  are  explicitly  marked.  The  several  transmissions 
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of  the  secondary  loop  are  calculated. 


toi (s)  AA(S)  <Ji(s)  = 0 0 


(4.2-2) 


The  "leakage"  transmission  t is  obtained  by  setting  ij>(s)  =0  in 
Equation  4.2-1.  There  is  no  change  in  the  saturation  level  B , so 

toi(s)=°  . 


t . (s)  = c?(s)  ...  = -A  (s) 

ci  J 4»(s)  = 0 me 


(4.2-3) 


A (s)  in  Figure  4.1-1  represents  the  measurement  dynamics  of  A 
me 

from  the  signal  x . The  technique  used  in  the  SOAL  is  a "peak  to  peak" 
detector  which  measures  the  extreme  values  of  x 


As  previously  noted,  this  technique  minimizes  the  measurement  error  due 
to  the  command  and  disturbance  signals.  This  measurement  is  done 


"almost  completely"  within  half  of  the  oscillation  period  T = — , so 

wo 


-.5T0s 

A (s)  = e = -t  . (s) 

me  ci 


(4.2-4) 


Next,  the  transmission  t from  'f  to  "C  , with  i|>(s)  =0  , is 

cs 

obtained  by  finding  the  change  in  the  measured  value  of  A (of  xq  ) 

due  to  a unit  change  in  B . This  transmission  is  composed  of  two 

parts.  The  first  one,  t ^ , represents  the  dynamics  of  the  limit 

cycle  build  up  in  the  primary  loop  (SOAS  case)  . The  second  part,  tcs2 

is  again  the  block  A . Hence 

me 


Sb*81  = tcsl(8)  X tcs2(s) 


W*0  = ~e 


-.5Tqs 


(4.2-5) 


t , (s)  , which  represents  the  dynamics  of  the  SOAS  limit  cycle  buildup, 

csl 

can  be  approximated  (Appendix  IV)  by 
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, , AA , , 

W8*  = Xi(s) 


SOAS 


x 


-J-(S) 


(4.2-6) 


1 + -y-(s) 


This  approximation  suits  the  synthesis  technique  used  and  is  valid 


for  L^(s)  of  type  1.  Equation  4.2-6  essentially  means  that  in  case  of 

I 4 

a change  of  AB  , the  final  value  of  AA  is  AB  — G, G.KP.  ( ju  ) and 

| it  1 2 h o J 

Lf (s)/3 

the  dynamics  of  this  build  up  is  very  close  the  dynamics  of  — — . 

1 + L (s)/3 

So 


(s)  = 

cs 


- G.G_KP,  ( j 0l>  ) 
tt  12  h J o 


Lf(s)/3  -.5Tqs 
— e 


1 + Lf  (s)/3 


(4.2-7) 


Finally  the  transmission  t (s)  which  is  the  transmission  from 

os 


£ to  the  output  AB  is  1. 


t (s)  = 1 
os 

Substituting  Equations  4. 2-2 ,4, 7, 8 in  1 , 


(4.2-8) 


AB.  . 

aX(s) 


-<l>  (s)  e 


-.5Tqs 


1 4 I — • 5TqS  ( Lf(s>/3  \ 

1 + I^Wh'n’h  (l  + Lf-(sV/37  ^S) 


where  the  secondary  loop  transmission  L 


Ls(s)= 


7GlG2Ph^o> 


K e 


-. 5T  s / Lf(s)/3 


/ UflSJ/J  \ 

Vl + Lf (s)/3/ 


<MS) 


(4.2-9) 


(4.2-10) 


In  Equation  4.2-10,  the  maximum  Lg  crossover  frequency  is  less  than 

ljkiQ 

approximately  — , because  of  the  existing  pure  delay  of  . 5Tq  , 

wo 

which  adds  -60°  at  — , and  one  must  allow  say  30°  phase  margin  - 

leaving  badly  needed  90°  for  |l  | reduction. 


It  is  very  interesting  to  see  how  the  uncertainty  in  the  plant  gain 
factor  K has  now  reappeared  in  the  secondary  loop  (after  being 
eliminated  from  the  main  loop),  by  the  presence  of  K in  the  denominator 
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in  Equation  4.2-9.  Hence  this  secondary  loop  must  be  designed  to  be 

stable  over  the  entire  range  of  variation  of  K . We  have  transferred 

the  problem  from  the  main  loop  to  the  secondary  loop.  This  means  that 

the  maximum  crossover  frequency  of  the  secondary  loop  Lg  is  limited 

by  — at  K=K  = K„  . The  maximum  crossover  frequency  at 
3 max  2 ^ 

“o 

K = K . = K.  may  be  much  smaller  than  — , and  it  means  that  at 

mm  1 1 3 

K=Kmin  = Kl  ' the  dynamics  of  the  adjustment  (secondary)  loop  can  be 
quite  slow,  i.e.  at  the  B adjustment  mechanism  is  relatively  slow. 


4.3  Mathematical  Relations 


Suppose  in  Figure  4.1-1  r(t)  = d (t)  = n ( t ) = 0 , K’  is  constant 

in  P(s)  = K 1 P (s)  , the  system  has  achieved  its  steady  state  limit 

h 

cycle  xq  = A sin(a)Qt)  , and  the  secondary  loop  has  also  achieved  its 
steady  state  value  with  the  relay  saturation  level  B'  . It  is  possible 
to  write  for  the  main  loop 
4R  1 

- ^*GlG2K,Ph(*V  = -1  (4-3'1> 

So,  the  loop  transmission  for  the  command  input  and  disturbance  signal, 
if  the  quasi-linear  constraints  are  not  violated,  is: 


L£(s) 


2B ' 
nA 


GlG2K'Ph(s) 


(4.3-2) 


Suppose  now  the  plant  gain  changes  instantaneously  from  K'  to  K"  , 
so  P(s)  =K"P^(s)  . There  is  a change  in  the  amplitude  of  oscillation 

at  x , which  is  sensed  and  corrected  in  a manner  explained  in 
Section  4.2.  After  a certain  transient,  the  relay  achieves  a new  steady 
stale.  Assuming  the  secondary  loop  Lg  is  of  type  1,  xo  (in  the  main 
loop)  ir;  back  to  its  previous  value  =A  sin(w  t)  , with  the  same 
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amplitude  and  frequency  as  before.  In  this  case 


4B" 

■ — • WVJ“o>  - -1 


So  the  forced  loop  transmission  is: 


2 R” 

Lf(s)  = —GlG2K"VS) 


Comparing  Equations  4.3-1  and  3, 
B'K'  = B"K"  = BK  = constant 
TTA 


BK  = 


l4GlG2Ph(^o) 


= constant 


(4.3-3) 

(4.3-4) 

(4 . 3-5a) 
(4 . 3-5b) 


Substituting  in  Equations  4.3-2  and  3,  we  get 
Lf (s)  = Lf (s>  " Lf(s)  = j. 


GlG2Ph(s) 

!GlG2Ph('3“o> 


(4.3-6) 


which  is  similar  to  expression  2.3-9  and  depicts  the  zero  sensitivity 
property  to  pure  gain  changes  of  the  plant. 

We  now  see  the  advantage  of  the  SOAL  over  the  SOAS.  Since  A is 

xr  i xd  i 

constant,  the  quasi-linear  condition  — (t)  <—  , — (t)  <—  is 

A a A a 

independent  of  the  gain  factor  K of  the  plant,  whereas  in  the  SOAS  A 

is  a function  of  K , so  A = A . = A,  , corresponding  to  K=K  . = K, 

mm  1 min  1 

was  used.  Also  in  the  SOAL 

|Co(j“o>  i = |G1(juo) | 

with  A and  | G ^ ( j ooq ) J constants,  so  there  is  no  "most  difficult"  case 
to  satisfy  in  ' ( t ) | f m , whereas  in  the  SOAS  the  "most  difficult" 


case  to  satisfy  was  at  K = K = k 

max  2 


<2 


These  two  constraints  resulted 
ippearing  in  Equation  2.5-37,  and  forcing  w to  be  "large"  in 


the  SOAS  (and  EEAS) . In  the  SOAL,  has  been  eliminated  as  a factor 

K1 

influencing  (f*)0^min  • 
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4.4  SOAL  Sensitivity  to  Uncertainty  in  Plant Dynamics 


If  there  is  uncertainty  in  the  plant  dynamics  (poles  and  zeros) , 
the  loop  transmission  and  the  oscillation  frequency  are  affected 
because 


Lo(j“o>  = iA-WV^o*  = -1 


(4.4-1) 


In  the  SOAL,  A remains  constant.  The  only  change  considered  in  this 
study  is  the  plant  overall  gain  K . A generalization  of  the  synthesis 
method  developed,  which  would  include  dynamic  changes  of  the  plant,  can 
be  developed  using  the  technique  of  Horowitz  and  Sidi  ^ for  linear 
systems. 


4.5  Synthesis  Procedure  for  the  SOAL 


A synthesis  procedure  is  developed  for  the  SOAL,  similar  to  that 


for  the  SOAS. 


4.5-1  Specifications  and  Data 


As  for  the  SOAS  (c.f.  2.5-1), 


Lf (jw) 


Tr(j“)  6 T(jw)  = G3(ju)r7T;M^) 


u)  * Bu) 
o r 


|Td(jw)|=  ^(jW>  = iV^T^y  < |Td  min(jw)|  : ^ 8^  (4 . 5-1) 


P(ju>)  = KPh(jw) 


K.  = K.  S K * K = K 
min  1 2 max 


The  extreme  command  input  ^(t)  -*•  Rg(s) 

The  extreme  disturbance  input  de(t)  -*■  De(s) 

In  addition  to  the  above  (same  as  for  the  SOAS)  specifications,  there 
exists  a complementary  specification  to  be  defined  for  the  SOAL.  The 


following  paragraph  covers  it. 


Dynamics  of  the  Adjustment  Looj 


The  dynamics  of  the  adjustment  loop  change  with  the  value  of  the 


gain  of  the  plant  (Section  4.2).  When  the  plant  has  its  maximum  gain 


the  maximum  possible  crossover  frequency  of  the  adjustment  loop  is 


When  the  plant  has  its  minimum  gain,  the  bandwidth  of  the 


about 


adjustment  loop  decreases  accordingly  to  the  variation  of  K.  and  the 


shaping 


Figure  4.5-1  - Secondary  Loop  Dynamics 


cs  max 


cs  mm 


9 dB/oct 


Assuming  the  loop  has  been  properly  designed  (in  the  frequency  range 


cs  max 


cs  max 


cs  mm 


Then  (see  Equation  4.2-9) 


cs  max 


cs  min' 


cs  max 


cs  mm 


Actually,  one  could  include  in  the  specifications  list  a minimum  dynamic 


1 
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response  for  the  adjustment  loop  (when  K = ),  and  this  might  require 

increasing  (joq  , as  per  Equation  4.5-2.  We  assume  there  is  no  such 
requirement,  and  the  adjustment  loop  should  have  its  possible  crossover 
frequency  maximized. 

4.5-2  The  Constraints 

The  constraints  on  the  SOAL  are  the  same  as  (c.f.  3.5-2)  on  the 
SOAS , i . e . , 


max  x . (t)  5 — 

t 1 1 ' a 


Wbi  * T 


i = r ,d  , a = 3 


i = r,d  , 6 = 3 


(4.5-3) 


If  satisfied,  the  dual-input  describing  functions  can  be  used. 

M, 


"f  ' Nd  ■ T ’ ft 


Also,  as  in  the  SOAS, 


C (t)  5m 

o max 


N = 2N, 
o f 


(4.5-4) 


(4.5-5) 


4.5-3  SOAL  - The  Synthesis  Procedure  (Smooth  Solution) 

The  SOAL  structure  is  pictured  again  in  Figure  4.5-2  with  loop 
transmissions 

(4.5-6) 


f 2B 

Vlu)  ■ WW*-’  ■ Nf  - x - x 


which  is  valid  for  u 5 -tt  for  both  the  forced  and  disturbance  signals 

p 


and 


L (ju>  ) = N G.G.KP.  (jw  ) 
o'o  ol2  hJo 


whic)i  is  valid  only  at  w = w 


M 

o A it  A 


(4.5-7) 


SSL,.  . 
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with  respect  to  the  "extreme"  input  r (t)  , with 

e 

Laplace  transfer  Re(s)  . Substituting  this  value  in  Equation  4.5-10 
and  comparing  with  Equation  4.5-12,  we  get: 


x yrw 

A MfG2KPh(ju)  * aw 


(4.5-13) 


| M^G^KP^  ( jw)  | >,  aw|ReTr  ( jw) 


(4.5-14) 


In  Equation  4.5-14  the  product  M-K  = — — is  invariant  (at  the  steady 

r ttA 

state)  since  A remains  constant,  and  the  product  B.  K (Equation  4.3-5) 
is  also  constant. 

Inequality  4.5-14  assures  non  violation  of  quasi-linearity 
constraints  due  to  cammand  input  signals  for  . An  expression 

similar  to  Equation  4.5-14  is  obtained  for  the  extreme  disturbance  input. 


Lf ( jw) 

ZdW  = ~D(  jw)  ! + L ( jai) 


(4.5-15) 


Zd(jfc>) 

Xd(jW)  = NfG2KPh(jw) 


But  = — , so 


-D(jw)  . Lf(jM) 
NfG2Kph(jw)  1 + Lf  ( j to) 


(4.5-16) 


J»(ju)  = — ' ^1  — 

A MfG2KPh(jw)  1+Lf(jw) 


(4.5-17) 


d , . , 1 1 

A ' J ' aw  3w 


L, (jw) 

IWV*"’  | » . .«  • |De(ju)  | lt  , 


(4.5-18) 


(4.5-19) 


Notice  that  in  Equation  4.5-19,  D(jw)  has  been  replaced  by  its 
extreme  value  D^fjw)  , which  is  the  only  parameter  that  can  assume 
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I 

1 


I 


I 


I 

i 

i 

i 


different  values,  so  its  "most  difficult"  value  has  been  chosen. 

Inequality  4.5-19  being  satisfied  assures  non  violation  of  quasi - 
linearity  constraints  due  to  the  disturbance  input.  There  is  a differ- 
ence in  the  explicitness  of  Equations  4.5-14  and  19,  because  is 

"almost"  completely  defined  (except  possibly  for  more  "far-off"  poles  and 
zeros  which  may  be  added),  while  L^(jw)  is  unknown. 


We  use  the  Nichol ’ s chart  technique  here  just  as  in  the  SOAS , as 
explained  in  Section  2.5-3,  Figures  2. 5-4, 5,  to  plot  the  bounds  on 
due  to  the  disturbance  attenuation  specifications. 


Next  we  will  develop  an  expression  assuring  the  plant  output 

oscillation  level,  to  be  $ the  limit  m . We  have 

lc0(3W0)l  = !M0G2KPh(jwo) I = I 2MfG2KPh ( ^^o5  I * m (4.5-20) 

Again,  there  is  no  "most  difficult"  case  to  satisfy  Equation  4.5-20, 

2BK 

since  M^K  = is  constant  (as  per  Equation  4.3-5). 

Recalling  Equations  4.5-13  and  17,  and  replacing  oj  by  , we  get 


xr 

j 

!T(j“o) 

= 

1 

x . 

T(jV 

= 

ReTr(ju)o) 

I H G KP  ( jw  ) 
z n o 


V^o’ 


MfG2Kph(ja,o) 


because  L^(juQ)  = - — , and 

Equation  4.5-20  into  21  and  22, 


Lf(jwo) 

De{j“o) 

1 + L ( joi  ) 
f o 

MfG2KPh ( jWo) 

(4.5-21) 


(4.5-22) 


1 + Lf (j  w) 


= 1 


Substituting 


X 

r 

ReTr(n> 

A (l<0o) 

= 

MfG2KPh(^o) 

Ki 

|x(l"’o) 

_ 

D (jcu  ) 
c J o 

■ 

MfG2Kph(jMo) 

* - It  R ( jcu  ) I 

m 1 r e J o 1 

* - Id  (jw  ) | 

m e J o 1 


(4.5-23) 


(4.5-24) 
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which  are  non-linear  functions  defining  (w  ) for  "smooth"  X (jw)  , 

3 o min  r J 

K2 

X , ( jto)  . It  is  important  to  note  that  the  factor  — present  in  the 

Q IvT 

SOAS  case  (Equations  2.5-29,30)  has  disappeared.  This  is  a great 

advantage  of  the  SOAL  over  the  SOAS,  but  the  price  paid  is  in  the 

dynamics  of  the  secondary  loop  which  becomes  very  slow  at  K = , 

k2 

especially  for  large  — — 


We  could  also  write: 


Z (jw  ) = R (jw  ) x T (jw  ) 
er  o e J o r J o 


(4.5-25) 


where  the  Zgr(jwo)  represents  the  extreme  plant  output,  due  to  the 


extreme  command  input  rg(t) 


Equation  4.5-23  becomes: 


f(^o) 


5 — Z (jw  ) 
m 1 er  J o 


(4.5-26) 


In  the  same  way 


lZed(ju)o>  I = De(ju)o}  * 


Lf(ja,o) 


1 + Lf (3W0) 


= De(ju.o) 


(4.5-27) 


where  ( jou)  is  the  extreme  plant  output,  due  to  the  extreme  disturb- 


ance d (t) 
e 


Equation  4.5-24  becomes 


d , . J 2 | , . . 

T(3“o)|  4 m I Zed  3uo 


(4.5-28) 


Equations  4.5-26,28  can  be  combined  into 

lx. 


* — Z (jw  ) 
m 1 e J o 


(4.5-29) 

where  Xf  represents  the  forced  signal  at  the  non-linearity  input,  and 
Zg(jw)  is  the  extreme  plant  output,  due  to  either  the  command  or 
disturbance  input.  Since  X^  is  "smooth" 


xf 

1 

2 i , | 

Xf 

— ( ju>) 

- Z jw  J 

— ( jw  ) 

A 

3w 

m 1 e J o' 

A J o 

1 

3w 


(4.5-30) 


< CQ 
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A typical  situation  is  shown  in  Figure  4.5-3,  and  compared  with  the 
SOAS  case  for  the  same  problem  specifications.  The  SOAS  value  is 
w . = 100  rd/sec  > to  _ (SOAL)  = 12  rd/sec 


Figure  4.5-3  - SOAL  - Oscillation  Frequency  Determination 

As  in  the  SOAS,  there  may  be  no  solution  for  some  disturbance 

signals,  but  for  the  same  specifications  and  constraints,  the  level  of 

K2 

disturbances  for  which  a solution  is  available  in  the  SOAL,  is/V  — 

K1 


times  as  great  as  in  the  SOAS.  This  can  be  a significant  factor. 


Tlu-  technique  used  for  definition  of  L^(jco)  and  toQ  is  identical 
to  the  one  defined  for  the  SOAS.  Only  ij>(s)  is  left  to  define,  in  the 
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secondary  loop.  As  noted  earlier,  iMs)  is  chosen  to  achieve  a 
satisfactory  secondary  loop  L , with  maximum  possible  crossover 

b 

frequency. 


4 . 6 Numerical  Example  - SOAL  - Smooth  Solution 


In  this  numerical  example  we  use  the  same  specifications  as  for  the 
numerical  example  used  for  the  SOAS,  Section  2.6,  but  with  one  differ- 
ence - mg  .01  , which  is  l/100th  of  the  value  used  in  the  SOAS,  where 

K1 

— is  also  equal  to  1/100.  In  this  way  the  same  functions  G,  (s)  , 

K 2 1 

G2(s)  and  G^(s)  appear  in  the  outer  loop  and  prefilter.  B , the 

saturation  value  for  the  ideal  relay,  is  a variable  with  B = B = 16  it 

max 

(as  in  the  SOAS) , corresponding  to  the  case  where  the  plant  has  its 

minimum  gain  K = K . =K=1  , (A  = .5V),  and  with 

min  1 


K1 

B = Bm . „ = B0  = x — = .16  it  , corresponding  to  the  case  where  the  plant 


min 


has  its  maximum  gain  K=K  = K„  = 100  ( A is  still  .5  V ) . iMs)  is 

nisx  2 

. 4 

chosen  = — , so  the  crossover  frequency  of  Lg  is  maximized  when  the 


plant  has  its  maximum  gain  K = 

- . 5T  s 1 1 

A = e ° , where  T = — = — — sec 

me  o a)  4 . 6 


In  this  design  the  block  defining 


o 

Computer  simulations  were  performed  on  the  CDC  6400  system  at 
University  of  Colorado,  using  "Mimic"  program.  Two  types  of  runs  were 
performed.  First,  the  gain  of  the  plant  is  changed  (with  different 
initial  gain  values)  and  no  command  or  disturbance  inputs  are  considered. 
The  dynamics  of  the  adjustment  loop  changes  accordingly.  The  gain  is 
changed  abruptly  from  to  . The  table  below  gives  the  list 

of  the  simulations  performed. 


ST***?* ' 


Table  1 - Simulations  performed  - Gain  changes 


In  these  figures  only  two  variables  are  shown;  sub  a depicts  x (t) 

o 

(the  input  to  the  non-linearity)  and  the  plant  ga~  K ; sub  b depicts 

B (the  non-linearity  saturation  level)  and  K . Note  in  Figure  4.6-1 

that  when  K.  =100  , K, . =50  , the  new  steady  state  is  achieved  after 

in  fi 

about  5 oscillations,  while  in  Figure  4.6-4,  K.  = 2 , K. . =1  , the 

in  f i 

dynamics  of  adjustment  are  very  slow. 

The  second  group  of  runs  shows  the  system  response  to  command  and 
disturbance  inputs,  when  extreme  values  are  applied  for  fixed  plant  gain 
values.  The  table  below  summarized  the  simulations  done. 


Figure 

K 

R 

D 

4.6-5 

100 

R 

e 

0 

4.6-6 

10 

R 

e 

0 

4.6-7 

1 

R 

e 

0 

4.6-8 

100 

0 

D 

e 

4.6-9 

10 

0 

D 

e 

4.6-10 

1 

0 

D 

e 

Table  2 - Simulations  performed. 

System  response  to  command  and  disturbances 


As  before,  in  each  graph  only  two  variables  are  shown;  sub  a 
represents  the  input  and  the  system  response  and  sub  b represents  the 
input  to  the  system  and  the  input  to  the  non-linearity. 
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It  is  worth  noting  the  differences  between  the  SOAL  and  the  final 

Minneapolis-Honeywell  (MH)  adaptive  flight  control  2-loop  design  used  in 
[211 

the  X-15  . The  basic  philosophy  of  the  MH  and  SOAL  designs  is 

identical,  i.e.  to  use  the  oscillating  component  amplitude  measurement 
in  order  to  change  the  gain  in  the  main  loop,  such  that  in  the  steady- 
state,  the  oscillating  amplitude  is  constant.  However,  the  means  of 
execution  are  radically  different. 

In  the  MH  design,  the  oscillation  is  measured  near  or  at  the  system 

output,  where  it  is  not  large  relative  to  the  other  signal  components. 

A band-pass  filter  centered  at  is  therefore  needed,  introducing  an 

[27] 

additional  very  significant  lag  in  the  second  loop  . Other  signals, 

such  as  command  and  especially  fast  disturbance  signals,  also  excite  the 

band-pass  filter.  If  its  bandwidth  is  therefore  decreased,  its  effective 

time  constant  is  increased  since  the  latter  is  inversely  proportional  to 
[27  ] 

the  bandwidth  . A disturbance  signal  may  thus  conceivably  be  inter- 

preted as  a plant  gain  change.  In  addition,  the  narrower  the  band  of 
the  filter,  the  more  sensitive  it  is  to  changes  in  the  oscillating 
frequency. 

Finally,  the  output  of  the  band-pass  amplifier  must  be  measured  and 
compared  with  a reference  - exactly  as  in  the  SOAL.  Thus,  the  band  pass 
filter  has  been  eliminated  in  the  SOAL.  In  the  SOAL,  the  measurement  of 
the  oscillation  is  made  at  a point  where  it  is  inherently  the  dominant 
signal.  Peak  to  peak  measurement  decreases  the  sensitivity  to  other 
signals,  providing  they  are  relatively  constant  over  half  a period. 
Finally,  it  is  important  to  note  that  no  analytic  design  theory  was  made 
available  by  the  Minneapolis-Honeywell  effort. 
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CHAPTER  FIVE 


S.O.A.N.L. 


(SOAL  with  non-linearity  in  the  secondary  loop) 


5.1  General 


Chapter  4 dealt  with  the  SOAL  structure,  which  has  some  important 

advantages  over  the  SOAS  structure.  The  SOAL  succeeded  in  eliminating 

K K 

the  factor  — from  the  outer  (adaptive)  loop.  But  the  factor 


K . 
min 


K . 
min 


reappeared  in  the  secondary  loop  which  determines  the  rate  of  adaptation 
of  the  system.  The  latter  is  given  by  Equation  4.2-9,  repeated  here. 

-.  5T  s 


AB 

AA 


(s)  = 


-Ms)  e 


1 + 


r G G,P,  (jio  ) 
tt  1 2 h o 


K e 


-.5Tqs  / Lf(s)/3  \ 

\1 * Lf (s) /3/ 


(5.1-1) 


^ (s) 


The  idea  of  the  SOANL  is  to  keep  the  existing  advantages  of  the 

SOAL,  i.e.,  elimination  of  the  factor  ( ) in  the  outer  loop,  but  to 

' min  ' 

eliminate  it  also  from  the  secondary  loop  by  eliminating  K in 

Equation  5.1-1.  This  could  be  achieved  by  introducing  a multiplier  in 

the  secondary  loop,  with  multiplication  factor  inversely  proportional  to 

K . But  K is  not  available  directly,  so  B^  , which  is  related  to 

K , can  be  used,  since  in  the  steady  state  case  KB  is  constant 

(Equation  4.3-5).  In  the  transient  situation  the  information  about  K 

existing  at  the  Ame  block  output  could  also  be  used  to  identify  K , 

because  if  B has  not  achieved  its  steady  state,  then  A *A 

m r 

(Figures  4.1-1,  5.1-1).  In  this  study  this  transient  information  has 
not  been  used,  but  can  be  considered  in  future  research. 


Figure  5.1-1  represents  the  SOANL  structure  which,  if  compared  with 


the  SOAL  in  Figure  4.1-1,  includes  a multiplier  which  multiplies 


Figure  5.1-1  - SOANL  Structure 


The  multiplier  changes  the  dynamics  of  the  secondary  loop,  which 


will  be  analyzed  in  Section  5.2.  Just  as  in  the  SOAL,  the  main  loop  is 


5 . 2 The  Stability  and  Dynamics  of  the  Secondary  Loo] 


SOANL 


As  in  Section  4.2,  the  fundamental  feedback  equation  will  be  used 


leading  to 


where  ^'(s)  is  the  equivalent  relation  between 


All  the  other  terms  have  the  same  meanings  as  in  Section  4.2 
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In  Figure  5.1-1  the  independent  input  is  considered  as  a change  in 
A , and  the  output  as  a change  in  B . The  several  transmissions  are: 


‘oi181  ■ I(s> 


(s)  = 0 


= 0 


t 


^(s)  -£(.) 


' (s)  = 0 _Ame(s)  6 


- . 5T  s 


(5.2-2) 

(5.2-3) 


We  can  also  write  (see  Figure  5.1-1) 
[e  (t)  B (t)  ] 


B (t)  = p(t)  * 


-L-1 


B . 
min 


where  y =JU  ijj  . Assuming  B(t)  does  not  change  "very  much"  from  its 


initial  value  B.  , 

l 

B (t)  = B + ABi(t) 


so 


so 


B (t)  « 


p(t)  * [e (t)  • Bi] 
mm 


B (s  ) = ip  (s ) x E(s) 


(i^) 

' min  ' 


Then 


mm 

tSG(s)  = I(S)  = (S)  = *(s) 


(i^) 

' mm  ' 


(5.2-4) 


(5.2-5) 


(5.2-6) 


t (s)  and  t (s)  are  identical  to  the  functions  found  for  the  SOAL. 
cs  os 


fccs(s)  lGlG2KPh<J“o,l  1 + L (s) /3 


t (s)  = 1 
os 


L (s)/3  -,5T  s 

— ° 


(5.2-7) 

(5.2-8) 


Substituting  Equations  5. 2-2, 3,6, 7,8  in  1, 

B. 


AB,  . 

-rr(s) 

AA 


' mm  ' 


• <Ms)  e 


- . 5T  s 
o 


1 + 


(5.2-9) 


KB.  -,5T  s 

7GlG2Ph(i“o)  * B~~'  6 ° * ^(s)  * IT 

mm 


Lf (s)/3 
Lf (sV/3 
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In  Equation  5.2-9  , the  initial  value  of  B(t)  , is  such  that 

if  the  steady  state  was  achieved  (Equation  4.3-5),  B^K  = constant  , and 
the  denominator  is  independent  of  the  gain  of  the  plant,  so  the  dynamics 
of  the  secondary  loop  are  independent  of  the  plant  gain.  The  assumption 
that  B (t)  does  not  change  "very  much"  once  the  steady  state  is  achieved, 
is  valid,  considering  the  dynamics  of  the  adjustment  loop  being  faster 
than  the  changes  in  K . More  than  that,  simulations  have  shown  that 
even  considering  bigger  changes  in  K and  consequently  in  B , the 
results  are  still  vlid  for  Equation  5.2-9. 

When  K = K (in  Equation  5.2-9),  B equals  B . , so 

max  ^ mm 

Equation  5.2-9  is  identical  to  the  equivalent  SOAL  equation  (4.2-9). 

This  means  that  for  K = K the  SOANL  has  the  same  dynamics  for  the 

max 

secondary  loop  as  the  SOAL.  Suppose  now  K = » so 

K 

B = B = B . x , and  again  the  denominator  of  Equation  5.2-9  has 

max  mm  K . 

mm 

the  same  function  as  before  ( K-K  . ),  and  the  dynamics  of  the  adjust- 

mm 

ment  loop  are  constant. 

SOANL  Sensitivity  to  Dynamic  Changes 

Section  4.4  for  the  SOAL  is  valid  for  the  SOANL. 


5 . 3 Synthesis  Procedure  and  Numerical  Example  (Smooth  Solution)  - SOANL 

The  synthesis  procedure,  the  specifications,  data  and  constraints 
are  the  same  as  for  the  SOAL,  Section  4.5. 

For  the  main  loop,  the  SOANL  is  identical  to  the  SOAL  both  for 
command  and  disturbance  signals.  The  advantage  is  in  the  secondary 
loop,  where  the  dynamics  of  adaptation  are  now  almost  constant. 
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The  same  specifications  of  the  SOAL  numerical  example  are  used  in 
the  SOANL  numerical  example.  As  a consequence,  all  the  SOANL  blocks  have 
the  same  functions  as  the  SOAL  blocks,  including  i|>(s)  • The  only  change 

is  the  multiplier  in  the  secondary  loop. 

A set  of  runs  were  made  at  the  CDC  6400  system  - C.U.,  using  Mimic 
language.  Two  types  of  runs  were  performed:  First,  changes  in  gain  of 
the  plant,  showing  the  dynamics  of  the  adjustment  loop.  Next,  the  system 
response  to  extreme  command  and  disturbance  inputs  for  several  plant 
gains . 

The  results  show  that  for  the  same  proportional  plant  gain  change, 
the  adjustment  loop  dynamics  are  constant  and  independent  of  K , which 
is  the  SOANL  main  advantage  if  compared  with  the  SOAL.  The  system 
response  for  both  the  extreme  command  and  disturbance  inputs  is  similar 
to  the  SOAL  response,  as  expected. 

Some  details  should  be  considered,  to  avoid  "pitfalls"  in  the  SOANL 
implementation.  Figure  5.4-1  represents  the  part  of  the  secondary  loop 
involving  the  multiplier. 


Figure  5.4-1  - Details  in  the  Secondary  Loop 


/**■*’. 
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I 


! 

i 


First,  if  for  any  reason  (due  to  dynamic  overshoot) , B goes  to 

zero,  both  loops  are  open.  A minimum  expected  value  for  B , based  on 

the  maximum  of  K , should  be  established,  and  if  necessary,  non-linear 

circuits  used  to  avoid  passing  this  minimum  value.  Second,  it  is 

advisable  to  have  an  integrator  in  iMs)  , so  there  is  no  steady  error 

in  A , the  limit  cycle  amplitude  at  x .It  is  also  recommended  to 

o 

avoid  pure  gains  in  <Ms)  . As  an  example,  suppose  we  choose 
Pi 

iMs)  = — + P2 


so  in  Figure  5.4-1 

s(t)  xf-^-xP2  + /fc  |-^-*e(t)dt  = B (t) 
min  0 min 

B . 

If  at  t^  , it  so  happens  that  e(t^)  =— — , then  the  above  equation 


becomes 

t 


1 B(t) 

J B . 

0 min 


x e (t) dt  = 0 


which  may  not  be  true,  i.e.,  there  is  no  solution  for  B(t) 


The  simulation  results  are  presented  in  the  next  figures.  First 
the  plant  gain  is  changed  and  no  input  applied.  This  is  shown  in 
Figures  5.4-2  - 10.  For  the  same  relative  plant  gain  change  the  adjust- 
ment dynamics  are  identical.  In  Figure  5.4-10  when  K , the  plant  gain, 
changes  instantaneously  from  its  maximum  value  (100)  to  its  minimum  value 
(1) , the  oscillations  are  temporarily  quenched. 

Next,  the  system  response  to  the  extreme  command  input  is  shown, 
for  several  plant  gains  (Figures  5.4-11  - 13).  Finally,  the  system 
response  to  the  extreme  disturbance  input  for  several  plant  gains  is 


represented  in  Figures  5.4-14  - 16. 


Secondary  Loop  Dynamics 


Figure  5.4-8b  SOANL  System  Response  Figure  5.4-8a 

Relay  Saturation  Level  Secondary  Loop  Dynamics  x 18  - Non-linearity  Input 
Plant  Gain  K - Plant  Gain 


I - Non-linearity  Input  K = 100  XC  - System  Output 

Applied  Command  R _ Appiied  Command 
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CHAPTER  SIX 
E .E .A.L. 

(EEAS  with  a secondary  loop 
controlling  the  non-linear  adaptive  saturating  level) 


6.1  General 


In  Chapter  3 it  was  shown  that  the  EEAS  has  some  advantages  and 
greater  flexibility  that  the  SOAS.  The  advantages  due  to  the  fact  that 
Lo(ja)Q)  need  not  be  equal  to  -1  are  in  the  bandwidth  of  the  loop  trans- 
mission and  on  the  maximum  disturbance  signals  levels,  which  do  not 


violate  the  quasi-linearity  constraints.  But  the  factor 


Kmax 

K . 
min 


still 


appears  in  the  expression  determining  these  properties  (Equations  3.5-30 
and  33) , and  still  there  exist  severe  limitations  on  the  possible  (less 
severe  than  for  the  SOAS)  disturbances,  if  quasi-linearity  conditions 
are  not  to  be  violated  (Equation  3.5-33) . 

In  Chapter  4 it  was  shown  that  there  is  no  "advantage"  in  the  SOAL 
flexibility  if  compared  with  the  SOAS,  since  Lo(jcoo)  = -1  , but  there 

is  an  advantage  in  the  system  loop  transmission  and/or  in  the  disturb- 


ances limitations  because 


max 

K . 
min 


disappears  from  the  expression  defining 


these  parameters  (Equations  4.5-23,24).  The  idea  of  the  EEAL  is  to  use 
the  advantages  of  the  SOAL  and  to  widen  the  EEAS  flexibility.  Several 
possible  structures  are  reviewed  in  Section  6.2,  and  the  specific  one 
shown  in  Fiqure  6.2-2  is  chosen.  This  structure  seems  to  be  the  one 
which  can  meet  most  of  the  above  EEAL  requirements. 


6.2  Possible  EEAL  Structures 


Figure  6.2-1  - Structure  Number  1 - EEAS 


In  this  structure,  special  care  should  be  taken  to  allow  a possible 
solution  (with  positive  B ) , for  all  values  of  the  gain  of  the  plant. 
Assuming  such  a solution  exists,  and  by  a similar  development  to  the  one 
used  in  Section  3.5,  the  following  equations  are  obtained: 


X . R T 

r , . . .1  e r , . 

— <3U  ) i - — (30)  ) 

A o m o 


I—  (iio  ) I a - Id  T.(jw  ) 

e d o 


(6.2-1) 


(6.2-2) 


which  are  similar  expressions  to  Equations  3.5-30,33,36,  where  the  factor 


have  been  achieved  in  this  structure,  but  no  more  flexibility  than  the 


EEAL  - Structure  Number  2 


Structure  number  2 is  represented  in  Figure  6.2-2,  and  is  the  one 


chosen  to  be  analyzed  with  more  detail 


Structure  Number  2 - EEAL 


The  main  difference  between  structures  number  2 and  1 (or  the  EEAS) 


is  the  place  where  the  external  signal  is  injected.  In  Figure  2 the 


output  of  G. (s)  ( v(t)  ) is  not  the  input  x(t)  to  the  non-linearity 


Quasi-linearity  constraints  must  be  satisfied  by  x(t)  and  not 
necessarily  by  v(t) 

In  this  structure  the  amplitude  V of  the  oscillating  signal  vo(t) 

is  kept  (or  at  least  tried  to  be  kept)  constant.  This  is  done  by  the 

"secondary"  loop  of  this  structure.  Figure  6.2-2.  Vme  measures  the 

amplitude  of  oscillation  of  VQ(t)  . The  result  of  this  measurement 

V is  compared  with  a reference  value  V . The  difference  (error) 
m r 

passes  through  a compensation  network  (s)  , whose  output  controls  B , 

which  affects  the  amplitude  V of  v (t)  . This  "secondary"  loop  is 

o 

not  an  "instantaneous"  loop  and  its  dynamics  and  stability  problems  must 
be  analyzed.  This  is  done  in  Section  6.5  in  a similar  way  to  the 
analysis  performed  for  the  SOAL  secondary  loop,  Section  4.2. 

Suppose  in  Figure  2 the  system  has  achieved  its  steady  state 
condition,  and  r(t)  =d(t)  = n (t)  =0  . Then 


F(t)  = A sin  a)  t , x (t)  = A sin(w  t+(p) 
o o o o 


v (t)  = V sin(w  t+cp+0) 
o o 


with  their  respective  phasors 

(ju>  t)  j (w  t +<p) 

I(jwo)  = Aq  e ° , X(ja)Q)  = A e 


(6.2-4) 


Also : 


I(juo) 

X(3<1)o)  ‘ 1 + Lo(ja)o) 


V*V  =TG1G2KV*V 


Combining  Equations  6. 2-5, 6, 7, 


A eDiP  + M G,G.KP.  (ju  ) e3tP  = A 
o 1 2 h J o c 


| A t KMoGlG2Ph(ja»0)|  = Aq 


(6.2-5) 


(6.2-6) 


(6.2-7) 


(6.2-8) 


(6.2-9) 
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Since  A is  constant  if  KM  was  kept  constant,  A remains  constant 
o o 

(as  K varies) . As  explained  before,  V is  kept  constant  (in  the 
steady  state)  by  the  secondary  loop. 


IV(ja>  ) I = MqK  |G1G2Ph  (jo)Q)  | = constant 


So 


| V ( j u>o ) 


— BK  = 2M  K = M K = T — _ , . , 

n f o lGiG2Ph(3“o> 


= constant 


(6.2-10a) 


(6.2-10b) 


The  factors  BK  , M.K  , M K have  constant  values.  The  loop  trans- 
f o 

mission  for  the  forced  and  disturbance  signals  has  the  property  of  zero 
sensitivity  to  pure  gain  changes,  since 
M K MfK 

= IT  GlG2Ph(jW)  = ~pr  GlG2Ph(jU>)  (6-2'n) 

is  independent  of  K 

Structure  number  2 seems  to  be  very  convenient  because  one  of  the 
links  existing  in  all  systems  studied  (SOAS,  EEAS,  SOAL,  SOANL)  , given  by 
A “ co|G1(ju»0)  | (6.2-12) 

is  broken.  A and  Cq  are  the  oscillation  components  at  the  non- 
linearity input  and  at  the  system  output.  In  the  EEAL  the  input  A to 

the  non-linearity  is  the  sum  of  and  of  V =C  |g, (jw  ) I 

001.0 

The  "real"  information  about  the  plant  gain  is  in  V ; that  is  the 
reason  why  v was  the  chosen  variable  from  which  the  amplitude  of 
oscillation  is  measured  ( x had  the  advantage  that  the  sinusoid  was 
easily  the  biggest  part  of  the  signal  and  peak  to  peak  measurement 
wiped  out  Xj.  ) . 

These  facts  are  also  part  of  the  EEAL  problems  and  limitations. 

If  the  injected  signal  3 = Aq  sin(a)Qt)  has  a much  larger  amplitude 

than  v = V sin(w  t+tp+9)  , |a  I >>  ! V ! , which  is  a common  case  for 

o 1 o 


sa 


-183- 


Figure  6.2-3  - Structure  Number  3 - EEAS 


At  x quasi-linearity  conditions  are  satisfied,  so  the  same  peak 
to  peak  measurement  system  used  for  the  SOAL  can  be  also  used.  The 
problem  in  this  structure  is  the  high  sensitivity  of  the  plant  identi- 


the  EEAL  (see  numerical  example,  Section  6.7) , then 

v(t)  =vr(t)  tv^Ct)  +vQ(t)  , does  not  satisfy  quasi-linearity  constraints 

and  |v  (t)  | or  |v.  (t)  | could  be  larger  than  |v  (t) | =V  , making 
it  d o inflx 

the  oscillation  amplitude  measurement  a difficult  process,  probably  with 
added  dynamic  lags.  This  structure,  due  to  the  non-existing  link 
expressed  by  Equation  6.2-12,  gives  the  expected  flexibility  for  the 
EEAL. 

EEAL  - Structure  Number  3 

Structure  number  3 is  similar  to  structure  number  2,  in  almost  all 
aspects.  The  only  difference  is  in  the  variable  chosen  to  measure  the 
oscillation  component.  Figure  6.2-3  represents  structure  number  3. 


F/B  9/A 
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(AQj  s 60  V , as  in  the  EEAL  numerical  example.  If  there  exists  an 
error  of  1.6  % in  the  measurement  of  A from  x(t)  , the  secondary  loop 
will  force  V to  have  an  error  of  at  least  100  %,  by  an  improper  B 
adjustment.  The  plant  gain  identification  is  incorrect.  In  case  XQ(t) 
and  vQ(t)  are  not  in  phase  (0°  or  180°)  the  sensitivity  can  be  even 
larger,  as  depicted  in  Figure  6.2-4  ( V is  the  amplitude  of  VQ(t)  > 
and  A the  amplitude  of  xQ(t)  )•  For  this  reason  structure  number  3 
is  not  chosen. 


Figure  6.2-3  - Small  Influence  of  V in  Plant  Identification 

( V1«2V1  , A2«A1  ) 


6 . 3 Mathematical  Relations 

Suppose  in  Figure  6.2-2,  the  chosen  structure,  r(t)  =d(t)  = h (t)  =0, 
the  system  is  in  its  steady  state,  the  linear  functions  and  the  injected 
signal  defined  and  known. 


First  the  oscillating  input  to  the  non-linearity  is  calculated. 
Based  on  Equation  6.2-6, 

K ju  ) 


X(jW0}  1 + L (jui  ) 


(6.3-la) 


or 


rt 

pT 


1 + L ( ju)  ) 
1 o J o 


(6.3-lb) 


where 


V><V  = T-Wh'^o’ 


(6.3-2) 


We  repeat  here  Equation  6.2-8 

A ejiP  + M G.G.KP,  (jo)  ) ej4>  = A 
o 1 2 h J o c 


o o 


(6.3-3) 

(6.3-4) 


A graphical  solution  for  Equation  6.3-4  is  provided  in  Figure  6.3-1 
giving  the  values  of  A and  g>  once  A , M K and  G.G_P.  (j<o  ) are 

O O 1 Z xl  O 


defined.  The  product  MqK  has  been  considered  as  a single  term  since. 


due  to  the  secondary  loop  (Equation  6.2-10),  it  has  a constant  value. 


Figure  6.3-1  - Graphical  Solution  - Equation  6.3-4 
In  Figure  6.3-1,  the  vector  OA  represents 


* ‘ A lGlG2Ph(J“o> 
o 


(6.3-5) 


where 


= e 


(6.3-6) 


O' 


-186- 


— , which  gives  the  solution  for  A , is  a real  positive  number  and 

o 

is  obtained  by  tracing  a parallel  line  to  the  real  axis,  with  origin  in 
A until  B , intersection  of  this  line  with  the  unit  radius  circle. 


(6.3-7) 


■+  — TCP 

OB  represents  e , which  gives  the  solution  for  the  phase  angle  ip  . 


-*■  -id) 
OB  = e 


(6.3-8) 


As  in  Section  3.3  for  the  EEAS , sometimes  there  are  "two"  solutions, 
and  sometimes  there  is  "no"  solution  for  Equation  6.3-4.  The  discussion 
of  the  possible  solutions  as  well  as  its  physical  interpretation  are 
similar  to  the  results  obtained  in  Section  3.3,  and  are  not  going  to  be 

Mo  3tt 

reviewed  here.  As  an  example  if  — |g..G  KP  (jw  )|  , and  — S0S  2 , 

A l 4 n o <c 

o 

as  represented  by  the  vector  OC  in  Figure  6.3-1,  there  is  no  possible 
real  positive  solution  for  A (and  ip  ) , so  the  system  self-oscillates 
at  w 


6.4  EEAL  Sensitivity  to  Change  in  Plant  Dynamics 


It  was  shown  that  the  EEAL  (the  chosen  structure)  has  zero  sensitiv- 
ity to  pure  gain  changes  of  the  plant  (Equation  6.2-11).  It  was  also 
shown  that  for  the  SOAS,  SOAL  and  SOANL , the  oscillation  frequency 
changes  if  the  plant  changes  its  dynamic  characteristics.  For  the  EEAS 
(Section  3.4-3),  since  the  oscillation  frequency  is  defined  by  the 
injected  signal,  there  is  no  change  in  the  oscillation  frequency,  but 


there  is  a change  in  the  L ( j u>  ) position  in  the  Nichol’s  chart,  as 

o o 

[x>r  Section  3.4-3.  In  this  research  the  only  case  studied  is  the  one 
where  the  plant  changes  its  overall  gain,  but  it  is  interesting  to  see 


L8' 


what  happens  to  Lo(jwo)  position  when  there  are  changes  in  the  poles 


vj<v 


and  zeros  of  the  plant.  In  Figure  6.2-2, 

V(juo)  Vj<V 

I(j<oo)  = 1+L0(3«0)  (6,4_1) 

If  there  is  a change  in  the  dynamics  of  the  plant  and  the  system  is 
oscillating  only  at  wo  , the  secondary  loop,  after  the  transient 
period,  assures  that  | V ( joo^) | returns  to  its  pervious  constant 

L0(j“0) 

amplitude,  so  in  Equation  6.4-1,  — — — = constant  , is  defined 

1 + L (jo)  ) 
o o 

in  the  Nichol's  chart  as  one  of  the  constant  closed  loop  contours.  Thus 
as  the  plant  changes  its  phase  angle  J P ( juiQ)  , Lo(ju)Q)  follows  this 

change  in  open  loop,  along  the  constant  closed  loop  contour  defined 
previously  for  the  EEAS , in  Section  3.4-3. 

The  idea  suggested  in  that  section,  for  "zero  sensitivity"  to 
changes  in  plant  dynamics,  is  valid  for  the  EEAL,  by  using  a dynamic 
non-linearity  instead  of  a static  one. 

6.5  Dynamics  and  Stability  of  the  Secondary  Loop  - EEAL 

In  Figure  6.2-2  the  function  of  the  secondary  loop  containing  i|>(s) 
is  to  adjust  the  saturation  level  B , so  that  the  amplitude  of  vq  is 
kept  constant  (in  the  steady  state) . The  dynamics  of  this  secondary 

loop  are  studied  here,  using  the  same  technique  as  for  the  SOAL, 

, [21 
Section  4.2,  based  on  the  Fundamental  Feedback  Equation 


^(s)  = t . (s)  + t . (s)  • t (s)  • t — r 
■f  . oi  ci  os  1 - <I>t  (s) 

cs 


(6.5-1) 


In  Figure  6.2-2  consider  the  independent  input  as  a change  in  V , 
and  the  output  as  a change  in  B . S and  ^ are  marked  in  this 


figure.  The  transmissions  are: 


AB 

fcoi(s)  = AV  ip  ( s)  =0  " 0 


There  is  no  "leakage"  transmission  in  the  loop. 


t . (s) 
ci 


■ ^ (s)  j'Hs)  = 0 = AV(S)  Ij/(s)  = 0 vme(s) 


(6.5-2) 


(6.5-3) 


V (s)  is  the  block  representing  the  measurement  characteristics 
me 

of  the  amplitude  of  oscillation  VQ(t)  from  v(t)  . The  technique  used 
for  the  EEAL  is  more  complicated  than  the  "peak  to  peak"  technique  used 
in  the  SOAL,  since  quasi-linearity  conditions  are  not  satisfied  at  v(t)  . 
The  technique  used  assumes  that  vf(t)  ( v(t)  = vq (t)  +vf(t)  ),  is  linear 
in  t for  each  cycle; (in  Figure  6-5-1,  xf(t)  is  considered  as  a 
straight  line  in  the  interval  of  time  t^  - tj  ) . 


t)  = XQ(t)  + (t) 


y w 2°  \y/ 


Figure  6.5-1  - Oscillation-Amplitude  Measurement  - Linear  Approximation 


-189- 


This  xf(t)  "linear  approximation"  for  each  cycle  seems  to  satisfy 
most  of  the  cases  (see  numerical  example  simulations) . A "parabolic 
approximation"  can  be  also  used  with  added  time  delay  in  the  measurement 
characteristics  and  with  better  static  accuracy. 


We  can  write  (Figure  6.5-1): 
v(t^)  = a + bt^  + V 
v(t2)  = a + bt2  - V 
v(t3)  = a + bt3  + V 

But  t2  - t^  = t3  - t2  = Tq/2  , where  Tq  is  the  oscillation  period.  So 

v(tx)  - 2v(t2)  + v(t3)  = 4V 
or 

V = i [vtt^  - 2v(t2)  + v(t3)J  (6.5-4) 


since  it  takes  one  period  to  almost  complete  the  measurement 

-sT 

t . (s)  = V « e ° 
ci  me 


(6.5-5) 


Next,  ip  ( s ) is  the  only  degree  of  freedom  in  the  design  of  the  secondary 

loop  (assuming  the  main  loop  has  been  already  designed  to  cope  with 

desired  specifications  and  existing  constraints).  Then  tcg(s)  , the 

transmission  from  ^ to  ^ , with  4*(s)  removed,  is  obtained  by  finding 

the  change  in  the  measured  value  of  V (of  v (t)  ) due  to  a unit 

o 

change  in  B . This  transmission  has  two  parts.  The  first,  tcsl  > 
represents  the  dynamics  of  the  limit  cycle  build-up,  without  the 
secondary  loop.  The  second  part,  fcCS2  • a9ain  the  block  Vme  , 
already  defined.  So 


and 


"cs (s)  = fccsl  * tcs2<8) 


-T  s 

tcs2(s)  5 6 ° 


(6.  ^ ‘6) 


(6.5-7) 
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I 

I 

i 

i 

! 


S 


In  this  case  tcsl(s)  is  easier  to  calculate  than  in  the  SOAL, 
since  the  input  to  the  non-linearity  is  "almost"  defined  by  the  injected 
signal.  This  is  due  to  the  tendency  in  the  EEAL  optimization  design  to 
have  |^0 (fc) |max>>|Vo(t) |max  ' as  exPlained  in  the  synthesis  procedure 
(Section  6.6) . In  this  case 


V (s)  = I (s)  * g.g.kp,  (s) 
o O A 7T  1 2 h 

o 

But  | I (s)  | «A  (as  mentioned  above)  , so 


"csl^ 


(6.5-8) 


The  last  transmission  t (s)  , from  to  the  output  &B £ 1 

os 


i 


Substituting  the  above  in  Equation  6.5-1, 


AB 

Z7(s) 


~4>{s)  e 


-T  s 
o 


1 * f GlG2Ph  *<=)  K e V 


(6.5-9) 


In  Equation  6.5-9  the  maximum  possible  loop  transmission  is  about  — , 

2tt 

because  of  the  inherent  delay  of  — =Tq  < in  the  measurement  process. 


o 

Again  it  is  interesting  to  see  how  "mother  nature"  reintroduces  the 
plant  gain  uncertainty  into  the  secondary  loop,  by  the  presence  of  K 
in  the  denominator  of  Equation  6.5-9.  Hence,  as  in  the  SOAL,  the 
secondary  loop  must  be  designed  to  be  stable  over  the  entire  region  of 
variation  of  K . At  large  K the  dynamics  of  the  adjustment  loop 
will  be  relatively  fast,  and  at  K = K ^ , the  dynamics  of  this  adjust- 

ment loop  will  be  relatively  slow.  A comment  regarding  the  "stability" 
of  the  main  loop.  It  has  been  noticed  that  for  the  EEAS  (and  EEAL)  it 
is  necessary  to  have  a minimum  gain  margin  to  avoid  self-oscillation. 

The  gain  margin  is  defined  by  L,.(ju  ) =^t!g,G„KP,  (jw  )|  . But  A is 

r tt  ttA  1 / h it  1 

the  amplitude  of  x()(t)  , the  sum  of  the  constant  injected  signal  ^(t) 
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and  v (t)  . If  there  occurs  an  abrupt  change  in  the  gain  of  the  plant 

o 

( K increases  its  value),  vQ(t)  changes  its  amplitude  but  xQ(t) 

(and  A ) does  not  follow  the  same  proportional  change.  While  the 
secondary  loop  has  not  achieved  the  steady  state  condition,  the  loop 
Lf(s)  increases  (or  decreases)  its  gain,  not  being  zero  sensitive  in 
the  transient  period.  It  is  therefore  recommended  that  an  extra  gain 
margin  to  the  existing  specification  for  the  EEAS  - 2.56  dB 
(Equations  3.2-8,  3.5-21)  , should  be  considered. 

6.6  Synthesis  Procedure  for  the  EEAL 

A synthesis  procedure  is  developed  for  the  EEAL.  Assumptions  are 
the  same  as  in  Section  3.5  for  the  EEAS. 

6.6-1  Specifications  and  Data 

It  is  assumed  that  the  set  of  specifications  and  data  listed  below 
are  known  or  can  be  estimated. 

- Desired  Command  Input  Response 

- Disturbance  Attenuation 

- The  Plant 

- Extreme  Command  Input 

- Extreme  Disturbance  Input 

The  meanings  and  comments  given  in  Sections  2.5-1  and  3.5-1  are  valid 
here. 

If  there  are  "abrupt"  changes  of  gain  (expected) , this  should  also 
be  specified.  The  reason  for  this  new  specification,  explained  in 
Section  6.5,  is  related  to  the  minimum  gair  margin  condition,  which 


guarantees  no  self-oscillation  in  the  loop. 
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While  in  all  the  structures  studied  until  now,  the  loop  transmission 
adjusts  itself  automatically,  the  EEAL  adaptive  properties  are  only 
effective  when  the  secondary  loop  achieved  its  steady  state.  Note  the 
difference  between  the  SOAL  and  the  EEAL.  The  SOAL  has  inherent  adaptive 
properties,  even  when  the  secondary  loop  has  not  achieved  its  steady 
state.  In  the  transient  situation  quasi-linearity  constraints  can  be 
violated,  if  oscillations  are  smaller  than  the  steady  state  oscillation 
and  an  extreme  input  or  disturbance  is  applied,  or  plant  output 
oscillation  level  might  exceed  the  specified  level. 

In  the  EEAL,  since  the  input  in  the  non-linearity  is  predominantly 

defined  by  the  injected  signal,  if  there  is  a change  in  the  plant  gain 

and  the  secondary  loop  is  not  able  to  "follow"  this  change,  adaptive 

properties  are  not  kept  during  this  dynamic  process,  and  the  minimum 

gain  margin  of  the  main  loop  should  be  kept,  to  avoid  self-oscillation, 

taking  into  account  the  possible  "transient"  changes  of  L^(ju)  . So, 

AK 

"abrupt"  changes  should  be  specified. 

Note  - There  is  a possibility  to  inject  a signal  ZT(t)  with  a variable 
amplitude  based  on  V (see  Figure  6.2-2) , which  would  partly  avoid 
this  problem,  but  no  studies  or  simulations  were  performed. 

Dynamics  of  the  Adjustment  Loop 

As  in  the  SOAL,  Section  4.5,  it  is  desirable  to  obtain  the  largest 
possible  crossover  frequency  for  the  secondary  loop. 

6.6-2  The  Constraints 

The  EEAL  has  to  satisfy  the  same  constraints  defined  for  the  EEAS, 
Section  3.5-2,  with  the  difference  that  self-oscillation  "quenching"  is 


S .v  turn 
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more  complex. 


Quasi-linearity  Constraints 


max  x. (t)  s — 
t 1 1 1 a 

u 


i = r,d 


. o 

“bi  * T 


i = r ,d 


a = 3 
6 = 3 

These  conditions,  if  satisfied,  permit  us  to  write: 
v = v + v + v. 


X = X + X + X 

o r 


y = y + y + y , 
1 Jo  r *d 


y . = N.x. 

i li 


Mf  2B 

= — ( = — for  ideal  relay) 


M 

N = ~r 

O A 


2M 


—■  ( = ~ for  ideal  relay) 
A A 


(6.6-la) 

(6.6-lb) 

(6.6-lc) 

(6.6-ld) 

(6.6-2a) 
(6 . 6-2b) 
(6.6-2C) 
(6.6-2d) 

(6.6-2e) 

(6 . 6-2f ) 


B is  the  output  of  »Ms)  , defining  the  saturating  level  of  the  ideal 
relay  (Figure  6.2-2). 

Self-oscillation  "Quenching” 

It  has  been  shown  (Section  3.5-2)  and  explained  the  problems 
involved  if  a self-oscillation  exists  in  the  loop.  For  the  EEAL  the 
problem  is  more  serious  since  the  system  is  not  "dynamically  adaptive", 
so  one  should  take  the  loop  transmission  gain  margin  constraint,  to  cope 
with  transient  changes  of  L or  errors  in  the  measurement  processes 
during  the  transient  period,  due  to  the  command  input  or  to  disturbances. 
A gain  margin  of  12  dB  seems  to  be  reasonable  to  cope  with  "errors"  in 
the  measurement  processes.  If  there  are  expected  abrupt  gain  changes 


of  more  than  6 dB 


ations  should  be  considered 


Maximum  Plant  Output  Oscillation 


max 


EEAL  - Smooth  Solution 


The  basic  EEAL  structure  in  repeated  in  Figure  6.6-1 


Figure  6.6-1  - Basic  EEAL  Structure 


-195- 


! 


I 

i 

I 

t 


u 

* 

U '! 

< 


f 


As  in  the  SOAL,  there  exist  two  types  of  loops.  One,  the  main  loop, 

valid  for  the  forced  and  disturbance  signals,  is  directly  related  to 

L (jw  ) and  valid  only  for  these.  The  second,  called  the  secondary 
o o 

loop,  is  responsible  for  the  adjustment  of  the  relay  saturation  level  B , 
assuring  to  VQ(t)  a constant  amplitude  V 


The  equations  below  follow  the  same  development  reasoning  and  have 
the  same  constraints  as  in  the  EEAS,  Equations  3.5-5  - 9. 


Lf(jw)  = NfG1G2KPh( ju) 

(6.6-5a) 

or 

Nc  = — (for  the  ideal  relay) 
f TTA 

(6.6-5b) 

VjV  = NoGlG2ICPh(>) 

(6.6-5c) 

4B 

Nq  = — (for  the  ideal  relay) 

(6.6-5d) 

c L.(jcb) 

Tr(jw>  = r(30»  = G3(:“)  • 1 +Lf(j(0) 

(6.6-6) 

Cr(jw) 

Xr(jUl)  " NfG2KPh(jw) 

(6.6-7) 

X RT  (jw) 

, i A MfG2KPh(juj) 

(6.6-8) 

1 j xr  11 

!!1  T(t)  ‘I-T 

(6.6-9) 

X 

Assuming  a "smooth"  — (jco) 

x 

r , . . „ 1 1 

T(3U)) 

(6.6-10) 

Comparing  Equations  6.6-8  and  10, 

Xr  RTr«“>  ,1 

A M,G_KP,(jw)  ' aw 

f 2 n 

(6.6-11) 

There  is  now  a difference  between  the  EEAS  and  the  EEAL,  since  there 
is  no  "most  difficult"  situation  for  Equation  6.6-11  to  be  satisfied, 
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as  regards  the  plant  gain  and  amplitude  of  oscillation,  since  MfK  and 

A are  constant  values  in  the  steady  state  (Equations  6.2-9,10).  The 

X 

only  variable  affecting  — (joi)  in  Equation  6.2-11  is  the  input  R(jw)  . 
It  is  most  difficult  to  satisfy  this  equation  when  the  command  input  has 
its  extreme  value  rg(t)  > with  Laplace  transform  (s  = jw)  , Rg  ( jtu ) 


Substituting  in  Equation  6.2-11, 

| MfG2KPh  ( ju'  * a • lo  • iRgT^jw)1 


(6.6-12) 


A similar  expression  to  Equation  6.6-12  is  obtained  concerning  the 


extreme  disturbance  input. 

Lf  (jo)) 

Zd(j“)  = "D(jW)  ' 1 + Lf W 
Zd( jw) 

Xd(jU)  = NfG2KPh(ja)) 

BUt  Lf(jw) 

Zd(jU)  = -D  ( jw)  • r;  L-(— r 


(6.6-13) 


(6.6-14) 


Substituting  these  two  equations  in  6.6-14, 


Xd 

T(j^) 


-D(jw) 


Lf ( jw) 


MfG2KPh(jt0)  1 + Lf(3w) 


(6.6-15) 


Considering  u "smooth"  — (ju>) 


d 11 

— (ju)  $ — = T" 
A au)  3w 


(6.6-16) 


Combining  Equations  6.6-15  and  16, 


V D ( ju))  . Lf(jt0)  < J_ 

A MfG2KPh(jw)  1 + Lf  ( jjo)  * aw 


(6.6-17) 


In  this  last  inequality  the  "most  difficult"  condition  to  be  satisfied 
is  when  the  extreme  disturbance  de(t)  is  applied.  So: 


L (ju) 

iMfG2KPh(j.  • I * a • u>  • |De(jcu)  | lVLfj ju)- 


(6.6-18) 
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Conditions  6.6-12  and  18  being  satisfied,  assure  amplitude  quasi-linear- 


ity constraints  are  not  violated  if  r(t)  $rg(t) 


d$d  (t) 
e 


There  is  less  explicitness  in  Equation  6.6-18  if  compared  with  12, 
since  T^fju)  is  "amost"  determined  and  Lf(jio)  is  unknown,  being  the 
main  purpose  of  the  synthesis  procedure. 


Next,  using  Equation  6.6-5, 


IWV^oH  * 2 


(6.6-19) 


Since  MfK  remains  constant  (Equation  6.2-10),  there  is  no  "most 
difficult"  case  to  satisfy.  Recalling  Equations  6.6-11  and  17,  and  re- 
placing w = u)q  , by  extrapolation  (Appendix  III) 


(6.6-20) 


I 1 

lx  1 

1 w 

ReTr(jwo) 

' 

1 i 

1 1 
l 

xd 

T(^o} 

CM 

MfG2KPh(jwo) 

Lf(jo>0> 


1 + Lf(joio) 


Substituting  Equation  6.6-19  in  20  and  21, 

lx 


(6.6-21) 


f(jwo> 


Xd 

T<jwo} 


* - T R (ju>  ) 
m 1 r e o 


a — D ( jw  ) 
m 1 e o 1 


Lf(j“o) 


1 + Lf  ( ju>Q) 


(6.6-22) 


(6.6-23) 


oscillating  frequency  for  a "smooth"  X^_  ( joo ) and  X^(ju) 


Equations  6.6-22  and  23  give  non-linear  functions  which  define  the 

These  two 

results  give  a big  advantage  compared  with  those  obtained  for  the  EEAS 
or  SOAL. 

Comparing  the  first  EEAL  result.  Equation  6.6-22,  with  the 

K2  ' 

equivalent  EEAS  one,  Equation  3.5-30,  the  factor  — 

K1 


has  disappeared 
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in  Equation  6.6-22,  which  is  already  an  advantage  for  the  EEAL,  and 
| L^.  ( j ) | has  also  disappeared  from  the  denominator,  which  could  also 
be  an  advantage  for  the  EEAL  if  |Lf(jwQ)  |«1  . There  is  no  difference 

when  comparing  Equation  6.6-22  with  the  equivalent  SOAL  result, 

Equation  4.5-23. 

Next,  comparing  the  second  EEAL  result,  Equation  6.6-23,  with  the 

k2 

EEAS  one,  Equation  3.5-31,  the  factor  — has  disappeared  and  D ( jui  ) 

6 O 

Lf<j(V 

is  multiplied  by  — — in  Equation  6.6-23,  instead  of  being 

1 + L^ 

multiplied  by  . — in  Equation  3.5-31,  and  this  gives  the 

l + Lf(;)u)o) 

flexibility  the  EEAL  was  looking  for.  In  the  EEAS  the  tendency  was  to 
have  !Lf(juQ) | »i  , to  cope  with  the  disturbance  input.  Since  , 

and  |Lf (jw^) | < 1 , the  EEAS  loop  transmission  asked  for  a relative 

"peaking"  at  it)Q  , with  all  the  problems  involved  in  it,  as  demonstrated 


in  the  EEAS  numerical  example,  Section  3.6. 


Lf(ju>o) 


In  the  EEAL,  D (jw  ) being  multiplied  by  ^ — — leads  to 

e o r 1 + Lf  (3“q) 

] L ( ju  ) I « 1 , so  w > w and  I L,  ( ju  ) I < I L_  ( jto  ) I <1  is  obtained 

't  o o it  ' £ o 1 ' t w 1 

2 

in  a more  "natural"  loop  shaping  manner,  as  shown  in  Figure  6.6-2. 

Comparing  this  same  EEAL  result.  Equation  6.6-23,  with  the  SOAL 


equivalent  one,  Equation  4.5-23,  the  same  relative  advantages  exist  as 


above , due  to 


yj'V 

1 + Lf  ( juQ) 


1 for  the  SOAL,  so  Equation  4.5-24  is  a 


special  case  for  Equation  6.6-23. 


In  all  structures  seen  before,  the  disturbance  has  been  the  main 
problem  to  handle,  especially  when  the  equivalent  applied  disturbance 

signal  at  the  plant  output  was  a step  function.  It  seems  there  is 


limit  on  how  small 


can  be  in  Equation  6.6-23.  The 


"physical”  explanation  for  this  fact  is  due  to  the  breaking  of  the  link- 
age between  cQ(t)  , the  plant  output  oscillation  signal,  and  xQ(t)  , 
the  non-linearity  input  oscillation  signal.  This  has  already  been 
discussed  in  Equation  6.2-12. 


But  there  is  a limitation  on  how  small 
as  per  Equation  6.4-1, 


V(juc) 

Lo(juo) 

Kjwo) 

1 + Lo(jwo) 

can  be  since, 


(6.6-24) 
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« 1 


then 


Lf(j“o) 


« 1 , and 


1 + Lf  (ju)Q) 


Vj“o) 


w !Lf (j“Q) 


Lo(j<V 


1 + L (ju  ) 
o J o 


« 1 


This  last  conclusion  means  in  Equation  6.6-24  that  |V(ju>o)  X ( j ) , 

so  V , the  amplitude  of  oscillation  at  v(t)  , is  also  very  small 
compared  to  the  injected  signal. 

(6.6-25) 


x (t)  = v (t)  + J(t)  « J(t)  = A sin(u)  t) 
o o o o 

The  quasi-linearity  condition  being  satisfied  at  x means 

lx. 


SO 


v(t) 


-(t) 


*3  ' 


could  achieve  large  values,  since  — «1  , as  per 

o 


Equations  6.6-24  and  25,  with  quasi-linear  conditions  probably  not 

satisfied,  making  it  difficult  to  measure  the  amplitude  of  oscillation 

of  v (t)  , and  "obtain"  the  real  information.  This  fact  gives  the 

o 


limit  on  how  small 


Lf(julo) 

1 + Lf  (ju>o) 


may  become. 


The  measurement  of  this  oscillation  is  more  difficult  to  obtain 
than  in  the  SOAL,  where  a "peak  to  peak"  detector  was  used.  The  used 
technique  and  approximations  have  already  been  described  in  Section  6.5 
(in  this  technique  the  information  of  the  phase  /Lp  ( j coo ) is  used)  . 


Finally,  if  Z (jw)  and  Z .(jw)  are  the  extreme  plant  output 
er  ed 


due  to  the  forced  signal  and  disturbance. 


Zer(jw)  = ReTr(jt0) 


De  Lf(jw) 
Zed(ju)  = 1 + Lf( ju) 


(6.6-26) 


(6.6-27) 


In  Equations  6.6-26  and  27  u = <i>o  , and  substituting  in  22  and  23, 


1 
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A final  expression  is  derived.  Combining  Equations  6.6-17  and  19, 

lx. 


3cu 


■T(jw o’!  4 m lDe(jwo) 


Lf(ju)o) 


1 + Lf ( juo) 


(6.6-31) 


Considering  the  two  extreme  factors  in  this  inequality: 

1 


| Lf(jwo) 


1 + ( j(jjQ) 


6oj  T D (j<i)  ) 
o e J o 


(6.6-32) 


Then  the  Nichol's  chart  is  used  as  in  the  EEAS.  Disturbance 
attenuation  constraints  and  self-oscillation  "quenching"  are  translated 
as  bounds  for  the  different  frequencies,  as  per  Section  3.5-3, 

Figure  3.5-4. 

The  bound  due  to  constraint  6.6-32  is  also  plotted  in  the  same 

Nichol's  chart  for  several  values  of  u>  , which  is  still  unknown.  The 

o 

technique  used  to  define  Lf(jii>)  is  very  similar  to  the  one  used  in  the 
EEAS,  Section  3.5-3. 

In  the  numerical  example.  Section  6.7,  the  synthesis  procedure  is 
followed  step  by  step. 


6.7  Numerical  Example  - EEAL  - Smooth  Solution 

The  chosen  EEAL  structure  is  repeated  in  Figure  6.7-1. 

Th<->  defined  set  of  specifications  and  constraints  in  this  numerical 
example  is: 

Desired  Command  Input  Response 

1 


Tr(s)  = 


, , 2*.5  . s 

1 + ~T~  5 + ^ 


X $>  (s) 


—7 




Figure  6.7-1  - EEAL  Structure 


where  4>  (s)  includes  far-off  poles  and  zeros  that  can  be  added  to  the 
"nominal"  desired  transmission. 

Minimum  Disturbance  Attenuation 

The  minimum  disturbance  attenuation  required  is  plotted  in 
Figure  6.6-2.  Note  there  are  no  specifications  for  u>  5 u>„  = .1  rd/sec 

H 

The  Plant 

P(s)  = | 1 $ Kx  5 K $ K2  = 100 

"Abrupt"  changes  of  6 dB  in  the  plant  gain  are  possible. 

Extreme  Command  Input 

r (t)  = 75  u(t)  -•>  R (s)  = 25. 
e e s 

Extreme  Disturbance  Input 

d (t)  = 20  u (t)  -*>  D (s)  = ^ 

e e s 

The  disturbance  response  should  not  be  oscillatory. 
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Quasi-linear  Constraints 

The  used  parameters  are: 

a = 3 

6 = 3 
so 


X 

i 

V , 

1 

0) 

o 

U) 

o 

T(t) 

<3  ' 

A (t) 

* 3 ' 

V‘T  ' 

wbd  * T 

Plant  Output  Oscillation 
|CQ(t) | * m - 1 

Limit  Cycle  Quenching 

In  Equation  6.6-3  we  use  p ' = , taking  in  consideration  an  extra 

6 dB  (if  compared  to  the  EEAS  problem)  to  cope  with  the  specified 
possibility  of  6 dB  "abrupt"  changes  in  the  plant  gain.  So, 

IVj(Vl  *2  ■■  lLf(ju)1r)!  * 4 = (‘12  d®) 

Excess  of  Poles  over  Zeros 

Lf(jui)  is  to  have  an  excess  of  5 poles  over  zeros. 

The  Solution 

In  the  Nichol's  chart  (Figure  6.7-3)  all  the  existing  bounds  on 
L^(jo))  are  plotted,  for  several  frequencies.  The  bounds  are  due  to: 

1 - Disturbance  attenuation  specifications,  as  per  Figure  6.7-2,  and 

valid  for  bjb  = .1  rd/sec  . w should  be  larger  than 
H O 

ui()  ; (,)('. (|  - . 3 rd/sec  . The  technique  used  is  the  same  cne  used  in 
Sections  2.6  (SOAS) , 3.6  (EEAS),  etc.  The  bounds  on  Lf(jw)  are 
plotted  for  a certain  number  of  frequencies. 


—v 


VK>>  | 

. . ■ ',■■■  ' \ = -42d 

t t l.(  ( v*'(  i 

1'  Itfilf  n 

i •.  r i 

1 1 . A 1 . 

Niimr i 1 • .i  1 Mx.implo  - Hounds  on 
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w (rd/sec) 

A (dB) 

.007 

-24 

.013 

-18 

.05 

-6 

.075 

-3 

.1 

0 

2 - Equation  6.6-32,  repeated  here,  provides  a bound 

yjy  < m 

l + Lf(jwo)  * 6u>o|De(ja>o)  | 

20 

In  this  numerical  example,  m = l , D (jto  ) = — , so 

e o co 

o 

Lf(jw  ) ^ 

r~- — : r $ ttTT  (=-41.56  dB) 

1 + L,  (ini  ) 120 

f o 

Note  that  in  this  numerical  example  the  resulting  bound  is  not  a function 
of  the  "possible"  oscillation  frequency,  but  a constant  value.  This 
bound  is  plotted  in  the  Nichol's  chart,  at  -42  dB. 

3 - Self-oscillation  "quenching"  - Gives  a gain  margin  condition  of 

12  dB,  so  L^(jw  ) <~  ( = -12  dB)  . This  bound  is  marked  as  the  point  Q 

f TT  4 

in  Figure  6.7-3. 


4  - Non-oscillatory  disturbance  response  requirement  - A good  engineering 

practice  which  generally  satisfies  this  requirement  is 
Lf  (jw) 

— — r < 2.3  dB  , V a)  . This  closed  loop  condition  is  also 

l+Lf(}(o) 

plotted  as  a bound  in  the  Nichol's  chart,  Figure  6.7-3. 


5  - L.^(jm)  has  to  satisfy  all  these  bounds  and  is  the  minimum 

frequency  allowing  Lf(jo))  to  satisfy  these  bounds.  It  has  been 


■li  ,i  iiiiwmi** 
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shown  1 that  the  optimum  loop  transmission,  in  the  sense  defined  in 
that  reference  ([5]),  is  the  one  that  "stays"  in  its  respective  bound 
for  each  frequency. 


The  chosen  is 


Lf(s)  = 


which  has  an  excess  of  5 poles  over  zeros  as  specified. 


The  oscillation  frequency  u)Q  is  the  one  corresponding  to  the 

intersection  of  L-(jw)  with  the  bound  — r =-42  dB  . So 

f l + Lf(3wQ) 

w = 1.375  rd/sec 
o 

Lf(jco)  is  plotted  in  Figure  6.7-3. 


Once  is  known,  it  is  possible  to  define  the  various  parts  of 

the  loop.  The  limits  of  iM^G^KP^t jw) | , due  to  Equations  6.6-18,12  , 
are  repeated  here, 

|MfG2KPh(ju>)  I *«.-»•  I Vju,)  I 1 (ji) 

as  curve  number  one  in  Figure  6.7-4. 

|MfG2KPh(jw)  | * a • a)  • |ReTr(jw)| 

as  curve  number  two  in  Figure  6.7-4.  In  the  first  expression,  a = 3 , 


|De(jco)J=—  , and  Lf(ju>)  = 


, and  this 


limit  is  completely  defined. 


.8  ' | s = ju> 


In  the  second  expression,  a=3  , 


T (ju>)  * 

r . A 2 x .5  a s^ 

1 + “71“  5 + ~2  . 

.1  s = ju> 


R ( jw)  = — and 
1 e 1 a) 


x4>(jw)  . 4>(jw)  represents  the  far  off 
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poles  and  zeros  which  have  not  been  yet  defined,  and  leaves  some  freedom 

to  the  designer  to  decide  what  can  be  added  to  the  nominal  T^tjw)  , 

"without  affecting"  the  system  response.  In  our  case,  we  choose 
Lf  ( jto) 

<Kjui)  = , -■  — — — — (this  makes  the  G„  calculation  and  implementation 
1 + Lf  (]w)  3 

very  easy) . 


The  transmission  T is 
r 


Tr(jw)  = 


1 + 2 x ‘5  s + — 
.1  .1 


Lf ( jm) 

1 +Lf  (ju) 


|s  = 3U 

As  noted  earlier,  it  is  up  to  the  designer  to  choose  the  far-off  poles 
and  zeros  to  be  added  to  the  transmission.  4>(s)  was  chosen  in  this 
case  for  reasons  of  easy  implementation  of  G. 


Since 


1 + Lf ( jmo) 


satisfies  Equation  6.6-32,  and  Im^G^KP^ jw) 


satisfies  6.6-18,  Equation  6.6-19  repeated  here  is  automatically  satisfied. 
lMfG2KPh(j“o)l  S I = ^ ( = ‘6  dB) 


One  can  see  in  Figure  6.7-4  that  at  u>  = oj  =1.375  rd/sec  , if 
| ( jmQ)  | in  curve  © is  equaled,  it  has  the  above  value  = -j=-6 


dB. 


In  this  numerical  example  the  limit  of  | M^G2KP^ ( jw) | due  to  the  extreme 

command  input,  curve  , at  wo  , is  lower  than  that  due  to  extreme 

disturbance  input,  curve  at  wo  , so  the  last  limit  prevails  and 

there  is  no  problem  with  the  maximum  plant  oscillation  constraint.  In 

case  this  is  not  so,  more  far-off  poles  could  be  added  to  , until 

|MrG  K.P  ( jw  )|.  . _ i |M,G_KP,  (ju>  )|.  _ ^ 

f 2 ho  due  to  R f 2 ho  due  to  D 

i 1 e 1 1 e 

In  case  more  far-off  poles  significantly  affect  the  desired  response, 
then  u)0  has  to  be  increased  until  the  above  inequality  can  be  satisfied. 


|MfG2P(jw>[ 


Figure  6.7-4  - EEAL  Numerical  Example 
Bounds  on  |m  g P(jio)|  and  Solution 
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We  choose 


28  l1+fJ 

M G,KP  (S)  = ' <- - 

Jy  + 2x-6  + s2  \2 

8 {1+~S  + ^2j 


which  satisfies  the  limits,  is  plotted  in  Figure  6.7-4  as  curve  3 


Now  we  can  calculate 


G1  Lf(s) 

T(S)  ' MfG2KPh(s) 


“(“¥•* 4)2  * f )‘ 


Since  M^K  remains  constant  (Equation  6.2-10) , when  the  plant  has  its 

minimum  gain  K = K = 1 , M,  will  have  its  maximum  value  M = M * M, 

I i i ii  i max 

2b 

But  M,  = — , and  if  we  choose  B,  = B = 30tt  M , , = 60 . When  K = K_  = 100 , 

f it  1 max  f 1 2 

K1  K1 

then  Mf2  = Mfmin=MflXi^=-6  ' and  B2  = Bmin  = BmaxX^=-3,T  * Since 
MfG2KPh(s)  =MflG2KlPh(s)  =Mf  2G2K2Ph(s)  ' 

MfG2KPh(s) 

G2  lS)  M.  G_K,  P,  (s) 
t l z l n 


V., 


Choosing  arbitrarily  A = 140  , 


. . 2 x .6 
1 + - — s + 

^ 


G (s)  \ .5 

G.  (s)  = — xA  = r ~ 

1 A / „ _ 2 \ 2 

( , . 2 x .7  . s 1 

^ ~ ' j) 


iL_ 

)Vt): 


The  amplitude  of  the  injected  signal  Aq  is  defined  by  (see  Equation  6. 3-1) 
A I 1 I I V^O*  I I 1 I 
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k 

i '! 

i 


In  the  Nichol's  chart.  Figure  6.7-3,  Lo(jii)o)  =2Lf(jwQ)  , and  the 
following  values  are  read: 


L ( ju>  ) 
o o 


so 


1 + L ( jo)  ) 
o J o 


A 

r W1 

o 


« -36  dB 


A w a = 140 

o 


L (jw  ) » -36  dB 
o o 


Since  o)q  = 1.375  rd/sec  , the  injected  signal  is 
y = 140  sin(l. 375  t) 


The  amplitude  V of  VQ(t)  can  be  calculated. 


_V_ 

A 


1 + L ( jw  ) 
o J o 


= -36  dB 


("&) 


v-W— » 

In  the  extreme  case,  when  quasi-linearity  conditions  are  most  likely  to 
be  violated. 


-(t) 


|xf (t) | « j = 46 


The  forced  signal  ( v^  or  v^  ) can  be  20  times  larger  than  the 
amplitude  of  oscillation  V of  v 


Due  to  the  far-off  poles  and  zeros  definition  in  T^fjw)  , 


G3(s)  = 


, . 2 x .5  s‘ 

1 + — I”  3 + ^2 


Finally  the  function  iJj(s)  is  defined  based  on  the  maximum  possible 

secondary  loop  bandwidth,  when  the  plant  has  its  maximum  gain 

(liquation  6 . 5-9)  . 

. , . .02 

it>(s)  = 
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6 . 8 EEAL  Numerical  Example  - Simulation  Results 

A digital  simulation  of  the  complete  system  was  performed  on  the 
CDC  system.  First  the  system  "response"  is  shown  when  the  secondary  loop 
is  open  ( ip  (s)  = 0 ).  The  results  show  mainly  the  oscillation  amplitude 
measurement  dynamic  properties  when  there  is  an  abrupt  change  in  the 
gain,  and  the  oscillation  amplitude  measurement  error  when  an  extreme 
input  or  disturbance  is  applied.  As  a reference,  the  system  response 
for  these  signals  is  also  shown  (Figures  6.8-1  - 4). 

One  can  see  that  for  abrupt  changes  in  the  gain  of  the  plant,  after 
one  cycle,  the  measurement  has  "almost"  been  completed.  When  the  extreme 
command  and  disturbance  are  applied,  the  change  in  the  measurement  is  an 
error,  since  the  plant  keeps  its  gain  value,  and  there  is  no  change  in 
the  oscillation  amplitude,  and  this  error  is  equivalent  to  a "disturbance" 
in  the  EEAL  secondary  loop. 

Next,  the  secondary  loop  is  closed,  and  (Figures  6.8-5  - 7)  the 
extreme  command  input  and  (Figures  6.8-6  - 8)  the  extreme  disturbance 
are  applied  for  several  values  of  the  gain  of  the  plant.  The  results 
are  close  to  the  previous  (secondary  loop  open)  ones  obtained. 

Finally,  the  change  in  the  dynamics  of  the  secondary  loop  is  depicted 
in  Figures  6.8-9  - 14.  The  results  shows  First,  the  lower  the  gain  of 
the  plant,  the  lower  the  dynamics  of  the  adjustment,  as  per  Equation  6.5-9; 
Second,  the  system  response  is  not  the  same  when  the  plant  gain  changes 
upwards  or  downwards.  This  is  due  to  the  reasons  explained  at  the  end 
of  Section  6.5.  The  results  show  the  capability  of  the  system  to  cope 
with  relatively  large  step  disturbances. 


. ..... 


Oscillation  Amplitude  Measurement  Oscillation  Amplitude 


Oscillation  Amplitude  Measurement 


Figure  6.%- 3b  - EEAL  System  Response 

Secondary  Loop  Open  { *<s>  -0  ) Figure  6.0-3a  - EEAL  System  Response 


Y50  - r - Applied  Command)  FZ  = z - Plant  Output 

Y50  =»  r - Applied  Coi 


Figure  6.8-3d  - EEAL  System  Response  Figure  6.8-3c  - EEAL  System  Response 


o 


o 

n 


FXG1  « v - Output  of  Gj  FZ  ” z - Plant  Output 

Y50  » r - Applied  Command  Y50  - r - Applied  Ca 


z - Plant  Output  FX » x - Non-linearity  Input 

Plant  Gain  K - Plant  Gain 


Figure  6.8-13b  - EEAL  System  Response  Figure  6.8-1’a  - eeal  System  Response 

Complete  System  Simulation  Complete  System  Simulation 


Measured  Amplitude  A * Plant  Gain 


Complete  System  Simulation  Complete  System  Simulation 

rnamics  of  the  Adjustment  Loop  Dynamics  of  the  Adjustment  Loop 
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CHAPTER  SEVEN 

E.E.A.N.L. 

(EEAL  with  a Non-linearity  in  the  Secondary  Loop) 

7 . 1 General 

Chapter  6 dealt  with  the  EEAL  structure,  which  presented  significant 
advantages  over  all  the  other  structures.  The  EEAL  was  able  to  cope  with 
almost  any  problem,  with  respect  to  the  extreme  command  and  (mainly)  dis- 
turbance response.  The  limitation  is  due  to  the  capability  of  measuring 
the  amplitude  V of  an  oscillating  signal  VQ(t)  , combined  with  an 
aperiodic  signal  vf(t)  , which  can  be  much  larger  than  the  oscillation 
amplitude.  The  main  disadvantage  of  the  EEAL  is  the  "reappearance"  of 
the  factor  K , gain  of  the  plant,  in  the  secondary  loop,  which  deter- 
mines the  rate  of  adaptation  of  the  system.  The  later  is  given  by 
Equation  6.5-9,  repeated  here. 

f(s)  » zillLeJl (7.1-1) 

AV  -Ts 

1 + 7 G1G2V'(S)  K e 

The  idea  of  the  EEANL  is  to  keep  the  advantages  of  the  EEAL,  and 
"eliminate"  the  factor  K from  the  secondary  loop,  in  Equation  7.1-1 
denominator.  The  increment  EEAL  -*■  EEANL  , is  basically  the  same  as 
the  SOAL  -*  SOANL  one. 


The  "elimination"  of  K in  Equation  7.1-1  could  be  achieved  by 
introducing  a multiplier  in  the  secondary  loop,  with  multiplication  factor 
inversely  proportional  to  K . But  K is  not  available  directly,  so 
B , which  is  related  to  K , can  be  used,  since  in  the  steady  state 
KB  is  constant,  as  per  Equation  6.2-10. 


L 


Figure  7.1-1  - EEANL  Structure 


The  multiplier  changes  the  dynamics  of  the  secondary  loop,  as  analyzed 
in  Section  7.2. 
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Most  of  the  transmissions  are  the  same  as  in  the  EEAL,  Section  6.5, 


and  repeated  here. 


fcoi(*>  = 0 

-T  S 

t (s)  « e ° 

Cl 

tcs(s>  “ 7 GlG2™h,Sl  e'T°S 
'os'8'  ■ 1 


(7.2-3) 


(7.2-4) 


(7.2-4) 


(7.2-4) 


The  only  function  to  be  defined  is  ifi'  (s)  . We  can  write  (Figure  7.1-1) 


B(t)  = y ( t)  * 


[(vr~  V * B(t)] 


(7.2-5) 


where  y 


= X>_1w 


Assuming  B(t)  does  not  change  "very  much"  from 


its  initial  value  B.  , 

l 

B ( t)  = B.  + AB. (t) 

i l 


B ( t)  & p(t)  * 


[(vr-vm)  • B.(t)] 


(7.2-6) 


B(s)  w iMs)  [VR(s)  - Vm(s)]  • (b^~) 

' min  ' 


ip'  (s)  « ip  (s) 


(4) 


(7.2-7) 


(7.2-8) 


Substituting  the  above  in  Equation  7.2-1, 


-(b-£-)<’(s)  e T°S 

' min' 


AV  KB.  _T  s 

1 + 7giW(s)  -r6  ° 


(7.2-9) 


In  Equation  7.2-9's  denominator,  B^  is  the  initial  value  of  B(t)  , 
and  assuming  the  steady  state  was  achieved,  B^K = constant 
(Equation  6.2-10) , and  this  expression  is  independent  of  K . The 
dynamic  of  the  secondary  loop  is  independent  of  the  plant  gain.  The 
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assumption  that  B (t)  does  not  change  "very  much",  once  the  secondary 
loop  steady  state  is  achieved,  is  valid,  considering  the  dynamics  of  the 
adjustment  loop  being  faster  than  the  changes  in  K . Simulations  have 
shown  (Section  7.3)  than  even  considering  bigger  changes  in  K and 
consequently  in  B , the  results  are  still  valid  for  Equation  7.2-9. 

When  K = K , then  B = B . , and  Equation  7.2-9  is  identical 

max  min  ^ 

to  the  EEAL  equivalent  equation  (6.5-9),  if  all  the  other  functions  are 

the  same.  This  means  that  when  K = Kmax  > t*ie  EEANL  secondary  loop 

dynamics  is  identical  to  the  EEAL  one.  Suppose  now  K = K . , then 

min 

K 

_ t,  maX 

B = = x v » and  the  denominator  of  Equation  7.2-9  has  the  same 

max  mm  K . ^ 

min 

function  as  before  (when  K = K )»  and  the  dynamics  of  the  adjustment 
loop  is  kept  constant. 

EEANL  Sensitivity  to  Dynamic  Changes 

Section  6.4,  EEAL  is  valid  for  the  EEANL. 

7 . 3 Synthesis  Procedure  and  Numerical  Example  - EEANL 

The  synthesis  procedure,  the  specifications,  data  and  constraints 
are  the  same  as  for  the  EEAL,  Section  6.6. 

The  EEAL  main  loop  synthesis  procedure  and  results  are  also  valid 
for  the  EEANL.  The  advantage  is  in  the  secondary  loop,  where  the  EEANL 
synthesis  procedure  is  the  same  as  the  EEAL  one,  but  the  results  provide 
the  dynamics  of  adaptation  be  almost  constant. 

The  same  specifications  of  the  EEAL  numerical  example  are  taken  in 
the  EEANL  numerical  example.  As  a result,  all  the  EEANL  blocks  have  the 
same  functions  as  the  EEAL  blocks,  including  i(/(s)  . The  only  difference 

is  the  added  multiplier  in  the  secondary  loop. 
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Some  implementation  details  mentioned  in  Section  5.3  for  the  SOANL 
are  also  valid  for  the  EEANL. 

EEANL  simulations  were  performed  on  the  CDC  6400  system.  Results 
are  presented  in  the  next  figures,  and  show  that  the  EEANL  structure, 
which  is  the  most  complex  of  the  studied  structures,  can  cope  with  very 
large  inputs  (command  or  disturbance)  and  severe  plant  gain  changes. 

Figures  7.3-1  to  7.3-6  show  the  system  response  when  abrupt  changes 
of  gain  (±6  dB)  are  applied.  The  dynamics  of  adjustment  is  the  same  when 
an  equivalent  gain  change  is  applied  independently  of  the  initial  plant 
gain  value. 

Figures  7.3-7  to  7.3-9  show  the  system  response  when  the  extreme 
command  input  is  applied. 

Figures  7.3-10  to  7.3-12  show  the  system  response  when  the  extreme 
disturbance  input  is  applied.  Notice  that  quasi-linearity  constraints 
are  satisfied  at  x and  not  at  v , and  that  some  errors  are  introduced 
in  the  measurement  processus  of  the  amplitude  V (obtained  from  v ) 
during  the  transient  state  due  to  the  applied  command  or  disturbance 


input. 


Plant  Gain  I'  - Plant  Gain 

Output  of  G^  Vm  - Measured  Amplitude  of  Oscillation 

B - Non-linearity  Saturation  Level 


Measured  Amplitude  of  Oscillation 
lon-linearity  Saturation  Level 


-linearity  Saturation 


K = 100 

Applied  Command  - Applied  Command 

Ion-linearity  Input  z - Plant  Output 


Applied  Command  rg  _ Applied  Command 

Ion-linearity  Input  z - Plant  Output 


i-linearity  Saturation 


filiation  Amplitude  Measurement  v - Output  of 

-linearity  Saturation  Level 


CHAPTER  EIGHT 


Conclusions 

This  research  has  provided  quantitative  synthesis  theories  for 
several  oscillating  system  structures,  the  most  innovative  of  which 
are  the  SOAL,  EEAL,  SOANL  and  EEANL . Some  recommendations  for  further 
research  are  herewith  listed. 

1.  Generalization  of  the  synthesis  techniques  for  general  plant 
uncertainty.  This  work  has  been  restricted,  in  the  synthesis  details, 
to  high  frequency  gain  uncertainty  - it  being  the  most  significant 
factor  contributing  to  large  loop  bandwidth  in  the  linear  time  in- 
variant feedback  solution  to  the  same  problem.  Such  a generalization 
was  performed  for  the  SOAS  and  EEAS  in  Refs.  9 , 10  and  the  same 
approach  can  be  used  for  the  SOAL,  EEAL,  etc.  However,  there  will  be 
more  complications  in  the  latter  because  of  their  smaller  values  of 

oj  , such  that  the  significant  plant  bandwidth  is  now  closer  to  the 
o 

value  of  0) 

o 

2.  Use  of  combinations  of  linear  and  nonlinear  elements  to  obtain 

more  complex  zero  sensitivity  properties.  For  example,  in  Fig.  8.1, 

with  P = — with  uncertainty  in  p , it  appears  that  zero  sensi- 
s+p 

tivity  to  p may  be  achievable. 


Figure  8.1  - SOAS  Generalization  - Example 
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3.  More  general  relationship  between  and  Lq  . If  N 

is  a static  nonlinear  element,  then  Lf  = Lq/2  . It  is  possible  to 
achieve  more  complex  relations,  by  using  dynamic  nonlinear  elements. 
This  might  permit,  for  example,  the  bandwidth  of  Lf  to  be  con- 
siderably smaller  than  that  of  and  thus  enable  one  to  tailor  Lf 

to  the  actual  feedback  needs.  Ideally,  one  could  then  tailor  the 
nonlinearity  to  the  specific  design  problem.  The  cascading  of  non- 
linear elements  separated  by  dynamic  linear  elements,  is  a possible 
approach . 

( 

4.  In  the  present  structures,  the  nonlinearity  appears  both  in 
the  feedback  loop  and  in  the  direct  transmission  from  input  to  output 

j (when  the  return  path  in  the  feedback  loop  is  open) . One  could  modify 

the  latter,  as  in  Fig.  8.2,  to  have  a parallel  path  bypassing  the 

i 

nonlinear  element. 
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Knowledge  of  K in  P = K is  available  from  the  value  of  B , the 
saturation  level  of  the  nonlinearity,  which  can  be  used  to  modify  . 
The  latter  could  be  used  to  speed  up,  if  necessary,  the  system  transfer 
function,  without  increasing  w . This  might  be  appropriate  in  the 
SOAL,  SOANL , EEAL,  EEANL  in  those  cases  where  the  ratio  of  desired 
system  response  bandwidth  to  , becomes  a significant  design 

factor. 

5.  Finally,  it  would  obviously  be  very  desirable  to  have  a 
mathematical  theory  for  the  dynamics  of  limit  cycle  build-up,  and  of 
limit  cycle  quenching  by  'large'  amplitude  external  signals. 

6.  It  is  interesting  how  the  factor  K /K  . reappears  in 

max  min 

various  guises  in  the  various  dithered  structures.  It  is  ba  lished 
completely  as  an  adaptive  problem  in  all  of  them.  But  in  the  SOAS  and 
EEAS  it  appears  as  a very  important  factor  in  the  quasilinearity 
problem.  In  the  other  structures  (SOAL,  etc.),  it  is  banished  as  a 
quasilinearity  problem  but  reappears  as  an  important  factor  in  the 
speed  of  adaptation  to  parameter  variations.  One  may  view  nonlinear 
feedback  synthesis  as  a means  of  obtaining  different  packages  of 
benefits  and  costs.  The  linear  time  invariant  feedback  class  presents 
one  package.  The  design  is  theoretically  independent  of  the  nature  of 
the  signals,  as  the  system  presents  the  same  'face'  to  all  signals. 

The  price  paid  for  the  benefits  of  feedback  is  well  defined  and  there 
is  overdesign  in  some  aspects.  One  can  decrease  the  price  by  turning 
to  nonlinear  feedback,  but  then  for  only  certain  classes  of  inputs. 

One  gets  a different  'package  deal'.  Ideally,  a nonlinear  feedback 
synthesis  theory  would  provide  one  with  the  means  of  choosing  the 


• 
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specific  nonlinear  structure  needed  to  obtain  the  ’package  deal' 
suitable  for  his  specific  quantitative  design  problem. 


II 


Figure  1-1  presents  a canonic  single-loop  structure  wherein  there 
is  ignorance  of  parameters  only.  Hence,  due  to  such  ignorance, 

A ln|T(ju)  | = A In  777  , T = 7777  , L = GP 

A In  L = A In  P (la-d) 

A specific  value  of  frequency  is  chosen;  say  07  rps.  The  values 
of  over  the  range  of  plant  parameters  are  calculated  and  the 

bounds  obtained.  The  procedure  is  illustrated  for  the  case 


P(S)  = ■,  . ; 1 S k S 10  ; ISaS  10 

s (s  + a) 

This  is  conveniently  done  in  the  plane  of  lnL(jto)  =ln|L|  + j Arg  L , 
the  abcissa  in  degrees  and  the  ordinate  in  decibels,  the  Nichols  chart. 
Thus,  at  (jl)  = 2 rps  , P(2j)  lies  within  the  boundaries  given  by  ABCD 
in  Figure  1-2.  Since  In  L = In  G + In  P , the  pattern  outlined  by  ABCD 
may  be  translated,  but  not  rotated,  on  the  Nichols  chart,  the  amount  of 
translation  being  given  by  the  value  of  In  G(2j)  . For  example, 

if  a trial  design  of  L(2j)  corresponds  to  the  template  of  P ( 2 j ) at 
A'B'C’D'  in  Figure  1-2,  then 

I G ( 2 j ) | db  = 11,(23)1^  - | P ( 2 j ) | db 

= (-2.0)  - (-13.0)  = 11.0  DB 
Arg  G (2 j ) = Arg  L(2j)  - Arg  P(2j) 

= (-60°)  - (-153.4°)  = 93.4° 

Hound:;  on  l.(jm)  in  the  Nichols  chart 

The  templates  of  l’^(ju>)  are  manipulated  to  find  the  position  of 
L(jrii)  which  satisfies  the  specifications  on  In  | T ( jo>)  | . It  is 
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recalled  that  these  are  in  the  form  T^(w)  £ |T(jw)  | . Taking 

the  ui=2  template,  one  tries,  for  example,  positioning  it,  as  shown  in 
Figure  1-2  at  A'B'C'D'.  Contours  of  constant  ln|L/(l+L)  | are  avail- 
able on  the  Nichols  chart.  Using  these  contours,  it  is  seen  that  the 
maximum  change  in  ln|T|  is  closely  (-0.49)  - (-5.7)  =5.2  DB  , the 
maximum  being  at  the  point  C',  the  minimum  at  the  point  A'.  Suppose  that 
the  specifications  tolerate  a change  of  6.5  DB  at  w=2  , so  the  above 

trial  position  of  | L ( 2 j ) | is  in  this  case  more  than  satisfactory.  The 
template  is  lowered  on  the  Nichols  chart  to  A"B”C”D",  where  the  extreme 
values  of  ln|L/(l  + L)|  are  at  C"  (-0.7  DB)  , A"  (-7.2  DB)  . Thus,  if 


Arg.  Lfl(2j)  =-60°  , then  -4.2  DB  is  the  smallest  magnitude  of  Lft(2j) 

which  satisfies  the  6.5  DB  specification  for  Aln|T|  . Any  larger 
magnitude  is  satisfactory  but  represents  over-design  at  that  frequency. 


The  manipulation  of  the  w = 2 template  is  repeated  along  a new  vertical 
line,  and  a corresponding  new  minimum  of  | L ( 2 j ) | found.  Sufficient 
points  arc  obtained  in  this  manner  to  permit  drawing  a continuous  curve 
of  the  bound  on  Lft(2j)  , as  shown  in  Figure  1-2.  The  above  is  repeated 
at  other  frequencies,  resulting  in  a family  of  boundaries  on  LA(2j) 


Disturbances  and  the  Single  High-frequency  Boundary 


Sooner  or  later  there  must  exist  a frequence  range  in  which 
sensitivity  increase  rather  than  decrease  is  achieved,  because  as  shown 
by  Bode  [2,26],  in  any  practical  feedback  system 

/ In  | S | dw  = 0 , where  S = 

is  the  sensitivity  of  T(joi)  to  uncertainty  in  P . Hence,  if  |s|  <1 
over  some  range,  the  primary  objective  in  a feedback  system,  it  must  be 
balanced  by  another  in  which  I S j > 1 . It  is  easy  to  live  with  this 


MMIMO 
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constraint  in  the  two-degree-of-f reedom  system,  because  at  sufficiently 

large  u>  , |t|  is  negligibly  small,  so  that  large  relative  changes  in 

| T | are  inconsequential.  It  should  be  noted,  however,  that  in  any 

practical  one-degree-of-freedom  system,  this  range  where  |s|  >1  , must 

lie  within  the  important  frequency  range  of  T(ju>)  . In  this  higher 

range  or  sooner,  the  plant  template  approaches  a vertical  line,  because 

at  large  a)  any  realistic  minimum-phase  P(jui)  -+K(ja>)  ^ , p being 

its  excess  of  poles  over  zeros.  Thus  in  the  previous  example 
-2 

P(jw)  ->-ka(jai)  and  the  plant  template  approaches  a vertical  line  of 

length  |Aln(ka) I = 40  db 
1 1 max 

As  noted,  in  the  above  relatively  high-frequency  range  |s|»l  is 
tolerable,  as  far  as  its  effect  on  T(jw)  is  concerned,  because  the 
prefilter  F in  Figure  1-1  attenuates  the  resulting  high  peaking  in 
L/(l+L)  . However,  the  disturbance  response  (Figure  1-1) 

C/D^  = (1+L)  ■*"  = S , is  then  also  very  large,  which  is  generally  not 
tolerable  because  there  is  no  equivalent  filter  available.  Even  if  the 
parameter  ignorance  problem  is  assumed  to  dominate  the  design,  it  is 
necessary  to  consider  the  disturbance  response,  at  least  to  the  extent 
of  adding  the  constraint 


This  constraint  leads  to  the  boundary  of  Figure  1-3 

effective  for  all  w , with  0,  = 180  - 0 a function  of  v , e.g. 

h m 1 

0 = 50°  if  y = 2.3  db  . At  low  frequencies  the  parameter  ignorance 
factor  dominates,  e.g.  B2,  B^  in  Figure  1-3,  so  these  boundaries  contain 
no  part  of  . There  is  an  intermediate  frequency  region  in  which 

part  of  the  boundary  contains  a portion  of  B^  , e.g.  B^  in  Figure  1-3. 


At  sufficiently  large  w there  exists  a frequency  to^  such  that 
becomes  the  complete  boundary  for  all  to  i to^ 


-«90* 


Practical  Suooptimum  L 


Figure  1-3  - Bounds  on  L and  Optimum  L on  Nichols  Chart 


The  Optimum  L(jio) 

It  has  been  shown  [1,5]  that  a realistic  definition  of  optimum  in 

single-loop  design  is  the  minimization  of  K , defined  by 
— 0 

lim  L(s)  = Ks  , where  e is  the  excess  of  poles  over  zeros  assigned 

s-**> 

to  L ( s ) 

It  has  been  proven  [5]  that  the  optimum  L lies  on  its  boundary 
at  each  and  that  such  an  optimum  exists  and  is  unique.  Most 

important  for  the  present  purpose,  is  that  in  significant  plant  ignorance 
problems  the  ideal  optimum  L has  the  properties  shown  in  Figure  1-3, 
i.e.  over  a significant  range  it  follows  along  UV  up  to  the  point  J 

at  which  it  abruptly  jumps  to  infinity  along  WW'W"  and  returns  on  the 
vertical  line  YZ,  whose  phase  is  -90  e deg.  Such  an  ideal  L(jw)  is, 
of  course,  impractical.  A practical  suboptimum  L is  shown  in 
Figure  1-3. 
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APPENDIX.  II. 


Derivation  of  Frequency  Response  Bound  from  Time  Domain  Bound 

Xf (jw) 


The  problem  of  finding  the  precise  bounds  on 


due  to  that 


on  max 
t 


xf  (t) 


, is  a very  difficult  one  in  general,  because  of  the 


integral  relation  between  xf(t)  and  its  frequency  response  Xf(jw) 
It  is  done  here  approximately,  and  the  results  confirmed  by  numerous 
experiments  and  actual  system  design  based  on  use  of  the  approximate 
relations . 


The  basic  point  is  that  the  search  need  not  be  over  the  set  over 
all  possible  transforms,  but  only  over  the  much  smaller  subset  of  those 
which  are  realistic  for  the  problem  under  consideration.  It  has  been 
shown  [9]  that  realistically,  the  search  need  be  only  over  the  subset  of 
"dominantly  first-order"  transfer  functions,  i.e.  those  which  over  the 
significant  frequency  range  are  effectively  first-order  (excess  of  only 
one  pole  over  zeros).  This  means  that  Re  Xf(jw)  >0  over  the  frequency 
range  of  importance.  In  that  study  (9] , the  simplest  first-order  system 


Xf 

was  ther.  used  — = — 2 — giving  max 

A s “b  t 


xf  (t) 


In  order  that  the 


latter  be  £ — , it  follows  that 

a 


Xf 

T(j*) 


= — . This  was  justified  by 

aw 


demonstrating  that  attempted  more  complex  forms  for  Xf(jw)  resulted  in 
basically  the  same  result.  Other  good  reasons  were  given  there,  why 
such  a simple  form  for  X^.(s)  was  desirable. 


Here,  we  accept  the  constraint  Re  Xj.(jw)  >0  and  use  a result  by 

I’afxjulis  I 27]  that  if  R(w)  = Re  F(jw)  i 0 , R ' (o' ) (0  V w , then 

max|s(t)|  £1.18  R(w)]  , where  s(t)  is  the  resulting  step  response. 

£ max 

Since  |r(w)|  < |F(jw)|  , we  are  conservative  in  using  |F(jw)|  , 

1,11  1 max 


F/0  9/4 


r 
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in  place  of  R(“)max  • We  simply  try  to  use  this  inequality  as  a means 


of  relating  the  bound  max 

t 


xf  (t) 


into  an  w-domain  bound.  This  gives 


Xf  ( ju>) 


1.18 

aw 


We  confine  ourselves  in  fact  to  "reasonably  smooth"  frequency  responses, 
which  are  here  defined  as  those  with  only  negative-real  poles  and  zeros, 
and  experimentally  check  the  above  inequality.  The  results  are  given  in 
Figure  II-l  and  the  table  below,  with  a =3  . The  most  complex  case 
has  4 zeros  and  9 poles.  These  experimental  results  are  in  very  good 
agreement  with  the  above  inequality. 


Case 


1 

2 

3 

4 

5 

6 

7 

8 


max 

t 


xf 

T(t> 


.29 

.30 

.28 

.33 

.30 

.33 

.33 

.33 


If  — (s)  is  not  "smooth",  the  condition 
A 


satisfies  max 
t 


xf 

T(t) 


Xf 

T(i“) 


(—  still 
am 


, but  it  is  too  severe.  For  example,  consider 


a simple  "band-pass”  filter  (twin-tee),  with  complex  zero-pole  pairs. 


1 + 


TT(s)  = 


2 z.  2 

1 . s 

s + — -r 

u>  2 

o U) 

o 


2 z 


1 + 


2 


9 + “T 


with 
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RESULTS  - EXTRAPOLATION  FOR  u = i 


In  all  structures  seen  in  this  study,  there  exists  a "strange" 

combination  of  two  loop  transmissions:  1)  L (ju  ) =— GGKP  (io)  ) 

o J o ttA  1 2 h J o ' 

valid  only  for  the  oscillating  signal  ( w = w ) , 

2b 

2)  L^(jiu)  = — G^G2KP^(jw)  , valid  for  the  slow  signal  components 
( u * uq/3  ) . 

In  the  synthesis  procedures  developed,  the  function  defining 
Xf  go 

— (jw)  , valid  for  ujy  , is  extrapolated  to  w = . This  is  shown 


in  Figure  III-l,  case  1 


(-]T( ju,))  * 
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and 


as  plotted  in  Case  1,  can  be  modified  in  such  a way  that 

lx«, 

[ — ( j uj)  V to  $ o 

3 


ss 

— (jw) 

xfl 

< 

~(^c 

) 

) 

fl(.  x 
— UU)) 


f 2.. 


as  represented  as  Case  2 in  Figure  Ill-i.  In  this  case  the  system  per- 
formance, with  respect  to  the  value  of  , is  improved,  and  the 

Xf2  Xfl 

maximum  value  — — (t)  w — — (t) 

Pi  Pi 


If  such  a solution  is  adopted,  extrapolation  of 


l“> 


is  not 


valid  for  w = w0  » but  a very  complex  system  will  be  achieved,  such 


that  the  amplitude 


f 2 


( jm) 


fl,.  . 


for  co  s • 


This  complex 


system  introduces  extra  delays  that  could  not  be  tolerated,  and 
oscillations  at  frequencies  not  necessarily  that  are  not  allowed, 

since  it  would  mask  the  real  information  about  the  gain  of  the  plant, 
existing  in  the  oscillation. 


These  conclusions  can  be  seen  in  Figures  III-2  and  III-3. 

lx 


Figure  III-2  depicts  a typical  case  where  10  different  o> 


T(jw) 


are 


plotted,  having  the  same  amplitude  for  co  < oij 
final  slopes  are  used.  The  maximum  value 


For  <i>  > , several 


I f i , 
max  — - (t) 

It | ! A 


in  all  cases  is  about  .27  . 


Figure  III-3  shows  the  step  response  for  Cases  la,  lc  and  li. 
tf  relatively  simple  systems  are  used,  as  in  Case  3,  Figure  III-l, 


Xfv 

“■  - ( i <ii  ) 

xfi 

— ~( ju»  ) 

A o 

A O 

■ I 

I 

I 

I 
i 


I 
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